Theory of deuteron stripping. From surface integrals to generalized f?-matrix 

approach. 
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There are two main reasons for absence of the practical theory of stripping to resonance states 
which could be used by experimental groups: numerical problem of the convergence of the DWBA 
matrix element when the full transition operator is included and it is unclear what spectroscopic 
information can be extracted from the analysis of transfer reactions populating the resonance states. 
The purpose of this paper is to address both questions. The theory of the deuteron stripping is 
developed, which is based on the post continuum discretized coupled channels (CDCC) formalism 
going beyond of the DWBA and surface integral formulation of the reaction theory [A. S. Kadyrov 
et al., Ann. Phys. 324, 1516 (2009)]. First, the formalism is developed for the DWBA and 
then extended to the CDCC formalism, which is ultimate goal of this work. The CDCC wave 
function takes into account not only the initial elastic d + A channel but also its coupling to the 
deuteron breakup channel p + n + A missing in the DWBA. Stripping to both bound states and 
resonances are included. The convergence problem for stripping to resonance states is solved in 
the post CDCC formalism. The reaction amplitude is parametrized in terms of the reduced width 
amplitudes (ANCs), inverse level matrix, boundary condition and channel radius, that is the same 
parameters which are used in the conventional i?-matrix method. For stripping to resonance states 
many-level, one and two-channel cases are considered. The theory provides a consistent tool to 
analyze both binary resonant reactions and deuteron stripping in terms of the same parameters. 

PACS numbers: 24.30-v, 25.45.-z, 25.45.Hi, 24.10.-i 
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I. INTRODUCTION 



Production of unstable nuclei close to proton and neu- 
tron drip lines has become possible in recent years, mak- 
ing deuteron stripping reactions (d, p) and {d, n) on these 
nuclei (in inverse kinematics) not only more and more 
feasible as beam intensity increasing but also a unique 
tool to study unstable nuclei and astrophysical (n, 7), 
(p, 7) and (p, a) processes. The deuteron stripping re- 
actions populating resonance states of final nuclei are 
important and most challenging part of reactions on un- 
statble nuclei. If for nucleon transfer reactions popu- 
lating bound states for about fifty years experimental- 
ists used the standard distorted waves Born approxima- 
tion (DWBA), an adequate theory for transfer reactions 
to resonance states yet to be developed. By standard 
DWBA I mean the approach in which the one-step trans- 
fer matrix element is evaluated with incoming and out- 
going distorted waves calculated by fitting the deuteron 
and proton elastic scattering with local optical potentials. 
The transition operator contains finite range effects as 
well as the full complex remnant term. The main idea 
of the DWBA is that the transition matrix element is so 
small that one can use the first order perturbation theory. 
Since the nuclear potential is quite large by itself (~ 100 
MeV), the smallness of the transition operator can be 
fulfilled only if the reaction is peripheral enough, so that 
the non-diagonal matrix element, representing the trans- 
fer reaction amplitude, becomes small. However, since 
the resonance wave function is large in the nuclear inte- 
rior and different channels are coupled in the nuclear in- 
terior, the character of the stripping to resonances can be 



quite different from the stripping to bound states. Nowa- 
days the standard DWBA is gradually being replaced 
by more advanced approaches like continuum dicretized 
coupled channels (CDCC) adiabatic distorted wave 

(ADWA) [1], coupled reaction channels (CRC) and the 
coupled channels in Born approximation (CCBA) avail- 
able in FRESCO code There are two main reasons for 
absence of the practical theory of stripping to resonance 
states which could be used by experimental groups. First 
one is the numerical problem of the convergence of the 
DWBA matrix element when the full transition operator 
is included. However, it is only a technical problem. The 
second pure scientific unsolved problem is what spectro- 
scopic information can be extracted from the analysis of 
transfer reactions populating the resonance states. Be- 
sides, since the standard DWBA is deficient to more ad- 
vanced methods like CDCC or ADWA, a new approach 
should go beyond of the DWBA. 

Majority of theoretical works devoted to the develop- 
ment of the theory of single-nucleon stripping into un- 
bound states of the residual nucleus have been published 
in 1970-s [6U2l|. Great interest in these reactions at that 
time stemmed primarily from the fact that they allow one 
to extract reliable information on the properties of nu- 
clear resonant states by means of the combined analysis 
of the data on stripping and elastic resonant scattering 
of nuclcons from the target nucleus || [TH, 
most of the cited works the theory of stripping into reso- 
nant states was developed within the standard DWBA 
by analogy with usual stripping to bound states. In 
this case the expression for the reaction amplitude ob- 
tained instead of the bound-state wave function for the 
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captured nucleon (form factor) contained a continuum 
wave function which leads to slow convergence of the ra- 
dial integrals or even to their divergence depending on 
the choice of this wave function. In Refs. [a ©, [HI the 
form factor was taken to be a scattering wave function, 
which described the resonant scattering of the nucleon 
from the target nucleus. This wave function was calcu- 
lated using a single-particle potential whose parameters 
were adjusted to give a resonance with the correspond- 
ing properties. The Gamov decaying-state wave function 
and the Weinberg wave function which are regular at 
the origin and purely outgoing at infinity were used in 
Refs. 53 and [14|, respectively. Various methods were 
suggested to calculate radial integrals practically with 
the above-mentioned form factors: (i) the introduction of 
the convergence factor exp(-ar) into the integrand [6|; 
the integral obtained was calculated for various a > 
and then its values were extrapolated numerically to the 
limit of a = 0; (ii) the method of contour integration 
in the complex r-plane (complex scaling) [3; (hi) the 
method based on the correct account of the boundary 
conditions in the three-body scattering problem [ll| ; (iv) 
the Zeldowich-Berggren method [13 of the regularization 
of integrals containing the Gamov function in which the 
convergence factor exp(— ar 2 ) was introduced [13; (v) 
the pseudo-bound-states method [T3|. The methods (ii) 
and (hi) were most convenient for numerical calculations. 
Although the above methods allow one to avoid formal 
difficulties, nevertheless all the methods are rather com- 
plicated because of cumbersome numerical calculations 
and carry on the shortcomings of the standard DWBA 
for stripping to bound states. 

Even if we put aside the technical problem of conver- 
gence of the matrix element for stripping to resonance 
states, there is more important question remains: the 
spectroscopic information which can be extracted from 
analysis of deuteron stripping reactions (and other trans- 
fer reactions) into resonant states. This is really a crucial 
question because the answer determines the reason why 
we measure nuclear reactions. For more than 50 years 
transfer reactions to bound states, and deuteron strip- 
ping in particular, have been used to determine the spec- 
troscopic factors, which measure the weight of the single- 
particle state in the overlap function of the initial and fi- 
nal nuclei. That is why there was always a temptation to 
develop a theory of stripping into resonant states which is 
fully similar to stripping to bound states. For example, in 
13] it was assumed that the spectroscopic factor could be 
extracted from deuteron stripping into resonance states. 
In this case the spectroscopic factor is the ratio of the 
observable and single-particle resonance widths. How- 
ever, the spectroscopic factor is not observable and de- 
pends on the single-particle potential used to calculate 
the single-particle width. In 22] it has been shown that 
spectroscopic factors are not invariant under finite-range 
unitary transformations and, hence, in exact approach 
nuclear reactions cannot be a tool to determine spectro- 
scopic factors. In [22| it was called separation of nuclear 



reactions and spectroscopic factors. However, there is a 
model-independent information, which can be extracted 
from deuteron stripping reactions. I mean the asymp- 
totic normalization coefficients (ANCs), which are the 
amplitudes of the tails of the overlap functions [13 and 
are invariant under finite range unitary transformations. 
The most model-independent definition of the ANC is 
that it determines the residue of the elastic scattering 
S matrix in the pole corresponding to bound, virtual or 
resonance states. For the resonance state the ANC and 
partial resonance widths are related [2J, [25J : 

[c F Ajl ] 2 = (-r/e 2 ^ fc ^«> ^-r xAjl . (i) 

fe a;A(0)j'i 

Here I and j are the orbital and total angular momentum 
of particle x in the resonance state F = (Ax), fi xA is the 
reduced mass of x and A, k x A(o)j i is the real part of the 
resonance relative momentum of x and A, <fij i(k xA ) is the 
non-resonant scattering phase shift, and T xA ji are 

the ANC and partial resonance width in the channel x+A 
with the quantum numbers I and j. Eq. ([!]) stands for 
narrow resonance, i.e. for k xA mji << k x A(o)ju where 
k x A(i)j i is the imaginary part of the resonance momen- 
tum k xA (R)ji = k xA (o)jl ~ik xA (i)ji, which determines 
the location of the resonance pole in the momentum 
plane. Due to relation ([T]), the resonance width is also 
invariant under finite-range unitary transformations and 
can be determined from the experiment. 

Nowadays, it is quite well understood that the ANCs 
can be determined from peripheral transfer reactions, 
see [26l - [3~fl ] and references therein. However, the ANC 
method has been applied only for transfer reactions pop- 
ulating bound states. It is well known that from bi- 
nary resonance scattering and reactions using the con- 
ventional i?-matrix approach one can determine the res- 
onance partial widths, which, as we have underscored are 
related to the ANCs. i?-matrix method is one of the most 
popular tools among the experimental groups worldwide 
because the approach is comparatively simple even for 
many-body, many-channel cases and deals with the for- 
mal partial resonance widths determined from the fit to 
the experimental data. These formal widths can be eas- 
ily related with the observable partial widths. Using the 
i?-matrix approach one can fit simultaneously data for 
all available channels. It allows one to control the consis- 
tency of the obtained physical parameters. The question 
is whether the theory of stripping to resonance states can 
be formulated in terms of the same parameters which are 
used in the i?-matrix analysis of the binary resonance re- 
actions. 

It is the purpose of this paper to address a theory of 
the deuteron stripping, which will solve all the above 
mentioned problems for the deuteron stripping into res- 
onant states. The delivered theory is based on the post 
CDCC formalism going beyond of the DWBA and sur- 
face integral formulation of the reaction theory. The 
CDCC wave function takes into account not only the 
initial elastic d + A channel but also its coupling to the 
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deuteron breakup channel p+n+A missing in the DWBA. 
The convergence problem is also resolved in this formal- 
ism. The reaction amplitude is parametrized in terms of 
the reduced width amplitudes (ANCs), inverse level ma- 
trix, boundary condition and channel radius, that is the 
same parameters which are used in the i?-matrix method. 
Thus the theory provides a consistent tool to analyze 
both binary resonant reactions and deuteron stripping in 
terms of the same parameters. 

The theory is based on the surface-integral formula- 
tion of nuclear reactions and valid for stripping to both 
bound and resonance states. First, just for demonstra- 
tion of the formalism, the transformation of the DWBA 
amplitude for stripping to the bound state is presented. 
The reaction matrix element is split into two parts: inter- 
nal (over the relative coordinate between the transferred 
nucleon and target) and external. The idea of such sep- 
aration is based on the fact that in the post formalism 
the main contribution to the stripping amplitude comes 
from the nuclear exterior while the prior form amplitude 
is dominated by the internal region. It will be shown 
that the dominant external post (internal prior) ampli- 
tude using the Green's theorem can be written as the 
dominant surface integral encircling the internal volume 
plus small addition from the prior external (post internal) 
part. Thus, both post and prior forms lead to the same 
reaction amplitude given by the sum of small internal 
post form, small external prior form and the dominant 
surface integral. The contribution of the post internal 
part can be minimized by a proper choice of the final- 
state optical potential, and the other two amplitudes are 
parameterized in terms of the reduced widths amplitudes 
(ANCs). After that the theory is extended to the CDCC 
formalism. Then the theory is applied for stripping to 
resonance states. First it is developed for the standard 
DWBA and then the post CDCC formalism based on the 
surface integrals is developed. One of the most impor- 
tant results of this paper is that the post CDCC form 
for stripping into resonant states can be written as the 
sum of the small internal (over the coordinate r n A) post 
form and the dominant surface part. The absence of the 
diverging (or poor converging) external part solves the 
problem of convergence of the matrix element for strip- 
ping to resonance state. 

In the developed approach the information about the 
resonance subprocess is contained in the scattering wave 
function of the fragments formed by resonance decay. 
This wave function is written in a standard i?-matrix 
form using its separation into the internal and external 
parts. It allows us to generalize the i?-matrix method 
for binary reactions to stripping reactions. Since the 
deuteron stripping into resonant states is 2 — > 3 parti- 
cles reaction, the excitation of the resonance occurs in 
the subsystem, while the third particle causes the distor- 
tion. The extracted partial resonance widths can be used 
for calculation of the (n, 7) processes. If the cross section 
for (n, 7) resonant capture is available, the simultaneous 
fit to the deuteron stripping and (n, 7) resonance cap- 



ture can be done. The method can be also applied for 
analysis of the Trojan Horse reactions [32]. Concrete cal- 
culations and the application of the theory for deuteron 
stripping and Trojan Horse reactions will be presented 
in the following up papers. In what follows we use the 
system of units in which H = c = 1. We also neglect the 
spins of the particles if not specified otherwise. 

II. SURFACE INTEGRAL FORMULATION FOR 
DEUTERON STRIPPING TO BOUND STATE. 

Before the theory of the deuteron stripping to reso- 
nant states will be outlined I will present a surface in- 
tegral formulation of the theory for stripping populating 
bound states. First, just for demonstration, I consider 
the DWBA and then extend it by including the CDCC 
wave functions. As it has been explained in Introduction, 
the transfer reaction matrix element will be split into two 
parts in the subspace determining the relative motion of 
the transferred nucleon and target: internal and external 
parts. After that replacing the potentials in the transi- 
tion operators by the kinetic energy operators and using 
the Green's theorem the matrix element in terms of the 
surface integral will be obtained. 

A. Stripping to bound state. Post form of DWBA. 

In this section we consider the post form DWBA am- 
plitude, which we split into the internal and external part 
in the subspace over the relative coordinate between the 
transferred n and A. Due to the choice of the transition 
operator in the post form, the internal part turns out to 
be small. The external part, which is parameterized in 
terms of the ANC, will be transformed into the dominant 
surface integral encircling the internal volume and small 
external prior DWBA amplitude. 

We start consideration from the exact reaction ampli- 
tude for the deuteron stripping to bound states 

d + A^p + F, (2) 

where F = (An) is the bound state. The post form of 
the exact reaction amplitude 

M^° st \k pF , k dA ) =< $(~ } | AV pF >, (3) 

where ^ is the exact scattering wave function in the 
initial state with the two-body incident wave d + A, 
&f ^ = X p f V*f i s the channel function in the exit state 
p + F, (fi is the bound-state wave function of nucleus 
i, x\j = x\^\ r ij) is tne distorted wave describing the 
relative motion of particles i and j with the relative mo- 
mentum k^; A V p f — V p a + V pn — UpF is the transition 
operator in the post form, V%j is the microscopic interac- 
tion potential between nuclei i and j, Uij is the optical 
potential between nuclei i and j; r,y is the radius- vector 
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connecting the center of mass of particles i and j. I 
remind that the exact wave function is fully anti- 

symmetrized but the channel wave function $^ ' is not 
antisymmetrized with respect to exchange of the exiting 
proton and nucleons in F. However, the internal wave 
function of F ipp in $^ ^ is fully antisymmetrized. The 
reason why we can drop the antisymmetrization in the 
channel wave function is the presence of the fully anti- 
symmetrized exact wave function in the initial state and 
fully symmetric transition operator what can be seen be- 
low when the transition operator is expressed in terms of 
the kinetic energy operators. 

To obtain the post form of the DWBA from Eq. © 

we replace by the channel wave function = 

,(+) 



'■Pd <Pa XdA m the initial d + A state: 



Then we use approximation 

tfF ~ I A <PA, 



s(+) 



> 



(4) 



(5) 



where 1% (r ra ^) is the overlap function of the bound state 
wave functions of nuclei F and A: 



Ia^ua) = (<Pa\ <Pf)- 



(6) 



Note that the integration in Eq. ^ is taken over all the 
internal coordinates of nucleus A. Then the transition 
operator in Eq. @ takes the form < <pa\& V pF \(fiA > — 
< <Pa\V p a\<Pa > +V pn - U pF . Potential < (pa\V p a\<Pa > 
is replaced by the optical potential U p a and we obtain a 
standard post form of the DWBA amplitude: 



M DW ( post \k pF , k, 



dA 



-< ^\AV pF \^ +) >, 



(7) 



where A V pF — U p a + V pn — U pF . Now we will transform 
this volume integral into the surface one. First, we adopt 
r n A and r pF as Jacobian variables and split the configura- 
tion space over r n A into the internal and external regions, 
while the integral over the second Jacobian variable, r pF , 
is taken over all the coordinate space. Splitting the reac- 
tion amplitude into internal and external amplitudes we 
get 



M DW ^ ost \k pF ,k dA ) 



■pF 



k dA ) 



■M. 



DW(post), 



<-pF 



(8) 



where the internal amplitude M?^^" 



is given by 



M 



DW(post) 
,(-) 



(k P F, Kia) 



-<X p F J lA 7 \&V pF \< PdXd + 2> 

Correspondingly, the external amplitude is given by 



(9) 



M^ ost \k pF ,\, dA ) 



<4- ) 7j|Ay pF |^xS > 



(+) 



r nA >Rn 



(10) 



Here, R n A is the channel radius similar to the one in- 
troduced in the i?-matrix approach, which separates the 
internal and external regions. 

The splitting of the amplitude into the internal a nd 
external parts in the subspace over the Jacobian variable 
r n A is natural and evident. The overlap function I^(r n A) 
is the only object in the reaction amplitude which pro- 
vides spectroscopic and structure information. In the ex- 
ternal region the overlap function has a standard radial 
shape given by the spherical Hankel function (for neu- 
trons) with the amplitude called the ANC (see below). 
To determine the behavior of the overlap function in the 
nuclear interior, which bring one of the main uncertain- 
ties in the analysis of the deuteron stripping, microscopic 
calculations are required 33] . In a standard approach the 
internal part of the overlap function is approximated by 
the single-particle bound state wave function calculated 
in the adopted mean field. The proportionality coeffi- 
cient is the square root of the spectroscopic factor. Due 
to the structure of the transition operator the external 
matrix element M®^^ ^ in the post form is dominant 
compared to a small contribution coming from the inter- 
nal part . This simple observation stems from 
the following. 

In the internal matrix element, r n A < RnA, due ab- 
sorption of the protons inside nucleus F, effective r pn ~ 
fpA ~ r pF > Rf, where R F is the radius of nucleus 
F. For the protons outside of F and neutrons inside or 
on the surface of A each nuclear interaction in the op- 
erator A V pF — U P A + Vpn — U pF is small. Potential 
U pF is arbitrary and often U pF is chosen to compensate 
for U p a so that the transition operator reduces to V pn - 
Since the DWBA is the first order perturbation theory, 
the minimization of the whole transition operator A V pF 
provides smaller higher order terms and, hence, better 
serves the theory. This choice is more preferable in the 
formalism presented here and we adopt U pF , which min- 
imizes A V pF = U p a + V P n — U pF at r n A < R n A making 
the contribution from the internal matrix element small 
compared to the external one. 

In the external matrix element (r n A > RnA), which is 
dominant, the overlap function 1^ can be replaced by 

its asymptotic tail. Although M®^^ ^ can be eas- 
ily calculated for stripping to the bound state, here we 
transform this matrix element into an alternative form, 
which has clear advantage in case of stripping to reso- 
nance states discussed below where convergence becomes 
a main impediment. 

Now we proceed to the transformation of the volume 
integral defining the external matrix element in terms 
of the dominant surface integral encircling the sphere at 
r n A = R n A and a small, due to the structure of the tran- 
sition operator in the prior form (see Eq. (]17jl). external 
volume integral in the prior form. Note that the trans- 
formation is exact within the DWBA formalism. 

To transform the external volume integral to the sur- 
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face one, we rewrite the transition operator as 



A V pF = U p A + Vp, 



u, 



pF 



-u dA ] 

-(U pA - 



- [U p f] 
UdA). (11) 



The bracketed operators are the right- hand- side opera- 
tors in the Schrodinger equations for the initial and final 
channel wave functions in the external region: 



(E - T) <p d Xd + 2 = (V P n + U dA ) <p d Xd + 2 (12) 



and 



(E-T)lZ XpF h = Up F lZxpF 



(-)- 



(13) 



To derive Eq. (|T3|) wc took into account that at r nA > 
RnA I a satisfies the asymptotic Schrodinger equation 
{snA — T u a) I a — 0, where is the binding energy of 
the bound state (ij) and Tij is the kinetic energy oper- 
ator of the relative motion of i and j. These equations 
imply the following connection between the external post 
form DWBA amplitude and the matrix element Mg W 
containing the surface integral: 



M. 



DW(post) 



(k P F, k, 



d A 



M° w (k pF: k 



d A 



-M. 



DW (prior) 



(V> k dA), (14) 



where 



M, 



DW (prior) 



(kpf, k d A) 



=<xU^\AV d A\'p d X d + 2 > 



r„A>RnA 



(15) 



the spherical harmonics and f = r/r is the unit vector. 
We assume that only one value of l n A contributes to ex- 
pansion ([T8|) . If the channel radius is taken larger than 
the range of the nuclear interaction, the radial overlap 
function can be replaced by its asymptotic term, 



*A j nA l nA ( H nA) ~ L 'Aj nA l ■« 



X K n A tl, (l KnA r nA ), (19) 



where h\ \ A (i K n AT n A) is the spherical Hankel function 
of the first order, C A ■ ; is the ANC of the overlap 
function, n n A — \J1 Jhui £nA is the bound state wave 
number. 

It is also useful to introduce the reduced-width ampli- 
tude used in the i?-matrix a ppr oach, which can be ex- 
pressed in terms of the ANC [25[ : 



JnA j nA l nA 



1 A j nA l„ A \ n nA) 



2fi 



ii A 



= \ ^^ A + 1 ^AC F A MlnA h^K nA RnA). (20) 

Correspondingly, the reduced width is 

RnA 



lnAj n A In 
RnA 



2/J>nA 



I 1 A j nA l„ 



,{RnA)Y 
a lnA hfJ A (iK n ARnA)] 2 ■ 

(21) 



and 



M 



g W (k P F, k d A) 



-< Y ( ~ ] I F 



-^\ VdXd + 2 > 



n A >Rn 



(16) 



Here, the transition operator in the prior form A V d A in 
the external region, where the nuclear n — A interaction 
disappears, takes the form 



A V d A = U p a — U d A- 
The overlap function is given by 



(17) 



I F (r nA ) = ^2 < J a M A jnA m jnA \J F M F > 

jnA m jnA m ^ A 

X < J„ M n l nA minA \JnA m jnA > 



x Y, 



l nA mi nA (*nA) J Aj nA U A ( r nA)- 



(18) 



Here, < ji mi jimi\jj,m^ > is the Clebsch-Gordan co- 
efficient, l n A {fni nA ) is the orbital angular momentum 
(its projection) of the relative motion of n and A, j n A 
(rrij nA ) is the total angular momentum (its projection) 
of n in the bound state F — (nA), Ji {Mi) is the spin 
(its projection) of nucleus i; I A i nA j nA { r nA) is the radial 
overlap function, which is a real function [23|, Yj m (f) is 



It is worth mentioning that, due to the presence of the 
channel radius R n A, the reduced width, in contrast to the 
ANC, is mo del- dependent. The dependence on the chan- 
nel radius becomes crucial with increase of the binding 
energy. We are going to use also the boundary condi- 
tion, which is the logarithmic derivative of the overlap 
function at r n A = RnA'- 



B n A = 



1 



d[r n Ahf\ (iK nA r n A)} 



hf' {iK nA RnA) 



dr 



Due to Eq. ([19]). the amplitude M { 



(22) 

DW (prior) 
ext Can De 

parametrized in terms of the ANC. We note that this am- 
plitude is also small. In the external region, r n A > RnA, 
the nuclear n — A interaction can be neglected. Besides 
in this region the overlap function exponentially fades 
away. Also, if the proton absorption is strong in the in- 
ternal region of A, the dominant contribution comes from 
TpA > Ra, where Ra is the radius of nucleus A. If the 
adopted radius channel R n A is larger than the n — A nu- 
clear interaction radius we can neglect n— A nuclear inter- 
action in the external region. In this region each nuclear 
potential U^ A and U^ A and their difference U p a — U d A are 
small. The Coulomb part U^ a ~U^ a k Z a e 2 R d /(2 R\), 
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where R d is the deuteron size and Za e is the charge 
of nucleus A, is also too small compared to the nuclear 
potential. Thus the dominant contribution to the post 
DWBA amplitude M®J (post) , Eq. JH]), and, hence, to 
the total post form DWBA amplitude M DW< - post ^ comes 
from the surface integral Mg W . Here and in what follows 
all the amplitudes with the transition operator 
are assigned the subscript S, which is abbreviation of 
" surface" , because the volume matrix elements of these 
amplitudes can be transformed into the surface ones in 
the subspace over variable r n A while over the second Ja- 
cobian variable r pF we always keep the volume integral. 

Now we express Mg in terms of the surface integral 
over variable r n A and the same technique will be used 
throughout the paper. The kinetic energy operator can 
be written as T — T pF + T n A ■ T pF is a Hermitian oper- 
ator in the subspace spanned by the bra and ket states 
in Eq. (fPo) . It can be proved if we take into account 
that at r p p —> oo the integrand in this equation van- 
ishes exponentially due to the presence of the bound state 
wave function ip d (r pn ) and the overlap function 7^(r n ^). 
Hence, integrating by parts twice the integral over r pF 
we obtain 



<X { JlA\^vF-^pFWdxt2 > 

=<X pF ] Ia\^ p f-T* p i 



r nA >R„ 



r„A>RnA 



(23) 



Then Mg w reduces to 



M 



--<X ( pF T A I TnA - ~^nA \ (fid X d + A > 



r nA >R n 



(24) 



We apply now Green's theorem to transform the vol- 
ume integral into the surface one, which encircles the 
inner volume over the coordinate r: 



dr/(r) 



S(r) 



r<R. 



= -i- j dS [g(v) V r /(r) - /(r) V r g(r)\ 



r=R 



1 



= R 2 I dn T 

2 /i 



Or 



dr 



r=R 



(25) 



Here, d S = R 2 d ft f , where is the solid angle. Note 
that the unit vector f is the normal vector to the sphere 
directed outside of the restricted by the surface volume. 
The integration in Eq. (|2"4")l over v n A is taken over the 
external volume restricted by two spherical surfaces: the 
inner surface with the radius R n A and the external sur- 



face with the radius i?_ 



A I 



DW 



oo, that is 



(k pF , k dA ) = -M^ nA (k pF , k dA ) 
+M^(k pF , k dA ). 



<dA). (26) 

The first term in this equation is the surface integral 
encircling the inner surface of the external volume at 
r n A = R n A while the second term is the surface integral 
taken at r n A = R n A ~ * 00 ■ A negative sign in front of 
the first term appears because the normal to the surface 
is directed inward to the center of the volume, i.e. op- 
posite to the normal to the external surface (at infinitely 
large radius). The second term vanishes because of the 
presence of the overlap function JjT, which decreases ex- 
ponentially at r n A — >• oo. Then for Mg W we get 



M° w (k pF , k, 



<IA 



(k pF , k 



d A 



2[i 



ti A 



dr pF X- + kL( r P f ) 



do. 



Vd(r pn ) X { kJ A (*dA) 



d K(v nA )Y 

dr nA 



u f, ^ 9 Vd( r m)x^ A ){*dA)- 

[J-A \ V nA)\ 







TnA 



- r rlA =RnA 

(27) 



Here we took into account that xL ( r ) = X-k( r )- 
voking Eqs. (fl8|) and (p~9|) we can rewrite Mg W in the 
form explicitly showing parametrization in terms of the 



reduced width amplitude (ANC) and boundray condi- 
tion, the quantities used in the i?-matrix approach: 
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M° w (k pF , k dA ) 



1 K nA R n A hfJ A (i K nA RnA) ^ < Ja Ma i nA I j f M f > 



jnA ™.j nA m i nA M n 



x < J n M n l nA m lnA \j nA m jnA >< JpM p J n M n \J d M d > C A • K 



<Pd(r pn ) xL A ( r dA) (B nA - 1) - RnA 



dr nJ 



r nA =RnA 

(28) 



' RnA 
2^nA 



^2 < J A M A j nA m jnA \J F M F >< J n M n l nA mi nA \j nA mjnA ' > 



jnA m jnA mi nA M„ 



x < JpM p J n M n \J d M d > ^ n Aj nA i nA J dr pF x { -k pF (r P F) J dfl rnA 



* m l n A 



JnA) 



dtpd{r pn )x^ A ){*dA) 



tPd{r pn ) Xk dA {*dA) (B nA - 1) - RnA 



dr nJ 



r nA = R. n 



(29) 



by 



Finally, the total post form DWBA amplitude is given B. Prior form of DWBA. Stripping to bound state. 



M DW ^(k pF , k dA ) = M™ {post \KF, k rf A) 



+M, 



DW (prior) , 
ext 



(kp F , k dA ) + M° w (k pF , k dA ). (30) 



M. 



DW{post) 

DW( P nor) wg rewrite Eq. (|30[) in a different form: 



Taking into account that Mg w = M" x " t 



M DW( P ost) {kp ^ ^ = M DW( P ost) {kp ^ k ^ 



M^ p " or \k pF ,k dA ) 



+ 



M. 



DW( P ost) 
ext 



(k pF , k dA )-M ( 



DW (prior) 



(31) 



Thus, the main result of this section is that the post 
form of the DWBA amplitude can be written as the sum 
of the peripheral parts, y[®w{ P nor) j^dw ^ anc j sma n 

internal term M^^^ 1 " 3 ^. The peripheral part itself con- 
sists of the dominant surface amplitude Mg W and small 

external prior form M^^^ prlor \ The peripheral part is 
parametrized in terms of the ANC (reduced width ampli- 
tude) , channel radius R nA and the logarithmic boundary 
condition, that is in terms of the parameters used in the 
i?-matrix fitting. The model dependence of these two pe- 
ripheral amplitudes is caused by the ambiguity of the op- 
tical potentials and channel radius R nA . The strongest 
model dependence comes from M^^^ post \ because, in 
addition to the ambiguity of the optical potentials, to 
calculate it one needs to know the behavior of the overlap 
function in the internal region. For peripheral reactions 



contribution of M, 



DW( P ost) 



can be neglected. 



In subsection A the post form of the DWBA amplitude 
has been considered. However, all the results hold also 
for the prior form 

M DW{p " or) (k pF , k dA ) = < Xp -J /H AV dA \ wxfl > 



fOWirrior)^^ kdA) + M DW( pn or ){kpFi ^ 



(32) 



where 



M. 



DW (prior) 



(k p _p, k d A) 



= <X { i F ] I 1 A\^V d A\VdX d + A ) > 



(+) 



r nA <RnA 



(33) 



and 



M. 



DW (prior) , 



=< x y i F A \^v dA Wdx { Zl> 



r nA >Rr, 



■',.,-/ (k,,; • k dA 

<.pF T 

with the transition operator 

AF dj4 = U p A + VnA ~ UdA 



(34) 



(35) 



The n — A interaction potential V nA —< ip A \V nA \ip A > is 
the mean field real potential supporting the bound state 
(nA). The splitting of the amplitude into the internal 
and external terms in the subspace over the coordinate 
r n A helps us to further transform the prior DWBA ampli- 
tude. Due to the structure of the transition operator the 
external matrix element M^^^^ in the prior form is 
small (see the discussion in subsection lll A[) and the main 
contribution in the prior form comes from the internal 
part '( prlor ) _ Since the internal part is given by the 
volume integral, its calculation requires the knowledge of 
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the overlap function in the internal region. The model 
dependence of the overlap function in the nuclear inte- 
rior (r n A < RnA) brings one of the main problems and 
main uncertainty in the calculation of the internal matrix 
element. However, using the surface integral we can re- 
distribute the internal contribution in terms of dominant 
the surface term (over variable r n ^) plus small internal 
part written in terms of the volume integral in the post 
form. With reasonable choice of the channel radius R n A 
the contribution from the internal volume integral in the 
post form can be significantly decreased compared the 
surface matrix element. The latter can be expressed in 
terms of the i?-matrix parameters - the observable re- 
duced width amplitude (ANC), boundary condition and 
channel radius. To transform j^P^\P rlor > j n ^ ^he sur- 
face integral in the subspace over variable r n A we rewrite 
the transition operator in the internal region as 

A VdA = U pA + V nA ~ U dA 

= [V nA + U pF ] + (UpA + V pn - U pF ) - [Vp„ + UdA]. 

(36) 

The bracketed transition operators are the potential op- 
erators in the Schrodinger equations for the initial and 
final channel wave functions. Hence, for the internal prior 
form of the DWBA we obtain 

<r (p ™ r) (v,k dA ) 

= MZT {p ° st) (KFi Ku) + M^ pHor \kp F , k dA ), 

(37) 

where 

M° w (k pF , k dA ) = - < *y /ji ^ - f i % Xd + 2 > 

= -< XpF /£ I ^nA - T*nA\ Vd X d + 2 > 

= -M^Jk pF ,k dA ). (38) 

Note that here is the surface integral encircling 

the border of the internal volume at r n A = R nA with the 
normal directed outward. Thus we have demonstrated, 
what should be expected from the very beginning, that 
M DW( P nor) = M DW{ P ost)_ Hence all the equations ob- 
tained in the previous subsection III Al are also valid in 
the prior formalism. 

It is worth mentioning that in the post formalism, in 
contrast to the prior one, we have obtained two surface 
integrals (in the subspace over r n A) with the radii r n A — 
R n A and r n A = R n A ~^ 00 an d then proved that the 
second integral is zero. From the equality of the post 
and prior DWBA amplitudes we could conclude that the 
surface matrix element over infinitely large sphere r n A — 
R nA oo, which appears only in the post formalism, 
vanishes. 

There is another interesting point to discuss which ex- 
plains the advantage of the above outlined formulation 
of the stripping. As we have discussed, due to differ- 
ent structure of the transition operators in the post and 



prior forms, the main contribution to the post (prior) 
form comes from the external (internal) part (in the sub- 
space over variable r n A)- Since both forms give identical 
amplitudes, that is, describe the same reaction mecha- 
nism and the same physics, such redistribution of the 
main contribution is possible only if the main contribu- 
tion to each form comes from the border between external 
and internal parts. In the post (prior) form this border 
attributed to the external (internal) form and can be ex- 
pressed in term of the surface integral. Let us rewrite 
equality M DW ^ nor ^ = M DW ^ ost ^> in the following form: 

M°J^"° r \k pF , k dA ) + M^ rior \k pF , k dA ) 

= M°J^{k pF , kcu) + M°J^(k pF , k«u). 

(39) 

In this form the dominant terms are ^/[P^ r (.P rior > anc [ 
M^Y^ 05 ^ while the rest two terms, M^t vt " prlor ^ and 
M i ^^ pos *' ) are smaller. From Eq. (j39| we get 

M°Y ipost) (kp F , k dA ) - M™ (p "° r) (k pF , k dA ) 
= M™^ ior \kp F , k dA ) - M^ {P ° St) (KF, ^a) 
= Mg w (k pF , k dA ) = ~M^ a (k pF , k dA ). (40) 

Thus the difference between the post and prior external 
amplitudes (or the prior and post internal ones) is the 
surface integral in the subspace over r n A- 

There is one more point left to discuss. When deriving 
the post form of the DWBA amplitude from Eq. ((4]) we 
used approximation ip F « I a 'PA neglecting the contri- 
bution from the channels n + A n , n > 0, where A n is the 
excited state of A. However, I will show now that the 
surface integral formulation doesn't require this approxi- 
mation. To this end let us split M(P° st '> into the internal 
and external parts in the subspace over variable r n A- In 
the internal part we use a standard DWBA approxima- 
tion (p F w /J[ ip a to arrive to the standard internal post 
DWBA amplitude. In the external part we rewrite the 
transition operator as 

A Vp F = V p A + Vp n - UpF 

= -[V A + U pF ] + [Vpn + V A + U dA ] + (VpA - U dA ) . 

(41) 

The bracketed operators are the right-hand side opera- 
tors of the Schrodinger equations 

(E - T) $| +) = (Vpn + V A + U dA ) *| +) (42) 

and 

(E - T) <S> { f } * = (V A + U pF ) $^ } *. (43) 
Hence, the external part of M^ post ^ reduces to 
M^ st \k pF ,k dA ) 

= M s{ext) (k pF , k dA ) + M^W, k^), (44) 
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where 



r nA >Rr, 



and 



Ms(ext)O i pF, k 



dA 



(45) 



r n A>R n A 

(46) 



In the matrix element M, 



(prior) 



we can use a standard 



DWBA approximation cpp rs I^fA which leads to the 
standard external prior DWBA amplitude. The matrix 
element Mgr ext \ can be rewritten as 



= < X p ~F <Pf\ T nA ~ ~^ nA \fd H> A X d + A > 
=<X ( Jll\^nA-^nAW d X { d + 2 > 



r n A>R n A 

r nA >R n j 
r nA >RnA 



MJT (k pF , k dA ), 



We took into account that < $ 



-L A 



(+) 



/ 



pF + T / 



> = 0, where Ta is the internal motion kinetic 



energy operator of nucleus A, and T n A <-Pa = <PA T n A ■ 
Thus Mgiext) can be transformed to the surface inte- 
gral over variable r n A encircling the inner volume with 
the radius r n A = RnA without invoking approximation 
(pp w /JT if a- It means that, when deriving the post form 
of the DWBA amplitude, the approximation ipp rs I A tp A 

is required only to obtain two small terms, M®^^ post ^ 
and M^ x Y^ pnor \ but not the dominant surface term 
—Mg^ . In this sense the surface integral formalism 
is an improvement of the DWBA. 



Deuteron stripping to bound states. Post 
CDCC formalism 



In the previous sections we succeeded to parametrize 
the DWBA amplitude in terms of the ANC except for a 
small term, M®^^ post \ The most serious shortcoming 
of the DWBA is that it neglects the coupling to open 
reaction and breakup channels. This coupling can be 
taken into account if an exact wave function in the initial 
or final states is used. However, the exact wave functions 
are not yet available (if they would be available in the 
whole configuration space, we don't need to calculate the 
matrix element because the asymptotic terms of the exact 
wave functions provide the reaction amplitudes in all the 
open channels). Here we use the CDCC formalism, which 
takes into account the elastic d + A and the deuteron 
breakup channel p + n + A in the initial state. 



In this subsection the surface integral formulation of 
the reaction theory will be applied to the post form of the 
CDCC amplitude for deuteron stripping to bound states. 
It will allow us to parametrize the stripping amplitude in 
the CDCC approach in terms of the i?-matrix parameters 
- the reduced width amplitude, boundary condition and 
the channel radius. To obtain the CDCC wave function 
describing the initial state of the stripping reaction, first 
the exact initial scattering wave function ^ is replaced 

by the three-body wave function ty? B ( + ' , which takes into 
account the coupling of the initial channel d + A and the 
deuteron breakup channel p + n + A [J-Q and satisfies 
the Schrodinger equation (in the three-body p + n + A 
model space) 



{E — T — U p a — U n 



3B(H 







(48) 



with the outgoing waves in the elastic channel d + A and 
the breakup channel p + n + A. A general solution of 
this equation with the d + A incident wave has outgo- 
ing waves in the elastic, breakup and two rearrangement 
channels, n+(p A) and p+ (n A) . To damp rearrangement 
channels in the asymptotic behavior of the wave function 
the optical potentials U p a and U„a with strong 



3B(+) 



(47) imaginary terms can be used [34|. \1/ 



3B(H 



\r d A, Tpn) = (fi d (r pn ) X ( Cl( r dA) 



is given by 



d PpnV'p^( r pn)Xp( + )(p p 



(49) 



Here, tpd(fpn) is the deuteron bound state wave function, 
V'ptn ( r pn) the p — n scattering wave function with the 
relative momentum p pn , X^] A {^dA) and Xp( Pj)n ) { v <ia) 
are the expansion coefficients, EdA — £ pn — P 2 /(2 ^dA) + 

P 2 pn/( 2 fJ-pn)- 

In practical application the wave function ty^ 3 ^ is 
replaced by the CDCC wave function, which is a solution 
of the projected Schrodinger equation 

(E — T — Ufr ] - U { n P A pn) - V pn ) ttf DCC{+) = 0. (50) 

Here, U^ pn) = P pn U pA Ppn, and 



P„ 



x y; 



E E 

Z p „=0 m lpn =-l p „ 



d ^r p „ Y[ pn m lpn (fpn) 



l l P ,^n lpn \ L pn) ( 51 ) 

is the projection operator, which truncates the number of 
the spherical harmonics Y\ m , n {^pn) in the coordinate 
r pn . Application of this operator to the three-body wave 
function suppresses the rearrangement channels in the 
asymptotic wave function. The CDCC wave function is 
taken in the form 



CDCC(+) 



(rpn,T d A) = Ppn £ ^(^J xf K+) (^), 



n=0 



(52) 
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where ippn( r P n) ~ Vdi^pn) is the deuteron bound state 
wave function, i[>pn (T pn ), n > 1, is the n-th discretized 
continuum state of the p — n pair obtained by averaging 
continuous breakup states in the n-th bin, Xi \ r dA) 
are the functions, which describe the relative motion of 
the center-of-mass of the p — n pair in the n-th state and 
A. Note that x\° ( r dA) asymptotically behaves as the 
incident Coulomb distorted d — A plane wave plus outgo- 
ing scattered wave, while x l ( " )(+) (i - dA) for n > asymp- 
totically do not contain any plane wave having only the 
outgoing scattered wave. 

To derive the post form of the CDCC amplitude from 
the exact one, first we replace the initial exact scattering 
wave function <fj +) by <£u#f B(+) . Note that $fB(+) ig 
the three-body model {p + n + A) wave function which 
treats nucleus A as a constituent particle leaving its inter- 
nal degrees of freedom intact. That is why the wave func- 
tion \f^ + ^ is approximated by the product of the bound 
state wave function if a and ^ B< > + \ Correspondingly, 
the transition operator AV p f — V v a + V pn — U pF is re- 
placed by A V p f = U p a + V pn — UpF ■ This replacement 
of the microscopic potential V p a in the exact post form 
amplitude by U p a is evident because the p— A interaction 
potential in the transition operator should be the same 
as the one in the Schrodinger equation for the initial scat- 
tering wave function ^\ B ^ + \ Potential V pn remains the 
same when we approximate the initial exact scattering 
wave function ny the three-body one. The final state 
optical potential U p p is arbitrary and we discuss the op- 
timal choice of this potential later on. These approxima- 
tions lead to the expression for the post form stripping 
amplitude in the three-body model in the initial state: 



M 3B ^ ost \k pF , k dA ) 

=<*y <p F \Av PF \<p A *i B (+) > 



< X pF 1 A 



I F A \AV pF \ i S> 



> 



the post form of the CDCC amplitude takes the form: 

M CDCC( P ost ){KFj ^ 

=<X ( p- F ] ll\AV P pT\^ DCC{+) >- (54) 

Now we split M CDCC into the internal and external parts 
in the subspace r n A- 



M CDCC ^ ost \kp F ,k dA ) = M i 



_ A *CDCC(post) 



+ M, 



CDCC(post) 



(kpF) k^). 



(55) 



The internal amplitude M^ GG ^ post ^ is given by 



M 



CDCC(post) 



in I 



(k P F, k 



dA 



=<x p f^|a<f^ cc(+) > 



(56) 



Correspondingly, the external amplitude is 



M. 



CDCC{post) 



(kpF, Kia) 



ext 

=<xtf/j[|AT^|*f DOO W> 



r„A>R n 



(57) 



I remind that the integral over the second Jacobian vari- 
able, r p F, is taken over all the coordinate space. Sim- 
ilarly to the DWBA case, the internal part is small if 
the channel radius R„a is not too large. Due to the 
strong absorption of the proton inside A, which is con- 
trolled by the imaginary part of the optical potential 
U the effective distances are t v a > Ra- Besides, 
in the internal region, r n A < RnA, and large r p A, where 

tP-o: 



TpA ~ r V n = \vp A - r nA \, U p l n + V P n can be well ap- 
proximated by a properly chosen optical potential U pF 
p 

minimizing A V pF and the internal matrix element. The 
next step is to transform the external matrix element to 
the surface one. To this end we rewrite the transition 
(53) operator in the form 



Thus, even if we treat the d + A collision in the initial 
channel in the three-body approach, the final state con- 
tains the overlap function, which is essentually many- 
body object. Eq. (|53[) is impractical to use because it 
requires the knowledge of the three-body wave function 
Eq. (|49p . which contains unknown expansion co- 
efficients Xk dA (rdA) and XP(p pn )( T dA)- In practical appli- 
cations the is approximated by the CDCC wave 

function ^ CDCCi * + \ which requires the knowledge of the 

finite number of the expansion coefficients. They can be 

found from the coupled equations. Correspondingly, the 

transition operator AV pF = U p a + V pn — U pF in Eq. (|53")l 
p p 

is replaced by AV pF = U + V pn — U pF . Note that 
only the potential U p a{y p a), where r p A = r d A + 1/2 r pn 
is affected by the projector P pn - Then the expression for 



u pF = [-U pF ] + [u^ n + v, 



v pn\ 



(58) 



The bracketed operators in (f58|) are the right-hand-side 
potential operators in the Schrodinger equations in the 
external region r n A > RnA, where the nuclear n — A 
interaction vanishes: 

(E - T) tff DCC{+) = {U p l n + V^) tff DCC{+) (59) 

and 

{E-T)j$*I%* = U pF £->lZ*. (60) 
Note that the second equation follows from 

(s nA - T nA ) I A =< <PA\V n A\ Vf > ■ (61) 
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In the external region, r n A > R n A, the source term 
on the right-hand-side disappears and Eq. becomes 
evident. Taking into account Eqs (|59l> and (f60|) we get 



M, 



CDCC(post) 



(k pF , k dA ) = M s 
<X pF 'lZ\T-T\V 



_ ,rCDCC(post) 



(kpF, k d A) 



> 



r„A>R„ 



(62) 



where T = T pF + T n A- Here, as in the previous section, 
for the surface integral we use the subscript " S" . Since 
the CDCC wave function doesn't propagate into the final 
state (its asymptotic terms have only elastic and breakup 
terms) the operator T pF is Hermitian, i.e. 



<x { Ji. 



pF 



pF 



pF 
2* 



pF\ 
T CDCC(+) 



r nA >R. 



r nA >R nA 
= 0. 

(63) 



It can be also shown explicitly taking into account that 
the volume integral over r p p can be transformed into 
the surface integral over the sphere with the radius 
r p p = R pF — > oo. Since the overlap function decays 
exponentially at r n A — > oo, the integration over r n A is 
limited. Hence, at r pF — > oo using Eqs (|B6[) we get 



l/r pF , see Eq. (|B16|) . Hence the surface integral vanishes 
at R nF -> oo as R^p/Rpp -> 0. 



IpF 

Then M 



CDCC(post) 



takes the form 



CDCC(post) 



M 



= <X P F ] lA^nA-^nA^ 



CDCC( + ) 



> 



r nA >R„ 



-M. 



CDCCipost) 



(k pF , k dA )+M^ CC{post) (k pF , k dA ). 

(64) 



Thus, the volume integral at r n A > R n A in the ma- 



trix element M, 



CDCC( P ost) 



can be written as the sum of 



two surface integrals encircling the external volume, the 
sphere with the radius r n A = R n A and the sphere with 
r n A — R n A ~~ ^ 00 • Note that the integral over r p p is 
taken over all the coordinate space. Evidently that the 
integral over the infinitely large sphere vanishes because 
the overlap function I~[ exponentially decreases. Hence, 



M, 



CDCC{post) 



(k pF , k dA ) = -M 



CDCC{ P ost) 



{k p Pi k, 



dA 



(65) 



that r d A ~ r p p — > oo and r pn ~ r pF — > oo. The first 
term of the CDCC wave function decays exponentially 
at t p f —¥ oo because of the presence of the deuteron 
bound state wave function. The terms with n > 1 decay 
as l/Vpp [35[. The distorted wave X p ~p*(Ppf) decays as j^CDCCipost) 



The negative sign in front of the inner surface integral 
appears because the normal vector to the inner surface 
is directed to the center, i.e. opposite to the direction of 
the normal to the external surface at r n A = R n A 00 • 
Now we can use equations from subsection A replacing 
the initial channel wave function by the CDCC one. For 



we get 



r?2 

n nA 
2 UnA 



R 



2/J.nA 



M 



CDCC{post) 



(kpF, k dA ) = -Mg®° c ( post \k pF , k dA ) 

[lI(rnA)}*&r nA -^r nA )*? DCC{+ \rp F , r nA ) 

CC + \r pF , r n A)~ 



dr pF£ + k ' M / dS!,„„ [*f DCC<+) (r,f> 



dr n 



r nA >R„ 



r nA =R nA 

(66) 



Natural Jacobian variables for ^ DCC ^ +S > are VdA anc j 
r„„, but here we use another set of Jacobian variables, 



r pF and r n A- Taking into account Eq. (JT5J) and (|19p wc 

get 
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M 



CDCC(post) 
S 



(k pF , k dA ) = -M Sr 



CDCC(post) 



k p F, k dA ) 



Rn 



2 fJ>nJ 



E 



< J a M A j nA m jnA \J F M F > 



3nA "lj„A mi „A 



x < J n M„ Z„a m/ nA | j„a mj nj4 > 7„a j„ a f nJl / d r pF X^ pf , (r pF ) / d fi r „ A Y", 



CDCC(+) 



{v pF , Y nA ) (B nA - 1) - RnA 



d^ DCC ^(r pF ,r nA ) 
dr nA 



rnA=Rn/ 



(67) 



Note that the CDCC wave function itself also depends on 
quantum numbers of p — n and d — A subsystems, which 
we don't specify here. It will be done in the following up 
paper where concrete calculations will be presented. 

Thus we have obtained a remarkable result: the post 
form of the CDCC amplitude, in contrast to the DWBA 
one, is given by the sum of only two terms: 



M CDCC( P os t){kp ^ ^ = M CDCC {p os t){K ^ ^ 

0. 

(68) 



. M CDCC (post){kpF)kdA) 



where the first term, which is the internal post form of the 
CDCC amplitude, can be minimized by a proper choice 
of U pF and the channel radius R nA , while the second 
term, which is dominant, represents the surface integral 
with the radius R nA , which encircles the internal volume 
in the subspace over the coordinate r nA . If the channel 
radius is larger than the n — A nuclear interaction radius 
the second term is parametrized in terms of the reduced 
width amplitude (ANC of the projection of the bound 
state wave function of F on the two-body state n+A) and 



the boundary condition at r nA = R nA . 
is small enough, 



If M 



CDCC(post) 
ini 



M CDCC ^ ost \k pF , k dA ) « -AT 



CDCC(post) 



(69) 



Thus we succeeded to parametrize the post form of the 
CDCC amplitude in terms of the i?-matrix parameters. 
Eq. (|68[) and parametrization of the surface term of the 
post CDCC amplitude in terms of the i?-matrix param- 
eters, Eq. (1671) . are one of the main results of this paper. 



r CDCC(post) 



can be min- 



Although it is assumed that M int 

imized so that the second term in Eq. becomes 

dominant, I would like to present a different form for 

Mint >CCt " P ° St \] i p F , kd A ), which leads to a different form 

for the whole amplitude M CDCC< -P ost ^{k pF , k dA ). To 

p 

this end, let us rewrite the transition operator A V pF " 



in M 



CDCC(post) 



k pF , k dA ) as 

UpF 



PpA + U nA + Vpn] ~ [VnA + U pF ] + V nA 



u 



(70) 



Here, V 71 a is the mean field potential supporting the 
bound state (nA) while U n v A is the projected optical 
potential describing the n — A interaction in the initial 
state of the reaction and entering the Schrodinger equa- 
tion for the projected CDCC wave function in the initial 
state. The bracketed potential operators are the right- 
hand-side operators of the Schrodinger equations in the 
internal region, r nA < R nA , 



and 



CDCC{+) _ (TT P pn 
— \ U pA 



(71) 



(E T) XpF } * I F A * = (V nA + U p p) Xp t I a * ■ (72) 



Replacing the bracketed potential operators [U pA n + 
UnT + V pn ] and [V nA + Upp] by E - f and E - T, 



correspondingly, we get for M 



CDCC(post) 



a new form: 



M anaa {post)(kpFjkdA) 



M 



CDCC{post) 



(k pF , k dA ) 



+ Mg° cc ^(k pF ,k dA ), 



M f„^«(k pF ,k dA ) 



(73) 



V pF 



Ppn^CDCC(+) 



> 



r nA <R n A 



AV^ A =V nA -U^ A 



(74) 



(75) 



Then the total post form of the CDCC amplitude can be 
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written as 

M CDCC( P ost) {KF) ^ = M CDCC {p ost) {K ^ ^ 



dA) 



rt nA 



CDCC(post) 



(k pF , k d A) 

M CDCC( P ost Hkp ^ ^ = M C S CC( pOS t) (kpF) kdA) 



r nA <R n A 



(76) 



Thus, we obtained another important result. The CDCC 
amplitude in the post form is equal to the inner volume 
integral over variable r n A with the transition operator 



u 



nA ■ 



This transition operator is the difference 



between the bound state potential V n A supporting the 
final bound state (n A) and the projected optical poten- 
tial describing the n — A interaction in the initial state. 
It is worth mentioning that Eqs (|68p and (f?6"|) are ex- 



act within the CDCC approach. If j^^ 00 ^ ^ [ s sma vj 
enough, then 



M CDCC( P o st){k 



■pF, K dA) 



M 



CDCC(post) 



(77) 



However, I prefer Eq. 



rather than (1751 . To calcu- 



late Mam CC ^ P ° st ^ one needs to know the overlap func- 
tion in the internal region, where the overlap function is 
model-dependent and requires microscopic calculations. 
In contrast, in Eq. (|68l) the dominant part is the surface 
integral, which is parametrized in terms of the reduced 
width amplitude (ANC). The model dependence of the 
surface part is related with the ambiguity of the optical 
potentials and the value of the cut-off orbital angular mo- 
mentum in the p — n subsystem in the CDCC approach. 
Comparison with experiment allows one to extract the re- 
duced width amplitude. The model-dependent internal 
part in Eq. (155)) is small. Eqs (f6"5f and (I75|) is prelude 
to the theory of the stripping to resonance, where the 
convergence problem of the external part is one of the 
main issue. As we have demonstrated in the post CDCC 
formalism the external part doesn't appear at all. It re- 
solves the convergence problem related with the external 
part. 



D. Deuteron stripping to bound states. Prior 
CDCC formalism 



A priori, the amplitudes in the post and prior forms of 
the CDCC formalism are not equal. That is why the ob- 
tained equations using the surface integrals are expected 
to be different in both formalisms. The prior form of the 
CDCC stripping amplitude is 



-< * 



CDCC{- 



(78) 



where 



A v dA - U pA 



+ V n A — UdA- 



(79) 



The projected CDCC wave function in the final state is 
a solution of the three-body Schrodinger equation 



{E-T-U^ 



V, 



nA 



Vp^)^ 000 ^* =0. (80) 



Here, 



PnA = 



E E 

1„a=0 mi =—l nJ 



r„A J lnA mi 



(fnA) 



X Y, 



InA mi 



(Ka)- 



(81) 



is the projection operator, which truncates the number of 
the spherical harmonics Yi nAmi (r n ^J in the coordinate 
r n A- 

Now, as usually, we split the amplitude M CDCC( - prior "> 
into the internal and external parts in the subspace over 
variable r n A- 



(V, k dA ). M ODaa <* rt °*Q hF , k dA ) 



= M 



CDCC(prior) 
int 



(k pF} k 



dA) 



M. 



CDCC(prior) 



(k pF , k d A), 

(82) 



where 

M CDCC ^° r \k pFl k dA ) 



=< f 



CDCC(~ 
f 



\UX+VnA-U d AWdX d + 2 > 



r n A<R n A 

(83) 



and 



M, 



CDCC(prior) 



{k p F, k d A) 



-< * 



CDCC(~) 
f 



\U^a A -UdAWdX^ > 



( + ) 



r„A>RnA 



(84) 



The external part of the prior amplitude (see discussion 
in subsection III Bp , due to the structure of the transition 
operator, is small and the dominant contribution comes 
from the internal amplitude. We will rewrite this ampli- 
tude singling out the surface integral over variable r n A- 
To do it we rewrite the transition operator 



pA 



Av p d r = u p p r 

+ v nA + V%*]- 



+ VnA — UdA 
[V pn +UdA 



]-+(Vpn-V p ^) 

(85) 



The bracketed operators are the right-hand-side opera- 
tors of the Schrodinger equations 



[E - T) * 
and 



CDCC(-)* 
/ 



pA 



V P nA ^iuCDCC(-)* 
' pn J 



(86) 



(E - T) ifd x ( d + 2 = [Vpn + U dA ) Vd XdA ■ (87) 
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Taking into account these equations we can rewrite 

CDCC(prior) , 
int 



MT„Y^^ yl " lu ' '(k pF , k dA ) in the form 



M, 



CDCC{prior) 



(k pF , k dA ) = M 

+Ae cc(pnor) (v,^), 

where 



CDCC{prior) 



(k pF , k, 



dA 



M CDCC( P rior) {KF) ^ 



-< 



CDCC(-) 



V. 



pn 



VdXdA > 



r nA <R nA 



(89) 



and 



CDCC(pr 



= - < * 



(k pF , k. d A 

CDCC( 

f 



T -T 



(+) 

Vd XdA 



> 



\<R n 



(90) 



Here, the kinetic energy operator T — T pF + T n A- In 
j^jCDCC(pnor) ^ e vo i ume integral over r pF can be trans- 
formed into the surface one taken over the sphere with the 
infinitely large radius: r pF = R p p — > oo. For r nA < R n A, 
due to the presence of the deuteron bound state wave 
function, the integrand goes to zero exponentially, that 
is this surface integral vanishes. Hence, only the sur- 
face integral encircling the inner volume with the radius 

TnA = RnA- 



M 



CDCC(prior) 



(k pF , k t 



dA 



< f 



CDCC(- 
f 



T nA - T nA | (p d x ( d + 2 > 



r nA <Rr, 



M CDCC { pos t){KF ^ dA) _ 



M, 



is 



given by Eq. 



(91) 

(EH). 

(k p F, k dA ) is an auxiliary internal part, 
which is small because at r n A < R n A and r p p > Rp due 
to the proton absorption in the nuclear interior, p — n 
nuclear interaction is significantly depleted, and so the 



CDCC(post) 

Sl *nA 

CDCC{prior) 



difference V, 



y pn 



Then 



M CDCC ^" or \k pF , k dA ) 

= MZ D x CC ^° r Hk pF , k dA ) - Ml DCC ^° st \k p p, k dA ) 

+ M^ CC{prlor \k pFl k d A), (92) 

Thus the total prior form CDCC amplitude consists of 
three terms, small auxiliary internal part, small external 
prior form and the dominant surface term. We can see 
that post and prior CDCC formalisms are not equivalent. 
In the approach used in the paper the configuration space 
over variable r n A was split into the internal and external 
parts. As it has been discussed in Introduction, such a 
splitting is natural because the main object of interest in 
the analysis of deuteron stripping is the overlap function 
I A of the bound states wave functions of the target A 
and final nucleus F. Its external part (r n A > Rua) is 



parametrized in terms of the observable ANC while the 
internal part is mo del- dependent. 

In the post formalism the external part is domi- 
nant. Invoking the post CDCC formalism allows us to 
rewrite the external CDCC matrix element in the form 
of the surface integral over variable r n A, which can be 
parametrized in terms of the parameters used in the R- 
matrix method for binary reactions, while the model- 
dependent internal part gives small contribution. Thus 
the volume part of the matrix element over variable r n A 
is transformed to the surface integral. For transfer to 
bound states such a transformation doesn't bring any sig- 
nificant advantages because the volume matrix element 
converges. However, for stripping to resonance states 
(see subsection MI C|) this transformation provides a de- 
cisive benefit because it solves the convergence problem 
of the matrix element. Here, the transformation of the 
post CDCC matrix element has been presented mostly 
for demonstration but the results will be used below in 
subsection MI CI for stripping to resonance states. 

The prior CDCC formalism would be preferable if we 
split the matrix element into the internal and external 
parts over variable r pn to separate the internal and pe- 
ripheral parts of the deuteron bound state wave function. 
But this wave function is well known and is not an ob- 
ject of study. That is why below, when considering the 
stripping to resonance states, we use only the post CDCC 
formalism. 



III. DEUTERON STRIPPING INTO 
RESONANCE STATES 

Now we proceed to the main goal of this paper, for- 
mulation of the deuteron stripping into resonance states 
using the surface integrals what will lead us to the gener- 
alized i?-matrix approach for the stripping into resonance 
states. Let us consider the deuteron stripping 



d + p + b + B. 



(93) 



We assume that the resonance formed in the system 
F = A + n can decay into channel B + b, which can 
be different from the entry channel A + n. We start from 
the post form and transform it to the surface integral 
following the method applied for the stripping to bound 
states. Now the application of the i?-matrix approach 
looks natural. Although we consider the deuteron strip- 
ping leading to a specific final channel d + A^p + b + B, 
there can be a few open channels coupled to the chan- 
nel n + A, which is formed after neutron is transferred 
to the target A. As in the previous sections, follow the 
i?-matrix approach, we split the integration region over 
r n A into two regions: internal and external. Internal re- 
gion is determined as the one where all open channels are 
coupled with each other, so that the transition from one 
channel to another can occur only in the internal region. 
The external region is the one where all the channels are 
decoupled. We obtain new forms for the DWBA and 
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then for the post form of the CDCC amplitude. For the 
DWBA both post and prior approach will lead to the 
same final expression. In the standard approach the post 
form of the DWBA amplitude is mainly contributed by 
the external part in the subspace r n A, where the con- 
vergence question of the DWBA matrix element, which 
contains the integration over r p p and r nA , becomes a 
main issue. In the prior form the main contribution to 
the DWBA matrix element mainly comes from the inter- 
nal region in the subspace r n ^, where a strong coupling 
between different open channels becomes an issue. In a 
new approach formulated below the DWBA amplitude 
(in the post and prior forms) is written as the sum of 
three amplitudes: small internal post and external prior 
forms, and the dominant surface integral in the subspace 
over r n A- This surface term is parametrized in terms 
of the reduced width amplitudes, resonance energies and 
boundary condition, that is the quantities used in a stan- 
dard i?-matrix approach. In the post CDCC approach 
the amplitude is given by the sum of the small inter- 
nal post form and the dominant surface term, that is, in 
contract to the DWBA, no external prior form appears 
in the CDCC method. This resolves the issue of the con- 
vergence for stripping into resonant states. 



A. Stripping to resonance states. Post form of 
DWBA. 



The post form of the DWBA amplitude can be ob- 
tained by generalizing the corresponding equation for the 
deuteron stripping to the bound state. As a starting 
point, we use Eq. (J7]) in which, to get the amplitude 
for the deuteron stripping to resonance states, we should 
replace the overlap function 1% by the exact scattering 

wave function with the incident wave in the channel 



b + B: 



M DW(post)^ k 



dA) 



,(+) 



where AV p f = U p a + V pn — U p f and 



OtS k bB —K-bB 



(94) 



(95) 



Since we consider the stripping to the resonance state, 
which decays into two fragments b and B, there are three 
particles, p, b and B, in the final state. Hence, the kine- 
matics of the final state of the reaction depends on two 
Jacobian momenta, for which we adopt the relative mo- 
mentum of two fragments b and B and by the momen- 
tum corresponding to the relative motion of the exit- 
ing proton and the center of mass of the system b + B. 
Thus the deuteron stripping reaction amplitude depends 
on the momentum P — {k p F, k^s}, which is the six- 
dimensional momentum conjugated to the Jacobian co- 
ordinates of the system p + b + B Y = {r p p, r^}. 



Then repeating the steps used in derivation of the ex- 
pression for the post form of the DWBA amplitude for 
deuteron stripping to the bound state we get 



M DW( P ost)^ ^ 

_l_ M DW (P rior ) ( p 1, 
+ M ext l-n k dA 



M^ w (P,k dA ). 



(96) 



Here, internal post amplitude M®Y^ post \p, k^) and 
external prior amplitude ^ w%or ^ {P , k^) are given 



by 

M™ [p ° St) {P, *dA 



^ X p F 



\- 1 J 

Tlf^\AV pF \^l> <n 07) 



and 



DW (prior) 



(P, kdA) 



z < X pF 



(-) r^(ext)(-) ] 



AV pF \<PdXdA > 



r„A>RnA 



(98) 



Here, T 



(*«*)(-) 



•.(ext) (- 



(r„A) 



\r nA ) =< ip A \y 



=< 

(exi)(-) 
bB 



(»"*)(-) 
bB 



> and 



>. 



The last term of Eq. (|96|) . which will be transformed 
to the surface integral, is 



M DW 



dA 



ApF L nA 



-?\<Pdx'&> 



l>Rn 



(99) 



Let us discuss the advantage of this new form of the 
DWBA amplitude for the deuteron stripping to reso- 
nance state(s). Since the internal part is 
given by the volume integral, its calculation requires 

the knowledge of ^S> bB in the internal region. The 
model dependence of this function in the nuclear inte- 
rior (r n A < Rua), where different coupled channels do 
contribute, brings one of the main problems and main 
uncertainty in the calculation of the internal matrix el- 
ement. However, as it has been discussed in subsection 
III Al this matrix element gives a small contribution to 
the total post form amplitude M c ^ p ° st ! due to the 
structure of the transition operator A V p f and constrain 
r nA < RnA- These arguments are also valid when con- 
sidering the stripping into resonance states. A proper 
choice of the optical potential U p f and the channel ra- 
dius R n A may significantly reduce the contribution from 
the internal post form DWBA amplitude. Due to the 
structure of the transition operator AV dA, which has 
been also discussed in subsection III Al the external ma- 
trix element M® x Y^ p " or ^ in the prior form is also small 
and in some cases with reasonable choice of the channel 
radius R n A even can be neglected. Note that in order to 
keep small M^j pos *^ the channel radius R n A cannot be 



DW {prior) 



cannot 



too large and in order to keep small M ( 
be too small. Thus with optimal choice of the chan- 
nel radius the dominant part is the surface part Mg W , 
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which contains only one volume integral over r p p. Eq. 

which presents a new form of the DWBA ampli- 
tude for stripping to resonance states, is quite important 
for analysis of the stripping to resonance In this sense 
the usage of the external prior amplitude M^^ p " or ' 
has clear benefit because it is small and better converges 
the external post form. Also small is the internal ampli- 



tude M 



DW(post) 



The main contribution to M DW( - post ^ 



comes from the surface term Ml? . Using the i?-matrix 
representation of the scattering wave function ^ bB we 
are able to express the total DWBA amplitude in terms 
of the reduced width amplitudes, level matrix, bound- 
ary condition and the channel radius, that is parameters 
used in a standard i?-matrix method to analyze binary 
resonant reactions n + A — > b + B, Since the reaction un- 
der consideration is the deuteron stripping, the presence 



of the deuteron in the initial state and exiting proton 
causes the distortions. That is why the reaction ampli- 
tude, in addition to the i?-matrix parameters describing 
the binary subprocess, contains additional factors - dis- 
torted waves in the initial and the final state. That is 
why we can call the obtained expression for the DWBA 
amplitude a generalized i?-matrix for deuteron stripping 
to resonance states. 

Now we proceed to the derivation of the expressions for 
each amplitude in the right-hand-side of Eq. (|96| and the 
total post form DWBA amplitude. Since the stripping 
into resonance states can lead to rearrangement, the exit 
channel b + B may differ from the entry channel n + A. 
To proceed further we now use the equations for V^ig 
obtained in Appendix [A] Taking into account Eqs. 
and (|A1[) we get 



J 



M, 



DW(post) 



(P, 



27T / kbB 

kbB V Mfcs 



N 



^ < sm s Ittii\Jf Mf > < J n M n J p M p \ J d Md > 



J F M f , l mi M„ 



Y; mi (-k b B) £ [T^BsuAEbB^lA-'U <X pF ) ^rnAn^V pF \^ dXd + 2 > 



r nA <R n 



(100) 



r 



In this equation we assume that the channel spin s and 
its projection m s in the exit channel c = b + B are fixed 
4l|. Jf is the resonance spin (Mp its projection) in the 
subsystem F = n+A — b+B and I is the b+B orbital an- 
gular momentum in the resonance state. The sum over 
Jf and I assumes that a few resonances with different 
spins may contribute to the reaction. The subscript c 
used in Appendix [A] for the channel b + B is replaced 
here by bB. Also S r 7 ^^ F =< tp A \X*' M » > is projec- 
tion of X^ F Mf introduced in Appendix [A] on the bound 
state ipA- The bound-state like wave function X ,1fMf 



describes the system F = n + A = b + B in the internal 
region. A priori, it can be calculated using, for example, 
the shell model approach (36|. In Appendix IA1 X^ F Mp is 
written as a nonorthogonal sum of coupled channels, see 
Eq. (|A4j) . If we neglect the contribution from the chan- 
nel c, then 5^^ F can be approximated by the internal 
part of the overlap function, see Eq. (p~8|) . Taking into 
account this equation (rewritten in L^-coupling scheme) 
we get 



1 . T DW(post),- D , \ 2"7T / k b B 

M mt >{P, k dA ) = — J — 

KbB V f^bB 



< s m s lmi\JF Mf > < s m s > I' mp \Jf Mf > 



E 

Jf Mf s' 1 1' m s r mi M n 

x < J n M n J A M A \s'm s , > < J n M n J p M p \J d M d > e~ iS ™> Y* mi (-k bB ) 

N 

(E bB )} 1/2 [A-X T <x P F ) Y; m Ar nA )I^ s , l , J Ar nA )\AV pF \^ d X dA ) > 



' > [1' ubB slJ F 



r„ A <R„A 



(101) 



Here we added the sum over the channel spin s' (its pro- bital angular momentum V . The sum over M n and s' 
jection m s ' in the entry channel c' = n+A of the resonant appears because the transferred neutron is intermediate 
subreaction n + A — > F — > b + B and over the n + A or- (virtual) . It is important that with a proper choice of 
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the optical potential U p p the matrix element M int 
can be minimized so that its model dependence wouldn't 
have impact on the total matrix element M DW ^ post ^ . 
To obtain the expression for prlor ^ we use for 

(ext)( — ) 

the external part ^Ib , which can be obtained from 
Eq. (|A31j) . assuming that the resonance contribution to 
this wave function is dominant. In the sum over Jp in 
Eq. (IA31[) we keep only those total angular momenta at 



which resonances contributing to the reaction occur. Let 
us consider two possible cases. 

(i) The exit channel c = b+B in the resonant sub-process 
n + A — > b+B is different from channel c! = n + A. In this 
case the external resonant wave function is given by Eq. 
(IA36[) and its projection on the bound state £ c < = ip A is 
determined by Eq. (|A37I) . Then {prior) recmces ^ 



M, 



DW (prior) , „ , n . 2 7T / VbB V"^ 

'(P, k dA ) = - 1 — J > 

KbB V V nA ^ 



i l < Imi s m s | Jp Mp > < I' mi' s' m s i \ Jp Mp > 



J F Mp s' I V m s i mi M n 



x < J n M n J A M A \s'm s , >< J n M n J p M p \J d M d > Yf mi (-k bB ) S J bB s lnAs , v 
x < xff OUkn r A : ^ V mi , (r nA )\ AV dA \if d X f2 > 



r nA 



r n A>RnA 
I 



(102) 



Here, V dA is given by Eq. (|35[) . In the external region 
V nA — and V dA — U pA — U dA . Also has been added 
the sum over the orbital angular momentum / and its 
projection mj (I' and mp) in the exit (entry) channel 
c = b+B (c 1 = n+A) of the resonant subreaction n+A -+ 
b + B, the sum over the channel spin s' and its projection 
m s i in the entry channel c' = n + A of the resonance 
subprocess n + A -+ b + B and the sum over M n because 



the neutron is transferred particle. The projections of the 
spins of the incident deuteron M dl the exiting proton M pi 
the channel spin s and its projection m a of the exiting 
particles b and B are fixed. We also use the symmetry 



of the S matrix: 



element S b g 



c' s' l';csl 



s 



c s l\c' s' V 



The matrix 



A i v is given byEq. (|A"4"5"|) . Substituting 
it into Eq. (fT02|) gives 



2. 

kbB V V nA 



E 



i < I mi s m s I Jf Mf > <l ! rnv s' m s > I Jf Mf > 



Jf Mf s' 1 1' m s i mi M n 

x < J n M n J A M A \s'm s , >< J n M n J p M p \J d M d > Y,* (-k bB ) 



xe -«fc. e -«; 



"„ (jp \il/2r A -li r r i t? yil/2 Ol'jknA, R nA ) 



(_) OjAknA, TnA) R nA m a T7 i (+) 

x <X pF > ^- '-j^jt j^Y; (r nA )\AV dA \^ d x dA > > 

TnA ^I'yKnA, -KnA) 



r nA >R-nA 



(103) 



Now we take into account that 



Oj(k 5 , R 5 ) = X /Ff(k d , R~ c ) + Gj{k~ Cl Rg) 



where in the absence of the Coulomb interac- 
tion Fi{p) = (^P/2) 1/2 J /+1/2 ( P ) and G t {p) = 

(-!)' ( 7r P/2) 1/2 ^-(/+i/2)(p), J±(i+i/2)(p) are Bessel 
functions. 



F?(k d , R z ) + Gf(k 5 , Rs) 



(104) 



which for the channel c = c' = n + A and / = I' takes the 
form 



Oi>(k nA , R nA ) = ^F^(k nA , R nA )+Gf,(k nA , R nA ) 



F,,{k A ,R A ) 

% arctan ^ > , ~ r 

X e I' ( nA ' R nA> 



Then using Eqs. (|A4ip and (|105p we get 



= JF{;(k nA , R nA ) + Gf,{k nA , R nA ) , 



(105) 



18 



M. 



DW (prior) 
ext 



(P, k, 



dA 



2tt 



[i-bB kbB Rni 



E 



i <lmi sm s \JF Mp > 



Jf Mjp s' I V m s i mi ■m l i M n 

x < I'mp s' m s ,\J F M F > < J n M n J A M A \s' m s > >< J n M n J p M p \J d M d > Y* (-k bB ) 



iV 



x e 



-is 



iA s' V J 



E { r -bBsljAEbB)} 1/2 [A-\ TlTn , 



U,T — 1 



(_) Of,(k n A, r n A) RnA 
X < Xp F — 



Of, (knA, RnA) 



r„A>Rn 



(106) 



(ii) If c = c', that is b = n and B = A. Here two cases 
are possible: non-diagonal transition for which s =/= s' 
or/and I ^= I' and diagonal transition with / = I' and 
s = s'. The amplitude for the nondiagonal transition can 



be obtained from (|102|) . Here we present the expression 
for the diagonal transition (elastic scattering) amplitude, 
which can be obtained taking into account Eq. (|A33I) : 



M. 



DW (prior) 
ext 



(P, k dA ) 



2n 



k n A Rn/ 



E 



i < I mi s m s | Jf Mp > < I m r s 77 V I Jf Mf > 



Jf A4f I rn s t mi m>/ M n 

x < J n M n J p M p \J d M d >< J n Mn J A M A \s m s > > Yf (-k nA ) 



-i _ qJf 

1 °( n A)sl;(nA)sl 



Ol(k n A, RnA) 



(-) Of(k n A, r n A) 
x < XpF — 

P TnA 



Rn/ 



■Y*{VnA)\AV dA \ VdXd + 2 > 



Of (k n A ) Rn 



r nA >Rn 



(107) 



Substituting the expression for the elastic scattering S- obtain 
matrix element S^ A) sl .^ nA) sl given by Eq. (|A43j) we 



M, 



DW (prior) 



ext 



(p, k dA ) = i 



2tt 



k n A RnA 



E 



i < I mi s m s I Jf Mf > < I tup s m s > I Jf Mf > 



Jf Mf I rn s / mi m^ M n 

x < J n M n JpM p \J d M d >< J n M n J A M A \sm s , > Y* m (-k n A) 



N 



X 



l_e- 2<4 »i«(l+t ^unAsljAEnA^^^XrlTrnAsuAEnA)] 1 ' 2 ) Ol(k nA , R u a) 



„ (-) Of (k nA , r nA ) R„a v * /- mat? i (+) 

X < XpF TFJl 7r^ Y l^nA)\AVdA\^dX y d A > 

•nA ^ i \ K nA, n n A) 



r nA >R„ 



(108) 



One-level, one channel case is the simplest one for which 

M DW{prior)^ ^ ^ dQwn ^ 
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M OW^ior) {p ^ dA)=i 



2tt 



k n A RnA 



E 



i < I mi s m s | Jp Mp > < I m\t s m s > \ Jp Mp > 



Jf Mp I m s i mi m ; / M n 



x < J n M n J p M p \J d M d >< J n M n J A M A \sm s , > Y t * {-k nA ) 



1 - e - 2lS *Asu F e 



2i<5„ 



(-) OfjknA, r nA ) 
x ^ Aj,_F 



Rn 



r nA Oi(k n A, RnA) 



Y*(v nA )\/\V dA W d X d + 2 > 



r nA >R„A 



Oi(k n A, RnA) 

(109) 



where 



JnA s I J F 



arctan 



nA s I J F 



(E nA ) 



2{E n A(o)si j F 



E 



nA(0) si J F 



> E n / 



(110) 



is the resonant phase shift, E n A(o) sij F is the real part of 
the complex resonance energy of the resonance with the 
quantum numbers s I Jf in the channel n + A. Now we 



derive the equation for Mg by transforming it into the 
surface integrals over variable r n A- We can repeat the 
discussion in Section II A. The integration in Eq. (|9"9")l 
over r n A is taken over the external volume restricted by 
two spherical surfaces: the inner surface with the radius 
R nA and the external surface with the radius R nA — > oo . 
As it has been shown in Appendix [B] after regularization 
the integral over the infinitely large sphere vanishes (see 
Eq. (HS|) and 



Mg (P, k, 



dA 



-M^(P,-k dA )=Rl A 



1 



2 [InA 

(+) 



dr„ 



A+) /_ UT (ext)(-), m d^VOXk^X 1 "^) . 
X-k pF \ r pF) [l nA \*nA)\ TT- . 



r nA =R n A 



(111) 



Here, —MP W (P, k d A) is the surface integral encircling 
the inner surface of the external volume at r n A = R n A- 
A negative sign appears because the normal vector to the 
surface is directed to the center of the volume, i.e. op- 
posite to the normal vector to the external surface (at 
infinitely large radius). For simplicity, we dropped the 
quantum numbers in Eq. (|111|) but they will be recov- 



ered below. Note that Eq. (|111[) can be obtained from 

Eq. (|2"7|) by substituting (T n A) f° r the overlap 

function 1% Jf Mp Ja Ma mnA (r nA ). 

For the exit channel c = b + B in the resonant sub- 
process n + A — > b + B different from channel c' = n + A 
using Eq. (|A37[) we get 



Ms° w (k pF , k dA ) = ~M™ (P, k dA ) 



2tt I VbB 



R 



■bB V v nA 



nA 



E 



< I mi s m s \JF Mp > 



Jf Mf I I' mi s' M n 



x <l'm v s m s , \JM>< J n M n J A M A \ s m a > >< J n M n J p M p \ J d M d > Y t * mi (-k bB ) s , v . bBsl J d r pF X { ^ pF ( r pF ) 



x / dO r „ A yj/ mi ,(f nA ) Vd(r P n)Xk, A (*dA) 



q OA^a) ,O v {k nA ,r nA )dip d {r^)x [ ^ A ){vdA) 



dr n/ 



x ( X F M 



r n A 



r n A 



r„A=R r , 



27T I VbB 1 



kbB V V nA 2 //„, 



i l < Imi s m s \JpM > < I' mi' s' m s > \ Jf Mf > < Jn M n Ja M a \s' m s > > 



J F Mil' mi ra i' s ' M n 



x <J n M n J p M p \J d M d > Y l * mi {-k b B)S^ sl . nAs , ll Oi,(k nA , R nA ) I dr pF x^ p (r pF ) / dCl InA Y v mi ,(r nA ) 



\n (r wWfr UR \\ R 9 ^^") ^£)(^) 1 

fd{r pn ) Xk d A r dA) \B nA — 1) - ri n A q 

O TnA J 



r n A=RnA 



(112) 
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Here. 



B n A = R 



nA 



dO t i (k n A,r„A) 




dr nA 


r n A=R,iA 



(113) 



Ol'(k n A, RnA) 

is the boundary condition. Sum over M n is a for- 



J 



mal because M d and M p are fixed. The coefficient 
< J n M n J p M p \J d M d > appears from the vertex d — > 
p + n and the product < V mi* s' m s i\JF Mf >< 
J n M n J a Ma\s' m s > > from the vertex n + A — > F. The 
matrix element S^£ sl As , v is given by Eq. (|A45p . Sub- 
stituting it into Eq. (|112l) gives 



M™{k pF , ^a) = -Mjf ' (P, k dA ) = f - J^- — 

k hB V VnA VnA 



V] i l <lmi si71 s \JfMf > 

Jf Mf 1 1' s' mi m^ m s i M n 



x < I'mv s'm s ,\J F M F ><J n M n J A M A \s'm s , ><J n M n J p M p \J d M d > YJ* m| (-£«,) e - M »«' e"* 5 '"' 

i/ fc.B s I Jp 

(E bB )}^ [A- 1 ] (£„A)] 1/2 Oi'(fc„A, P„a) Ur pF x% F (r P F) / d^Y^ 



<Pd{r pn ) Xk dA (*dA) (BnA ~ 1) ~ PnA 



<9^(r p „)xLti)( r ^) 



r 



Taking into account Eq. (|A41[) and Eq. (I104j) we arrive at the final form for Mg W (k p p, k dA ) 
I 



(114) 



M° w {k pF , k dA ) = -Mg^JP, k dA ) = 



— — i l < I mi s m s I Jp Mp > 

VbB^nAkbB . ,„ , 

Jp Mp 11 s mi m ( / m s i M n 



x < V m v s' m s , I J F M f > < J n M n J A M A \ s' m s , >< J n M n J p M p \ J d M d > F z * mi (-k 6B ) e~ 1 ' 

jv r r 

v bB s I Jp 

(E bB )] 1/2 [A \ T jrnAs'i'j dr pF x ( -2 pF (r P F) / dCl rnA Yi> 



■ , x (+), wo ,N p gyrf(^n)Xktj)( r ^) 
</Wjm) Xk d i(l"dA) (-B nj4 - 1) - i?nA g r ^ 



r„A=R^ 



(115) 



Now let us consider the diagonal transition csl csl, be replaced by ^csTro -cs i m + ^cslm ^-c L m given by 
where c = c' = n + A To get once again we start Eqs . (TAlO)) and (fAlt)) . Then the equation for surface 

from Eq. (jlllj) . Now in this equation Y^*'*- - ' should matrix element for the diagonal transition takes the form 



Mg w (k pF , k dA )=i 



fX n A k 



nA 



E 



i < I mi s m s | Jf Mf > < I mi" s m s " \ Jf Mf > 



Jf Mp I mi min m s ii M n 

x < J n M n J A M A \sm s „ > < J n M n J p M p \J d M d > Yi* mi {-k nA ) 



N 



1 - e 



-i2 52 



v nA s I Jp {E nA )] 112 [A- 1 ] 

vt J- t nA s I Jp (E nA )} 1/2 )\ Oi(k nA , R nA ) 

U,T — 1 



d r pF X ( -k„ (r p F) / d O r „ A Y x „ H „ (r nA ) (fi d (r pn ) x^CJa ( Y dA) (B nA ~ 1) - RnA 



(+) 



d^ d {r p n)xC dA )^dA) 
dr n A 



r„A = RnA 

(116) 
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Summing up all three amplitudes M int (pos ' (P, k dA ), 

M°J {pn ° r \P, k dA ) and M° w (k pF ,k dA ) 
-M^ A (k pF ,k dA ) we get the total post DWBA 

for the {d,p) stripping. 



(i) Resonant reaction n + A -+ b + B, that is 
c = b + B ^ c' = n + A. Then the total post 
form of the DWBA deuteron stripping amplitude is 



HbB h 



bB 



E 



i 1 < s m s lmi\Jp Mp > < s f m s > V my \Jp Mf > 



< J n M n J A M A \s'm s ,><J n M n J p M p \J d M d > e~ iS ™i *T mi (-k«j) £ [T ubBsl j F (£ 6B )] 1/2 [A" 1 ], 



N 



X | <X p Jl% s >l>jA r nA)\AV pF \ VdXd + 2 > 

(-) OljknA, r nA ) R nA 



, 2 U n A 

+ \ -5 JrnAs'l'J 

r„A<RnA V tinA 



X < 



X P F 



yYir mi ,(f„ A )|Ay pF |^x^ > 



r nA 



Oi,(k nA , R nA 



r„A>Rn 



2 flnA 



mil (*V. 

< \ (+) 1 wo in o 5 Vd(^n)Xkti)( r ^) 

¥>d (r P n ) Xk d l ( r <M ) (B nA - 1 ) - i?„.4 — 



<9r ny 



r nA =R„ 



(117) 



Assuming in this equation b = n and B = A, that is 
c = c' but Z 7^ /' and/or s ^ s' we get the expression 
for the DWBA deuteron stripping for the non-diagonal 
transition in the resonant subprocess (n + A)i s — > F — > 
{n + A) v s ,. 

Equation (| 1 1 7[) is very instructive for understanding 
the difference between the stripping to resonance states 
and on-shell binary resonant reactions. As we can see, the 
transfer reaction amplitude contains the resonance fac- 
tors determining the resonant subprocess n + A — > b + B, 
the partial width amplitude [T„i,Bsij F (E bB )} 1/2 of the 
level v for the decay to the exit channel b + B, the matrix 
elements of the inverse i?-matrix level matrix [A _1 ]„ r 
and the reduced width amplitude j T nA s * v j of the level 
r for the entry channel n + A rather than the correspond- 
ing partial width amplitude which would present if we 
consider the corresponding on-shell binary resonant re- 
action n + A -+ b + B. The difference is crucial because 
the partial width amplitude [T^bBsi j F (-EbB)] 1 ^ 2 contains 
the penetrability factor, see Eq. (|A41I) . which is missing 
in the reduced width amplitude and, hence, in Eq. (Ill 7[) . 
The lower is the energy of the resonance, the stronger is 
its suppression due to the barrier penetrability in the en- 



trance channel in the on-shell binary resonant reaction 
n + A — > b + B. Besides, if a few resonances do con- 
tribute with the different I', then the higher is I', the 
stronger is its suppression. However, it is not the case 
if one tries to populate low-energy resonances with dif- 
ferent I' using transfer reaction. Missing penetrability 
factor in the entry channel of the subresonance reaction 
n + A -+ b + B in the transfer amplitude makes it pos- 
sible to populate low-lying resonances. Moreover for the 
same reason, the resonances with higher I' are not sup- 
pressed in the stripping. Hence, when a few resonances 
are populated in the transfer reaction, the measured ex- 
perimental spectrum of the fragments b and B can be 
quite different from the one measured using the on-shell 
binary resonant reaction. The missing penetrability fac- 
tor in the entry channel n+A of the resonant subreaction 
n + A — > b + B in the transfer reaction explains the power 
of the Trojan Horse method as indirect technique in nu- 
clear astrophysics (see [HI, [13] and references therein). 

(ii) Diagonal transition in the resonant subprocess {n+ 
A)i s -» F -> (n + A)i s , that is, c = d , I = I' , s = s'. 
The total post form of the deuteron stripping DWBA 
amplitude is 
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M DW( P ost) (p ; k dA ) = 2 7T il < s TO « 1 m i I J F M F ><s m s „ I m V > \ J F Mp > 

Jf Mp I m s n nil m ;" M n 

x < J n M n J A M A \sm s „ ><J n M n J p M p \J d M d > e -' lS ^< y* (-£„a) 



\ J ^7— Y,^^ sl j F (EnA)} 1/2 [A- 1 U <X { ; F ] I F A sl j F {r nA )\&V pF \v d X d + 2 > 

V VnA K n A y 
N 

i _ e -i 2 , ^ + , ^ [r„ nA a , j, (£„a)] 1/2 [A- 1 ],. r r nA , , Jf (£„a)] 1/2 ' 



r„A<-R„ 



■ (M^nAj tt-nA) < X pF 



Rn 



k n A RnA r nA 0*(k nA , R nA ) 

dtpd{rpn)x ( ^ A ){*dA) 



Y*{i nA )\&V pF \ Vd x d + 2 > 



r„A>RnA 



^UnA 



¥>d(r p n) Xktl( r ^) ( B nA - 1) ~ RnA 



dr n/ 



r nA =R 




(118) 



B. Stripping to resonance states. Prior form of 
DWBA. 



Here we show that starting from the prior form we are 
able to obtain the generalized DWBA i?-matrix ampli- 
tude for the deuteron stripping to resonance states, Eq. 
(|96p . much easier than from the post form. The prior 
of the DWBA amplitude for deuteron stripping to reso- 
nance states is 

M DW(prior)^ ^ 

=<X$*$\\*V dA \<Pd<PAX < & >, (119) 

where A V dA is defined by Eq. (|35l) . As usually, we split 
the amplitude into internal and external parts 



+M°Y iprior) (p^ dA ) 



(120) 



with 



M t D J ip "° r \P, k dA ) 



< XpF 



(-) $ M)H 



\AV dA \ip dXd + 2 > 



IB 



r nA <R n 



(121) 



and 



M, 



DW (prior) 



=<x& ) *sr )( - , iA7 <U | WX £i ) > 



r nA >Rn 



(122) 



The splitting of the amplitude into the internal and ex- 
ternal parts in the subspace over the coordinate r nA is 
necessary to rewrite the prior DWBA amplitude in the 
generalized i?-matrix approach for stripping to resonance 



states. As we have discussed in subsections III Al and 



IIII Al the external matrix element M 



DW (prior ) 



the 



prior form is small and in some cases with reasonable 
choice of the channel radius R nA even can be neglected. 
It is important for analysis of the stripping to resonance 
states because the external part in the post form doesn't 
converge. In this sense the usage of the prior form in 
the external part has clear benefit. The main contribu- 
tion to the prior form amplitude M DW ^ rior ^ comes from 
the internal part M^ /< " p " or \ Since the internal part 
is given by the volume integral, its calculation requires 

(int)( — ) 

the knowledge of ^ bB in the internal region. The 
model dependence of this function in the nuclear interior 
{fnA < Rua), where different coupled channels do con- 
tribute, brings one of the main problems and main uncer- 
tainty in the calculation of the internal matrix element. 
Using the surface integral we can rewrite the volume inte- 
gral of the internal matrix element in terms of the volume 
integral in the post form and dominant surface integral 
taken over the sphere at r nA = R nA . With reasonable 
choice of the channel radius R nA the contribution from 
the internal volume integral in the post form can be min- 
imized to make it significantly smaller than the surface 
matrix element. The latter can be expressed in terms of 
the i?-matrix parameters - the observable reduced width 
amplitude (ANC), boundary condition and channel ra- 
dius. Repeating the steps outlined in subsection III Bl we 
get 



M, 



DW [prior) 



(P, k dA ) 



„ r DW(post) 
= M int 



(P, k dA ) + M° W (P, kdA). 



(123) 
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Here, M^ r ^ post ^ has been previously considered and is 
given by Eqs. ([TOT) and ([TUT]) while M^ w takes the form 

rDW 



M^ w (P,k dA ) 



(-) r (ext)(-), t; 7^1 



r nA <R nA 

(124) 



where ?(<**)(-) =< *^ t)(_) |^ >. The fact that the 
volume integral in this equation is the internal one makes 
transformation of this volume matrix element to the sur- 
face one much easier than for the post form. The tran- 
sition operator T = T p p + T nA . Since r nA < R n A at 
r p p — > oo the integrand in Eq. (|9Q[) vanishes exponen- 
tially due to the presence of ip d . Hence, the operator T p p 
is Hermitian, that is, applying the integration by parts 
over r P F twice we get 



< XpF 1 I J ~ J I fd X dA > 

= < x£ > < > | 2* - 2* I WxS > 



0. 



r nA <R nj 

(125) 



Thus Mg w {P, k dA ) reduces to 



2* 



nA \VdX [ dA > 



r nA <Rn 



(126) 



Using the Green's theorem we can transform this vol- 
ume integral into the surface one. Note that the volume 
integral over r nA is constrained by the sphere with the 
radius r nA = R nA . Hence, only one surface integral ap- 
pears with r nA = R nA . Here we see an important ad- 
vantage of using the prior form versus the post one. In 
the post form transformation of the external volume in- 
tegral to the surface one led to two surface integrals at 
r nA = R nA and r nA = R nA — > oo. It required an elabo- 
rate proof, which included regularization, to demonstrate 



that the surface integral at r nA = R 



oo vanishes. 



After transformation to the surface integral we get 
Mg W (P, k dA ) = -M™ (P, k dA ), 



(127) 



Eqs (jllip . (|112p and (|115l) determine this surface inte- 
gral. 



Stripping to resonance states. Post CDCC 
formalism. 



The CDCC approach for stripping to resonance states, 
which takes into account the deuteron breakup in the 
initial channel, definitely has advantage compared to a 
standard DWBA. The application of the surface formu- 
lation of the reaction theory for the DWBA has been 



done mainly for demonstration, but our main goal is the 
CDCC. 

Here we present the derivation of the post form CDCC 
amplitude using the surface integral. This amplitude is 



M CDCC ^ ost \k pF , K-dA 

=<x^ ) *fe t)( ~ ) IAF P ' 



ypF *bB 



pF 



^CDCC(+) > 



(128) 



This equation is an extension of the post CDCC am- 
plitude for stripping to bound states, see Eq. (|54[) . 

p _ 

A V pF " is defined by Eq. ([7TJ|) . Now, as usually, we split 

j^jGDCC(post) j n £ fa e m ternal and external parts in the 
subspace r nA : 



M CDCC(jpost)rp^ k| 



dA 



rCDCC{poat) , 
int 

rCDCC(post) , 
ext 



M^?^^>(P, k dA ) 



■M^^" uab) {P, k dA ) 



(129) 



The internal amplitude ^^^ cc vP° st ) j s given by 



CDCC(post) 



(P, k dA ) 



z < Xp F * bS 



r nA <R„ 



(130) 



Correspondingly, the external amplitude is 



M. 



CDCC(post) 



(P, k dA ) 



^ XpF ^bB I 



AV ^CDCC {+)> 



r nA >RnA 

(131) 



Now we repeat the steps outlined in subsection lll CI Tak- 
ing into account Eqs (|55j). (151?)) and ([BTJ|) we arrive at 



M. 



CDCCipost)^ = M CDCC{post) {p ^ ^ 



<X ( p F ) <T ){ - ) \T-^f DCC{+) > 



r nA >R nA 

(132) 



where T = T p p + T nA . It is shown in Appendix [Cl that 



j^CDCc{post) ^ e reduced to 

T CDCC(post) i 



""""'(P, k dA ) 

v (_) \n( ext )(.-)\^r . _7* , hi/ 
_< XpF ^bB I 1 nA J- nA\ ^ 



CDCC(post) 



M 

CDCC(+) 



(P k d A) 



> 



r nA <R nA 

(133) 



This integral can be directly transformed into the surface 
integral with r nA = R nA encircling the internal volume, 
while the integral over r p p is taken over all the coordi- 
nate space. Thus we have shown that the post CDCC 
amplitude for stripping to resonance states is given by 
the difference of two terms, internal post CDCC ampli- 
tude and the surface integral: 



M CDCC(post) f p^ ^ 



dA 



CDCC(post) 
int 



-M 



CDCC(post) 



(P, k dA ) 

(k pF , k dA ). (134) 
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The internal amplitude ]\/[CDCC(post) can ^ e m j mm j zec j 
by a proper choice of U pF and the channel radius R nA , 
while the surface integral is dominant. If the channel 
radius is larger than the n — A nuclear interaction ra- 
dius the second term is parametrized in terms of the re- 
duced width amplitude and the boundary condition at 
T n A = RnA- Thus we succeeded to parametrize the post 
form of the CDCC amplitude in terms of the i?-matrix 
parameters. It is one of the main results of this paper. 
Eq. (|134[) is the most important result of this paper. 
Due to the absence of the external term, which is present 



in the DWBA and which causes the convergence issue, 
the convergence problem in the post CDCC approach is 
resolved: the integration in the surface matrix element 
is performed over the full coordinate space only over one 
coordinate r pF rather than over two coordinates, r pF and 
r n A- 

Expression for ^[Gi)CC(%iost) ^ different cases can 
be obtained from Eq. (jlOip by replacing the 
initial channel wave function (pd(r pn ) Xit A )( T dA) by 
V i ; (r pF , r nA ): 



M. 



CDCC(post) 
int 



(P, k dA ) 



2vr k, 



bB 



E 



i < sm s I mi\J F M F > 



kbB y fJ-bB 

x < s'm sl l'm v \J F M F > < J n M n J A M A \s'm s , > e -*'& ' Y t * mi (-k 6B ) ^ [T vbB si jA^b)] 1 ' 2 [A" 1 ] 



Jf Mp s' I V m s f mi m l / M n 
-iSH 



N 



< Xff I F As'Vj F {r nA ) \AV pF \*f DCC{+ \r pF , r nA ) > 



r nA <RnA 



(135) 



Note that the CDCC wave function itself also depends on 
quantum numbers of p — n and d — A subsystems, which 
we don't specify here. It will be done in the following 
up paper where concrete calculations will be presented. 
Natural Jacobian variables for ^ DCC ^ +S> are VdA anc i 
Tpn, but we use here another set of Jacobian variables, 



r pF and r nA . 

To write down explicitly Mg^ >CCt ' post \'k pF , k^) in 

terms of the surface integral we can use Eq. re- 
placing the initial channel wave function by the CDCC 
one: 



M 



' (J'TlA 



(k pF , k dA ) = -M^ CC ^ ost \k pF , k dA ) 



CDCC(post) 
S 

/! '" ! |dr pFX W F (r pF )|d^ 



dr pF X- + kL( r P F ) / dn r nA f nA 



,(ext)(-) 



CDCC(+) 



r„A<R n A 



*(exi)(— )* 



^cdcc(+) ( \ 9T nA 



(eat) (-)« OW i '{r pF ,Y nA ) 

TlA 



dr n/ 



(136) 



We can extend corresponding equations from subsection in the resonant subreaction d s 1 1' — > cs I, where c = b+B 
IIII Al bv replacing the initial channel wave function by the and d = n + A, we get from Eq (| 1 1 5[) 
CDCC one. In particular, for the nodiagonal transition 
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M 



CDCC(post) 



(k pF , k dA ) = _ M ™ CC ^° St \P, k dA ) = 7T 



2 R nA 



V-bB UnA kbB 



E 



Jf Mf I V b 1 mi m,ii m B i M n 

x < I mi sm s \J F M F >< I' m v s' m s >\J F M F >< J n M n J A M A \s' m s > > 

n r r 

X^f-kiBje-'S' [^^BsUF^EbB^^^-XrlrnAs-VJ J d V pF X ^ (r pF ) j d fi r „ A Y,, m ,, (f nA ) 



f ,T — 1 



*f DCC(+) (V, r«A) (BnA - 1) - RnA 



d*f DCC ^\r pF ,r nA ) 
dr nA 



r nA =R„ 



(137) 



r 



Correspondingly, the surface integral for the diagonal transition c si ^ c si can be obtained from Eq. (|116l) : 
I 



M 



CDCC(post) 
S 



(k pF , k dA ) 



l^nA k nA 



E 



i < I mi s m s | Jf Mp > < I ra\n s m s n \ Jp Mp > 



x < J n M n J A M A \sm s „ > Y* (-k nA ) 



N 



1 - e~"" nAl 
dr pF X-L F ( r pF ) / dn rnA Y lmi „(r nA ) 



1 LH 1 + '(r pF , r nA ) (B nA - 1) - R nA 



Ol(k nA , RnA) 

d^ DCC ^(r pF ,r nA ) 



dr nA 



r nA =R nA 

(138) 



I 

Summing up two amplitudes M^ l f CC ^ post \p, k dA ) (i) Resonant reaction n + A -» 6 + B, that is c = b + B ^ 
and Mj^kj^, k dA ) = -M]?™ (k pF , k dA ) we get the d = n + A. The total post form of the CDCC deuteron 
total post CDCC amplitude forthe (d,p) stripping. stripping amplitude can be obtained from Eq. {HZ}: 

I 



M CDCC( P ost)^ k ^ )(pj ^ = 2n 



1 



[i-bB kbB 



i l < s m s lmi\J F M F > < s m s / 1' mp \ J F M F > 

Jf Mf s' I V ra s t mi m^ M n 
N 



<J n M n J A M A \s'm s , > e " lS ^i Yl * mi (~-k bB ) £ [T v bB . I j F {E bB )] 1/2 [A" 1 ], 



x { < X ( pF } lL>v J P {r nA ) I A V pF \yf DCC(+ \r pF , v nA ) > 



r nA <R nA 



1 RnA 
%HnA 



IrnAs'l'J J ' dr pF X^f^pf) J&^ nA Yf m,, (^tia) 

*f DCC ( + \r pF) r nA ) (B nA - 1) - R nA **< ^ r ^> 

O r nA 



r nA =Rn 



(139) 



Assuming in this equation b = n and B — A, that is c = d but I ^ V and/or s ^ s' we get the expression for 
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the post CDCC deuteron stripping for the non-diagonal (ii) Diagonal transition, c = c', I — I' , s = s'. The 

transition in the resonant subprocess (n + A)i s — > F — > total post form of the CDCC amplitude is 
{n + A) Va ,. 



M CDCC ^ ost \p : k dA ) = 2tt i l <sm s lmi\J F M F ><sm s „lmi„\J F M F > 

Jf Mf I m 8 n mi Ttii/t M n 

x < J n M n J A M A \sm s » > e~ iS ^ Y t * (-k nA ) 



\\l 1 , l r vnAsij F (E nA )} 1/2 [A^U < XpF } I F A sl j F (r nA ) | A V pF \9 

[ V Hn A K nA ' 

N 

v n A s I Jf ut T T n A s I Jf 



T7 ui,CDCC(+) 



(r pF , r nA ) > 



r nA <R„ 



l_ e -«0« 



l —— [ dr pF x^ pF (rpF) / dO rnA y imi „(f„A) 

^ ^n A K nA J J 

xSi cdcc(+) ( w„ i \ n> 0*f DCC(+) (r P F, r„ A ) 

* 4 l r P-F) r nAj ~ J- J ~ -KnA 



dr nA 



(140) 



r 



ables Td A ~ r P F — ^ oo and r p „ ~ r p j? — ¥ oo. But then 
^f DCC ^(r pF ,v nA ) 



Eqs (|139p and (|140p are the final and main result of 
this paper. Both matrix elements consist of only terms, 
the internal post CDCC and the surface term. The in- 
ternal term contains the integration over r nA in the in- 
ternal volume r nA < R nA . Hence, at r pF — ¥ oo vari- 

oo. 

r pF E35l the integral over v pF 
does converge. The same conclusion is true for the sur- 
face integral in which r nA = R nA . Hence, in this matrix 
element also tyf DCC ^ + \r pFl r nA ) ~ r~p and integral 
over r pF converges. Both amplitudes are parametrized 
in terms of the parameters used in the conventional R- 
matrix approach and providea tool to analyze the strip- 
ping into resonance states using generalized i?-matrix 
approach. Finally, both amplitudes, (|139l) and (I140p . 
don't have penetration factor in the entry channel n + A 
of the resonance formation in the resonant sub-reactions 
n+A — > b+B and n+A — > n+A. That is why stripping to 
resonantstates provides a powerful tool to measure reso- 
nances in the subsystem n+A very close to the threshold, 
which can be suppressed in the on-shell binary resonant 
reaction but not in the stripping to resonance states. 



IV. SUMMARY 

The theory of the deuteron stripping populating bound 
and resonance states based on the surface integral for- 
malism is presented. To demonstrate the theory I first 
develop it for the DWBA. Since the DWBA is outdated 
and, definitely, deficient compared to the CDCC, the the- 
ory is extended to the CDCC formalism. The theory is 
applied for stripping to bound and resonance states. The 



eventual goal of this paper is to deliver the theory of the 
deuteron stripping to resonance states within the CDCC 
formalism using the surface integral formulation of the 
reaction theory [38j . Transformation of the volume ma- 
trix element to the surface one (in the subspace over r nA ) 
and i?-matrix representation of the scattering wave func- 
tion of the fragments formed by the resonance decay al- 
lows one to parametrize the reaction amplitude in terms 
of the 7?-matrix parameters used in the analysis of the 
binary resonant reactions. Since the reaction under con- 
sideration is the deuteron stripping, the presence of the 
deuteron in the initial state and exiting proton causes 
the distortions. That is why the reaction amplitude, 
in addition to the i?-matrix parameters describing the 
binary subprocess, contains additional factors - CDCC 
wave function describing the d — A scattering in the ini- 
tial channel (coupled to the deuteron breakup channel) 
and the proton distorted wave in the final state. Hence, 
the approach can be called a generalized R matrix for the 
stripping to resonance states. The advantage of the ap- 
proach is that the reaction amplitude for stripping to res- 
onance states in the post CDCC formalism doesn't have 
convergence problem and is parametrized in terms of the 
same observables as binary resonant reactions. Hence, 
the formalism provides experimentalists a consistent tool 
to analyze binary resonant reactions and stripping re- 
actions populating resonant states extracting the same 
observable parameters, namely, reduced widths (ANCs). 
The power of the method has been demonstrated in the 
analysis of the Trojan Horse reaction 19 F(d, n a) 16 [32l |. 
The numerical application of the method will be demon- 
strated in the following up papers. 
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Appendix A: b + B scattering wave function "J" 



In this Appendix we consider the representation of the 
scattering fl^g wave function used in the i?-matrix ap- 



proach for binary resonance processes 39, 4' 



We consider only two coupled channels c = b + B 
and d = n + A. Also X Jp Mp is an eigenf unction of 
the Hamiltonian describing the compound system F = 
n + A = b + B in the internal region excited to the dis- 
crete level t with the total angular momentum Jp and 
its projection Mp t 42]. A separable form for 
reflects the fact that we consider the b + B interaction 
proceeding through resonance states. The entry channel 
of this scattering is the channel c = b + B. The inverse 
level matrix contains contribution from all N resonance 
levels. In a simple one level case it reduces to the well- 
known Breit-Wigner resonance propagator. All the open 
channels coupled to c contribute to X Jf Mf and deter- 
mine possible exit channel contributions into resonance 
scattering. Hence, in the internal region, where different 
open channels are coupled, X Jf Mf can be written as a 
nonorthogonal sum of these channels [39] : 



1. Internal scattering wave function 

A general equation for the internal wave function con- 
tains the sum over total angular momentum Jp and its 
projection Mp. Since we are interested in a wave func- 
tion describing a resonance in the system F = b+B, 
we consider only the internal wave function at given Jp, 
at which resonance occurs. In the internal region in the 
state with the total momentum Jp, channel spin s (its 
projection m s ) in the initial channel c = b + B the wave 
function ^bB can be written as [3!| 



X 



Jf M f 



E 

c s I m s - 



l 



— W T r, 

r, 



TCJ 



cslms csljp j 



}, 

(A4) 



Kc V ^ c Ml mi 
N 

xr; mi (k c ) Y, [T V csuAEo)] l ' 2 [A- l } VT x J T 

U,T=\ 



Here, E c = E^g and k c = k^ are the relative energy 
and momentum of particles b and B, fi c — fibB, ^uc(E c ) 
is the formal (i?-matrix) partial resonance width of the 
level v in the channel c = b + B, A is the i?-matrix level 
matrix, N is the number of the levels included, a c ; is the 
Coulomb scattering phase shift in the channel c and the 
partial wave I, <5'f is the hard-sphere scattering phase 
shift in the channel c given by 



where £g is the product of the internal bound state wave 
functions of the fragments in the channel c, c = c, c', 
u ds j Jf Atc) is a set of the radial wave functions of the 
relative motion of the fragments in the channel c with 
the channel spin s, orbital angular momentum / and to- 
tal angular momentum Jp in some adopted potential, 
4>~~y F j where rris is the projection of s, is the channel 

wave function (in LS'-coupling). Also A is the antisym- 
metrization operator between the nucleons of the frag- 
i < sm s lmi\Jp Mp >ments in the channel c. We consider only two coupled 
channels, c = b + B and d = n + A. Thus the initial 
channel c can propagate into two final channels c and d 
via the intermediate resonances. Although Eq. (|A4I) con- 
tains the sum over all channel spins s and projections in 
(Al) each open channel, in what follows consider the contribu- 
tion to X Jf Mf only from the channel with fixed channel 
spin and its projection. 

First, let us consider the contribution of the channel 



c s" m s 



a into X Jf Mf . In this channel £ c = ip b ip B and 



x F/» mj „ (r c ) </> c 



(A5) 



srhs 
°cl 



-ui c i + arctan 



Fi{k c , R c ) 

G ; (fc c , R c y 



(A2) 



where F;(fc c , r c ) and Gi(k c , r c ) are regular and singular 
Coulomb solutions of the radial Schrodinger equation, 



t >cs"m s „= Y, < Jb Mb Jb Mb\s" m a H > 
M b M B 



X 1p J„ M b ~>Pj B M E 



(A6) 



UJl 



arctan — , 



(A3) 



ct c i is the Coulomb scattering phase shift in the Z-partial 
wave, r\ c is the Coulomb parameter for the scattering of 
the fragments in the channel c. 



Here, (f> cs n m ,, is the channel spin wave function in the 
channel c s" m s » , Mi is the spin wave function of 
particle i, I" (mi») is the relative orbital angular mo- 
mentum (its projection) of the fragments in the channel 
c, r c = TbB is the radius-vector connecting b and the 
center-of-mass of B. We adopt the channel radius R c 
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large enough to neglect antisymmetrization between nu- 
cleons of b and B at r c = R c , that is 



-1/2 



r c =R c 



r c =R c 



(A7) 



where N c = (^)' 

Assuming that the overlap of the channel c at the chan- 
nel radius R c with the channel c' is negligible we get for 
the component of X^^'m ,, projected on £ c = ipbifiB at 
r c = RcM 



~J F M F / t-> ~ \ , > I y~JfMf . 



r c =R c 



(A8) 



where 



ires",/" J F (^ C ) = A^ c ^ Wrq «cs",i" J F j(r c ). (A9) 



At r c = i? c by definition [3£ 



-cs"i"^(-Rc) = V 2 ^c RcItcs" i" j f , (A10) 



where 7 TCS » /" ,/ F is the reduced width amplitude of the 
level r in the channel c s" I" Jf- I remind that the system 
of units H = c = 1 is being used throughout the paper if 
not specified otherwise. Then 



Z J r F J> F ms ,,(Rch) = ^ E V^cRcJtcs" l>> Jf 



iJ f M f 
c s" I" m a 



(All) 



Thus we can express the component Sj^^/^ „{ r c) 
taken at the channel radius r c — Rc in terms of the 
sum of the reduced width amplitudes, where the sum is 
taken over all allowed in the channel c partial waves I" 
at given Jp and s" . Then the component of ^cfm™ in 
the exit channel cs" m s n projected onto £ c = ifb^B at 
r r = R r takes the form 



1 csm 3 :cs" m,„\ rLcLc > ~ ^ <ic\^ csms S 



^ J F (mt)^ 
m s ;c s' 

2 7T / k r 



N 



kc Rc y /^c 



2tt 



E e" 15 - z< < sm s l mi \J F M F > Yf^fa,) £ [T l/csi j F (£ c )] 1/2 [A" 1 ]^ S^ fl ^ m> „ (R c f c ) 



M, i mi 



iV 



e _i < sm s lmi\J F M F >< s" m s " /" m/» | Jf M F > Y* mi (k c ) 



M, 1 1" mi m^r 

E cs I Jf i E c)] 1/2 [A" 1 ] ^ y^^K lr c s" I" J F Y l" m,„ (?c) 0c 

I 



(A12) 



IAT— 1 



Here, s is any channel spin value in the channel c = b+B The diagonal component, I" — I and s" — s, which is 

allowed by the spin and angular momentum conservation needed to determine the elastic scattering amplitude (see 
law. In particular, s" may coincide with s, that is s" — s. below) is 



T -J F (mt)(+) 
1 cslm s :cslm„„ K^dc) 



2?T Ik, 



c e~ iS °° i l 



k c R c y /Lt, 

FvcsUriEc)} 1 / 2 [A-% T y/2^ c R c lTCsl j F Y t mv , (r c ) 



< s m s lmi\J F M F > < s m s » Z m/» | Jf Mp > Y* mi (k c ) 



iV 



(A13) 



l/.T — l 



A similar consideration can be applied when we consider the contribution of the channel d s' ra s > into 
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x j f m f _ In thig c h anne i £ c , = p A an d 

^m., = E < s ' m *' l'mv\J F M F > 

(A14) 



6c' S 'm s , = E < JnM n J A M A \s'm s i > 

M n Ma 



(A15) 



Here, C < s ' m s , is the channel spin wave function in the 
channel d with the channel spin s' and its projection m s i , 
I' (mil) is the relative orbital angular momentum (its 
projection) of the fragments in the channel d, r c > = v n A 
is the radius- vector connecting n and the center-of-mass 
of A. We adopt the channel radius R c / large enough to 
neglect antisymmetrization between n and nucleons of A 
at r c i — R c i, that is 



A{w T c > j Cc'0/ s / / msl u c , a , v Jf j} 
-1/2 



(A16) 



where Av 



f (A+i)\ 



(A + l)- 1 ' 2 . 



Assuming that the overlap of the channel d at the 
channel radius R c i with the channel c is negligible we get 



for the component of A^,^^ ( projected onto £ c . 



fA 



at r c i = R c i 

^7 E ^i' F s^l' F m e , U r c> a> V J F (Rc>), 



r c ,=R c . 



(A17) 



where 



irc's'i' j F ( r c) = N rJ ^ w TC 'jUc' S 'i'j F j(r C ')- (A18) 



At TV = Rd 



Ur c> s> V J F (Rc< ) = \/2/V Rc' T 



r c' s' V J p i 



(A19) 



where /v = fi n A, 7tc's'1'J f is the reduced width am- 
plitude of the level r in the channel d s' I' J F . Then 

S r c'^ m., f c /) = — v/2 He R c > f TC > s'l'Jp 

c ' I' 

x<i> J /Ji> F m ,> ( A2 °) 

that is it can be expressed in terms of the sum of the 
reduced widths amplitudes in all allowed partial waves V 
in the channel d at given J F and s' . Then the component 
*caml!c't' L,, projected on £ c , = ^ at r c > = E c < takes 
the form 



T J F (int)(+) ,j3 ~ \ 
1 c s m s ; c' s' m,/ V -""<=' c ' ' 



2 7T / fc r 
fc n V fj,, 



N 



e- i4 c?i' < S m s ;m ; |J F M F > y; mi (k c ) £ [r„ C8ijF (£c)] 1/2 [A- 1 ]„s*f m8 ,(fl c ,v) 



M, i m, 

2tt /fc, 



y f-^c 



^ e lS c° i l < sm s lmi\J F M F > < s' m s * I' mi>\J F M F > 



M 1 1' mi m ( / 



JV 



X Y,* mi (k c ) ^ [T„,l, r (Bc)] 1/J [A V V 2 /V #c' 7r c' a' I' Jf *l' n»,, (?C ) « 

I 



'c s' m a / • 



(A21) 



z/,r— 1 



The component T^f™*^, F m ,{R C ' ? c ') is given by 



,J F (int)(+) 



s ;c' s' l' m s ,(Rc' *c') — 



2 7T fc, 



c e iS " S ' i l ^2 < sm s lmi\J F M F > < s' m s > I' mp\J F M F > 



M mi m v 



N 



X^* mi (k c ) ^ [T„ c , , J F (S c )] 1/2 [A- 1 ],, ^2 /V Rc? It c> s> V J F Yl> m ,, (%> ) • 
i/,r=l 



'c' s' m / • 



(A22) 
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2. External scattering wave function *t,B 



Now we proceed to the expression for the \E'i + ' ) in the 

external region, where r c > R c or r c ' > Rc' ■ In the 

external region the wave function ^csmi with fixed 

channel spin and its projection in the incident channel c 

can be written as 

v&(e»*)(+) _ \rj(°) _l vTr( ea =*)(+) ( \W\ 

csm s csm s 1 csm s ;r ) yji^tjy 

where the first term is the incident wave and the second 
term is the sum of the outgoing waves in all the open 
channels. The incident term is 

* M 0) = 4 7T £ c E i l <sm s lmi\J F M F > 

Jf Mf I mi Ta s ii 

Fi(k c , r c ) 



x < sm s n lmi\J F M F > Y i * mi (k c )e I 

xY lmi (t c )<j> csmB „, (A24) 

where the subscript c means that the incident wave is in 
the channel c. The sum over m s » is a formal because 

< s m s I mi \J F M F > < s m s " lmi\J F M F > 

Jf M f 

— Sm s m s n ■ (A25) 

Note that here we use the incident wave with the unit 
amplitude rather than with the unit flux density. The 

component ^{^.^j m „ , which corresponds to the exit 
channel c s I m s >> and fixed J F , projected on £ c reduces to 

M mi 

x < s m s Imi \ J F Mp > < s m s » I mi\Jp Mp > Y* m (k c ) 



x e 



^^%^Y lmi (r c )^ 



Now we take into account that 

e iu "Oi(k c , r c ) - e-^> Ii(k c , r c ) 



Fi(k c , r c ) 



2i 



(A26) 



(A27) 



Here, Oi(k c ,r c ) and Ii(k c ,r c ) are the Coulomb Jost sin- 
gular solution of the Schrodinger equation with outgoing 
and ingoing asymptotic behavior (we follow the defini- 
tions used in 39]): 

Oi{k c ,r c ) rc Z°° e^^"-^ in(2*.r„)-W2+*„o] j ( A28 ) 
and 

h{k c ,r c ) rc 3°° e -^ c - Vc M2 fce r c ) 1 -iV2+.co]. (A29) 



Then we can rewrite T 



J F (ext)(0) 
csl m s ;c s I m s 
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in the form 



T J F (ext)(0) , s _ . . 

1 cslm B ;cslm s „\ L c) — 1 , l 
r*>e 1 c 

x < s m s I mi \Jp Mp > < s m s n lmi\Jp Mp > 

M mi 

x F ; * mi (k c ) [l,{kc, r c ) - e l2 ^'Oi(k c , r c )] Y lrril (r c )0 csm 

(A30) 



Thus the incident wave is the pure Coulomb scattering 
wave function in the incident channel c. The second term 
in Eq. (|A23|) is given by the sum of the outgoing waves 
in the open channels |40l |: 



cs m s ; r 



2?r x - 

l T~ E 



^7 l 



E 



■i 



x < s m s ImAJp Mp > Y l 



I in 



(kc) [ 



Jf Mf 1 1 mi m ; - 



Jcc^ss^n 



s 



Jf 



c s l;c s I 
X ^c s m J ■ 



Oi(k c , r c ) < sm s lmj\J F M F > Yj (r 5 ) 

(A3I) 



Here, £g is the product of the bound state wave func- 
tions in the channel c - - - ' 



c, c 



S^fr , is the S'-matrix 

csl;csl 



element for transition csl — > csl. Note that we consider 
the outgoing waves in the channel with given total angu- 
lar momentum J F , initial channel spin s (its projection 
m s ) and final channel spin s (its projection m§). Since 
only two open channels are taken into account here, we 
will write explicitly the outgoing waves in both channels. 
First consider the elastic scattering, that is the outgoing 
channel c = c = b + B and the channel spin and orbital 
angular momentum coincide with the incident channel 
values, that is s = s and 1 = 1. The component of the 
outgoing clastic scattered wave (csl — > csl) is 



(ext)(+) 

csl m s ;cs I m" s 
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Sc 



E 

Jf Mp rni m l h 



x < s m s Ityii\Jf Mf > < s m s I mi" \ Jp Mf > 

x i^* mi (k c )[e i2 ^« - S J / al . csl ]Oi(k c , r c )Yi mi „(r c ) 

x cj> csm: . (A32) 



Hence, the projection of ^ , 

2tt 



(ext)(+) 
c s I m s : c s 



, „ on f c leads to 



«(ea!t)(+) 
csl m s ;c slm 



(rc) 



E 



Jf Mp mi m l n 

x < s m s lmi\JF Mf > < s m s I mi» \ Jf Mf > 
x z l Yf mi (k c ) [e 42 ^ - S J J lcsl ] Oi(k c , r c ) Y lnii „ (f c ) 
X (t>c sm- (A33) 



Correspondingly, for the inelastic scattering, c — c but 
either s ^ s or / ^ I or both differ from the entry values, 
we get 



(ext)(+) 

csl m s ; c m H 



■ 2lT £ 

= ~ Z I ^ c 



E 

Jf Mp mi m^, 



x < Imi s m s | Mp > < l N mi>> s" m s >> \ Jf Mp > 
x i l Y?fc)S%, r c )Y l//mi „(r c )4> cs "n 



(A34) 
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Then the projection of ^cslm^cs' 1 I" 



on £ c is 



<<;<< m ,„( r c) =-* 



2tt 

for. Tr 



E 



*(ext)(+) 
'csl m s \c 8' 

Jf Mf Tni m;/, 

x < I mi s m s \ Jf Mf > < I" mi>> s m s *i\ Jf Mf > 

x i l Y l '$ c )S%, l „. c . l O l »(k c , r c )Y lflmil/ (r c )^ cs , ms/l . 

(A35) 

Finally, for the outgoing scattered wave in the reaction 
channel c=</ = ji + 4we have 



(ext)(+) 

c s I m s ;c' s' I' m s 



l yy~ l \hv~^ c ' E 

n c I c' V L c 7 » r 



x < I mi s m s | Jf Mf > <l' my s' m a > \ Jf Mp > 



(A36) 



It leads to its projection on £ c 



,(e:rt)(+) 



,(rc) 



2 7T l~V 



'csl m s : c' s' V 

Jf Mf mi Tn>i> 

x < I mi s m s | Jp Mf > <V rrii> s' m s ' \ Jf Mf > 

x Y l * mi (k c )S J c r s , l , csl O l ,(k c ,, r c ,)Y Vmi ,{i c >)4> c > s > me ,. 

(A37) 

Now we can derive the expression for the matrix ele- 
ments of the 5* matrix. Since the wave function \&c is 
continuous using Eqs. (|A13|) . (|A30[) and (|A33|) we get 
the equality 



r^J F {mt) , p - ^ _ ~ 

L csl m.icslm.,, K-^d-c) — J- 



J F (ext)(0) ( f? r } 

cslm s ; cslm,„ kP-c l c) 



T 



(ext)(+) 



- cslm s ;cslm s /i Fc )' 

which boils down to 

N 



(A38) 



-I 



V.T—1 



= i [h{k c , R c ) - .s/iv. ^o, .:/.-,. R c )] . 

Taking into account that [39| 

Ijjkc, Rc) = Gi{k c , R c )~iFi{k Cl R c ) l2uJci = 
Oi{k c , R c ) G,(fc c , R c )+iFi(k c , R c ) 6 

and 



(A39) 



^tcsUf( E c) - 2Pcl(E c , R c )1tcsUf' 



where 



P cl (E c , R c )- ^ Rc) C +Gf(k c , R c ) 



(A40) 
(A41) 

(A42) 



is the Coulomb-centrifugal barrier penetrability, we get 
the elastic scattering matrix S- matrix element: 

AT 

S^ ;csi = e - 2 ' i5 '"(l+* E [r,c S UF(Ec)} 1/2 [A-\ T 

ZV,T — 1 

x [TrcsuAE^] 1 ' 2 ). (A43) 
From equality of Eqs. (|A22I) and (IA37|) at r d = R c > 

yJ F (int) (r> £ \ _ y(ext)(+) (r> f. \ 

1 csm s -c' s' m s /V- rt C 1 c' ) — csl tci s\c' s' V Tfi s i \^^c ! *■ c r ) 

(A44) 

we obtain the reaction matrix element: 

N 



SZ^s^ie-iKe-*** E [^csIJf(Ec)} 1 



/2 



I/,T— 1 

x [A-V [T rc , a , VJp {E c ,)) 1 / 2 . 



(A45) 



Both obtained matrix elements coincide with the corre- 
sponding matrix elements from [40(. The only difference 
is in the definition of the solid scattering phase shifts. 
The obtained matrix elements of the S matrix confirm 
that the relative normalization of the internal and exter- 
nal wave parts of are correct and we can use them 
to calculate the reaction amplitude of the deuteron strip- 
ping proceeding through resonance states. 



Appendix B: Matrix element Nig 

Let us consider the DWBA surface (in the subspace 
over v n A) matrix element 

M° W {P, k dA ) =< X « T£f <">| ^ - 2* | 



Xi PdXdA> 

r nA >R n 

mF^ f) (P, *dA) + Ms£a)Q<pF, k dA ) 



where T' nA =< ipA^^n >, 



(Bl) 



bB 



M™ F) (P, k dA ) = J dv nA J dr.Fxit T ( :f ( ->*(r nA ) 

r nA >R n A 

x[ j f p F-f pF ]ipdX$> (B2) 
and 

M™ A) (P, k dA ) = J dr pF J dr nA Xp -> T^->(v nA ) 

r nA >R n A 

x[^ nA -f nA ] VdXd + J>. (B3) 



^s{pF) can ^ e wr itten as 
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Mg { J F) (P,k dA )= Jdr nA J dr pF #>?tt*( T nA)[K F -?pF]<PdX ( & 
- J dr nA J dr pFX £* T^->(r nA ) [f pF - f pF ]tp d x <# > 

r nA <R nA 

= J dv nA J dv pFXp -J* T^->(r nA )[f pF ^ pF ]<fdX d + 2- (B4) 



We took into account that for any finite volume r nA < This matrix element can be rewritten in the form, in 
R nA the matrix element containing T pF - ~f ' pF vanishes which the integral over r pF is transformed to the surface 
as it has been discussed in Section II A for deuteron strip- integral: 
ping to bound states. To estimate Mg™ F ^ we need equa- 
tions connecting different variables: 

r dA = 1/2 v pn + r nA , (B5) 
r pF = A/(A + l)r nA +r pn . (B6) 



Now in the matrix element (|B4p we replace the variable 
*nA by r pn . Then we get 



m dw 

M S(pF) 



(P, k dA ) = -(^) 3 J dr pn J dr pFX ( -/(r pF ) 

x r^->(^[r pF - v pn ]) [T pF ^ pF ] 

. , (+) .A+l A + 2 , 
x MrpnlXdA t p f 2A~ Tpn >- ( 7 ) 



M 



DW 



A + l 3 _4_ 

1 A ' A + l a-^ioo 1LpF 



x [ dv pn Lp d (r pn ) I dn TpP [XdVi 



lim Rt 



2 /ipp 



M ,A+1 A + 2 . dx { pF ] *{v pF ) T^->{^[r pF - r pn }) 



A 



■ r pF 



2A 



(+)(A+1 



X pF (r pF )T 



A+2 



A P F 1A P n > 



ii A 



A 



dr pF 



'r pF 



r pF =R pF ^oo 



(B8) 



Due to the presence of the deuteron bound state wave 
function the integration over r pn is limited. At r pF — ¥ 
co and r pn < oo we can replace the distorted waves in 
the initial and final channels by their leading asymptotic 
terms: 



( + ) I \ r dA~ 

XdA ( r dA) ~ 



a ik dA r dA +ir) dA \n(k dA r dA -'k dA -r dA ) 



(B9) 



and 



(-)* r pF ~>oc iU. pF -r pF +ir] pF ln(k pF r pF +k pF -r pF ) 



X pF 



(BIO) 



Here, rjij is the Coulomb parameter of particles i and j 
in the continuum. Note that r dA — ^j^- r pF — r pn, 
and at r pF — > oo and r pn < oo r dA — ¥ oo. Then 
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Q e ik dA -r dA +ir] dA ln(k dA r dA -k dA -r dA ) 

dr pF 



r pF ^oo . A+l ^ _ eikdA .(A+i .r pF -4^r J> „)+i )7( i J 4 ln(fc dA r dA -k dA -r dA ) 



.4 



(Bll) 



and 

^g — ikpF'fpF+llIpF ln(fe p F fpF + kpFTpp) 



(B12) 



For "f^*^ _ ^*(^x^[rj,^— r pn ]) we can take only the exter- 
nal part, which contains the resonant S matrix element, 
see Eq. (|A37I) . Neglecting all the spin-dependent and 
angular parts and leaving only its radial part, which is 
O n A (f n A )/r n A j we get for its leading asymptotic term: 



The leading term of its derivative at r p p — > oo is 



dOnAjknA, r nA ) /r n A r pF ^oo . 1 
— t I K n A 

or P F r pF 



x e 



x e 



-i[itn A ln(2 ^1 fe„ A r pF )+l nA tt/2- cr„ A0 ] 



(B14) 



OnA(fc«A, »*7iA) rpp-4-oo ^4 1 f a+i 



X e 



r„A A + 1 r pF 

-i[r)n A \n(2k nA r nA )+l nA -x/2—cr nA a] 



(k„A r pF -k nA r pF -r pn ) 

(B13) 

I 



Then Mg^ F ^ reduces to 



M S(%F) (P, kdA) ~ lim R pF [ dr pn (p d {r pn ) e+ l ^ * dA ' rpn [ dn rpF [(— ^ k <M + k pF ) • r pF 

RpF— >oo J J A 



A + l 
A 



i (^Ji k dA -k pF )-f pF _R pF t fc„ A -RpF g'fdA In (fc<i A r dA (R pF )-k. dA -y dA (R p f)) +i VpF ln(fc pF flp F +kp F -fpF R pF ) 



g" v A — aj\ — pr l -pr ---pr g - A •"/isi - ~p r g 

x g-*[ ^ fc„A fpFT p „+l)„A ln(2 ^±1 fe„ A _R p f)] 



(B15) 



Taking into account the asymptotic behavior of the plane 
wave 



e iq-r pF rpRjao 27T 



iqr pF 
e -iir pF S (q + r pF )] 



[e*« r "5(q- V) 



(B16) 



where q = A ^ ^n dA — k pF we obtain that the matrix 



element 



M^ F) (P, k dA ) ~ lim AC^) e»«^* 



h(RpF) e 



— t 9 -RpF] 



Thus Mg^fpj has no limit at R p f — > oo but regularization 
of this matrix element by integrating the matrix element 
over an infinitesimal bin in the momentum plane leads to 



disappearance of M DM 



S(pF)- 



9+e 



Te I dqM^ F) (P,k dA ) 



lim 



sin (e R p f) 



R pF —>co e R p F 
q-e 

[e iqR * F fi(R pF ) - e~ iqR * F f2{R pF )} = 0, (B18) 



where e << q. 

A similar prove can be applied to estimate 
given by Eq. (|B3|) . Since the integral over r n A is taken 
over external volume with r n A > R n A the transformation 
of the volume integral into the surface one leads to two 
surface integrals: 



(B17) M^ A) (P, k dA ) = -M^ a (P, k dA ) + Mg» (P, k dA ). 

(B19) 



The first term is the surface integral encircling the in- 
ner surface of the external volume at r n A = R n A: while 
the second term is the surface integral taken at r n A = 



34 



R A — > oo. A negative sign in front of the first term 
appears because the normal to the surface is directed in- 
ward to the center of the volume, that is opposite to the 



normal to the external surface (at infinitely large radius). 
The surface integral over the infinitely large sphere in the 
subspace over r nA is 



M^(P, k, 



dA 



, 2 1 

, llm RnA 7T. 

-R . — f oo ^ fJ-nA 



dr pn (pd(r P n) [XduiO-/ 2 ) 

(+) 



dy$\A/(A + 1) r nA + v pn ) T^f^jr^) 



dr n/ 



X { p t(A/{A + 1) r nA + r pn ) T^->*(r nA ) **** ^ + ^ 



(B20) 



Here, the Jacobian variable r pF is replaced by r pn by 
r P n = r P F — A/(^4 + 1) r n A- The disappearance of the 
matrix element (|B20|) can be proved similarly to the proof 
for Mgfopy Due to the presence of the bound state 
fd(r pn ) the integration over r pn is limited by finite dis- 
tances. Hence, r pF — > oo and r dA — > oo at r nA — » oo. 
Replacing the distorted waves by their leading asymp- 
totic terms (IB9[) and (|B10I) . singling out the plane wave 
containing v nA and using the asymptotic representation 
of this plane wave, see Eq. (|B16p . integrating over ft rnA 
we eventually arrive at 

M™(P,k dA )~ Hm [e^^ gi (R nA ) 

R nA -+oo 

+ e -WR nAg ^R' nA )]. (B21) 



Regularization of this matrix element by integrating it 
over an infinitesimal bin in the momentum plane q' leads 



to disappearance of MP W , that is 



M™J P ,k dA ) = -M™(P,k, 



Appendix C: Matrix element M s 



Id A 



CDCC(post) 



(B22) 



(P, k dA ) 



Here we show how to transform M 



CDCC(post) 



into the 



surface integral over the coordinate r nA . ]\/[C DCC (P° st ) 
can be written as 



Mg DCC (P,k dA )= J dr nA J dr pFX { ^ pF (r pF )T 

r n A>R n A 

= Mg° t cc (P, k dA ) - Mg°f c (P, k dA ), 



nA (Y nA )[i - i \w i (r pF , r„ A j 



(CI) 



where 



M 



CDCC 
Stot 



(P, k dA )= J dr nA J dr pF x { ^ pF (r pF ) 



Ti-/(r„ A )[^ - f}^f DCC ^\r pF , v nA ) (C2) 



and 



M 



CDCC 
Sint 



(P,k dA )= J dr nA f dr pF x { -k pF (r P F) 



r nA <R n 



xT<£ t)( ->*(r BA )[?T-2*]* 



cdcc(+) 



(r pF) r„ A ). (C3) 



Note that in the matrix element the integration is 

carried over r p j? and r nA in all the coordinate space while 



in Mg^ t the external region in the subspace over r nA is 
excluded. Let us first consider Mg^. The CDCC wave 
function is given by Eq. (f5*2"j). If we substitute the first 
term, n = 0, which contains the deuteron bound state 
wave function, the transformation leads to the surface 
integrals with r pF = R pF — > oo and r nA — R nA — > oo. 
Both surface integrals vanish and the proof is similar to 
the one presented in the previous section. For the rest of 
the CDCC wave function corresponding to the sum with 
n > 0, which we call ^ff CC< " + \ transformation to the 
surface integrals gives 
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M™ [C \P, k<M) = J dr nA J dr^ X W F (r P F)Ti-/(r IlA )[T pF -^ pF ]^ cc W(r pF , v nA ) 
+ J dr nA J dr pF X { +l F (r pF )r ( -2*(r nA ){T nA -T nA }* 



CDCC(+), , 



Hit) —^E- I dQ r I <lr„,j 



RpF — >oo 2 /ipF 



Y (+) (r „\r(->(r a) — 



CDCC{+)( s 



dr P F 



>r pF 



+ lim R " A 



di} TrlA I dr 



pF 



v (+) /_ \ T (-W- ^ ( r pF, r wA ) 



v (+) (r v)M> CDCC(+) (r „ r ,,1 9T nA*( r ^) 



(CM) 



Let us, first, consider the first term, in which R p f — > oo. 
Let us divide the integration region over r Mj 4 into the 
region r nA / R pF — > and the region where r nA > R pF — > 
oo. In the first region we get that ~ R p f — > oo and 
r pn ~ i? pF —¥ oo. Taking into account the asymptotic 

behavior of ^f c DCC(+) {r pF , v nA ) ~ r"j and Eq. ([B16| 

we get that the first term goes to zero as R~p — > 0. In 
the remained region r nj 4 ~ R pF — > oo and we consider it 
later. The second term of Eq. (|C4| . in which R n A — > oo, 
we also separate into two regions: r p p/R nA — > and 
r pF ^ -RnA — > oo. In the first region r pn ~ i? Mj 4 — > oo 



0. To 



and r^A ~ -RnA -> oo and ^> ic ( \r pF , r nA ) 
Hence the matrix element goes to zero as R~ A — 
consider the behavior of the first and second terms of 
Eq. (|C4[) in the second regions, where r nA , r pF — » oo, 



it is more convenient to introduce the hyper-spherical 
coordinates in the six-dimensional hyper-space: 



VnA 2 , VpF 2 

' nA ' oF' 



i\i A = p\ sin a, 

PnA 

< a < n/2. 



r P F = P \ cos a, 

VpF 



(C5) 



Here, m is the scaling mass parameter, for example, 
the nucleon mass. Then Mg^ t in the region, where 
r n A, r P F — > oo, can be written as the integral over the hy- 
persphere encircling the volume integral with the radius 
of the hyper-sphere p — > oo (38[ : 



m 



1 

2 {pnAP P FY /2 p"^°° 

d_ 

dp 



lim p 5 / dr 



7r/2 

r 

dr nA / da sin 2 a cos a 
o 

d 



X p f V c pF) l nA V CnA )~n^ ! i [r P F, r nA ) -W i [YpF, r nA ) — x pF (r pF )l nA (r r , 



(C6) 



Here the hyper-radius p is the parameter going to infin- 
ity. The integrand is contains highly oscillating (actually 
infinitely oscillating) functions. The behavior of the inte- 
gral at p — > oo depends on the asymptotic behavior of the 
integrand. The integration over dr pF can be performed 

directly using the asymptotic form of X p ~f* ( v pf)- It i s 
given by the Coulomb distorted plane wave, but for sim- 



plicity, what does not affect the final result we neglect, 
as in the previous section, the Coulomb effects. Then 
the asymptotic form of the plane wave is given by Eq. 
(IB16[) and, hence, integration over r pF using (5-functions 
is trivial leading to r pF — ±h pF . After performing the 
integration over dr pF only two integrals are left. From 
Eqs (|B5]), p6]) and JC5]) we get for 



3G 



.4 2 



(A + l) 



2 T nA 



I m „ A j m 

P \ / cos a T — 7 \ — 1 : 

V PpF A+l y H pF V PnA 



m A 2 
z sin 2 a + 



(A + l) 2 fin. 



sin a 



(C7) 



and 



4 r P F + 



A 



2(A+l 



TpF T n A 



(A + 2) 2 
4(A + 1) 2 



.4 



/l m 

V 4 ^ 4(4 + 1) 



m 



z sin 2 a 



(A + 2) 2 to 



4(4 + l) 2 pL n A 



sm a . 



(C8) 



Here, z = Y nA -\t p F- We recall also that in Eq. (|52[) at n > 

ippn{ v pn) at r p „ — > oo contains the asymptotic terms 

±ik pn r pn («)(+)/ \ e » fc dA '"dA , p 

s — p- , while x\ ijdA) ~ £ — ^ — , where we, for 

simplicity, neglected the Coulomb distortion. Then after 
integration over df p p the leading asymptotic form of the 
integrand with omitted Coulomb effects is a product of 
highly oscillating at p — > oo exponents: 

p ±ik pF r pF p ik n Ar n A p±ikp n r pn pikdATdA 



p->oo J_ ipg(a, Z ) 



(C9) 



Thus we need to estimate a highly oscillatory integral: 

1 

Ji~ lim / dz / dasin 2 a cos 2 a e lp9(z ' a) . (CIO) 
-i o 

Evidently that this integral and, hence, Mg^.(P, k^) 
vanishes at p — > oo, whether a stationary phase point 
does exist or not, because the integration brings p to the 
denominator. 

Now we proceed to Mg^ t (P, k^). We rewrite is as 

MgZiP, = j dv nA J dVpFX^M 

r nA <R n A 

T (mt)(-),, , rfc 7^ I ff CBCC(+) f ^ 

x I nA (r„A) [1 p f - 1 P F\ *i [r pF , r nA ) 

J dr nA J &YpFXpF ] * \ v pf) 



r„A<Rr, 



CDCC( + ) 



(r pF , r nA ). 
(Cll) 



Let first consider the first matrix element containing T p p ■ 
It is easy to show that this matrix element vanishes. Af- 
ter transforming it into the surface integral over r p p we 
get 



J dr nA J dr p FX ( p F ) *( r pF)T { 1 l lA t)i h (r n A) 



r n A<RnA 



J pF — J- pF 
1 



, Tl CDCC( + ), \ 



2 ppF RpF 



lim i? 2 f / dO rpJ 



CDCC(+) 



(r pF , r nA ) 



/ dr„,Ti-/(r 

r-„A<H„A 



dr P F 



XpF*(rpF) 



^f DCC(+) (r pf ,w) 
dr pF 



r v p —R v 



(C12) 



The matrix element containing n = term of the CDCC 
wave function vanishes because in the subspace r nA < 
R n A at r p i? — » oo the deuteron bound state wave func- 
tion exponentially fades away. The terms of the CDCC 
wave function with n > 1 also produce vanishing matrix 
element because the CDCC wave function corresponding 
to these terms in the subspace r nA < R nA at r p p oo 
decays as l/r^p, that the whole matrix element vanishes 

as lim i? 2 p/i??, p — > 0. Thus we arrive at 
RpF^oo pr pe 
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Mg w (P,k dA ) = -M!?Z(P^<lA) 



dr. 



'pF 



r n A<RnA 



' UnA 



CDCC(+) 



(r pF , r nA ) — — T^*(r nA ) - T y nJ >*{r nA ) — — 
or nA or nA 



CDCC(+) i 



TpF, YnA) 

(C13)' 
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There are two main reasons for absence of the practical theory of stripping to resonance states 
which could be used by experimental groups: numerical problem of the convergence of the DWBA 
matrix element when the full transition operator is included and it is unclear what spectroscopic 
information can be extracted from the analysis of transfer reactions populating the resonance states. 
The purpose of this paper is to address both questions. The theory of the deuteron stripping is 
developed, which is based on the post continuum discretized coupled channels (CDCC) formalism 
going beyond of the DWBA and surface integral formulation of the reaction theory [A. S. Kadyrov 
et al., Ann. Phys. 324, 1516 (2009)]. First, the formalism is developed for the DWBA and 
then extended to the CDCC formalism, which is ultimate goal of this work. The CDCC wave 
function takes into account not only the initial elastic d + A channel but also its coupling to the 
deuteron breakup channel p + n + A missing in the DWBA. Stripping to both bound states and 
resonances are included. The convergence problem for stripping to resonance states is solved in 
the post CDCC formalism. The reaction amplitude is parametrized in terms of the reduced width 
amplitudes (ANCs), inverse level matrix, boundary condition and channel radius, that is the same 
parameters which are used in the conventional i?-matrix method. For stripping to resonance states 
many-level, one and two-channel cases are considered. The theory provides a consistent tool to 
analyze both binary resonant reactions and deuteron stripping in terms of the same parameters. 
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I. INTRODUCTION 



Production of unstable nuclei close to proton and neu- 
tron drip lines has become possible in recent years, mak- 
ing deuteron stripping reactions (d, p) and {d, n) on these 
nuclei (in inverse kinematics) not only more and more 
feasible as beam intensity increasing but also a unique 
tool to study unstable nuclei and astrophysical (n, 7), 
(p, 7) and (p, a) processes. The deuteron stripping re- 
actions populating resonance states of final nuclei are 
important and most challenging part of reactions on un- 
statble nuclei. If for nucleon transfer reactions popu- 
lating bound states for about fifty years experimental- 
ists used the standard distorted waves Born approxima- 
tion (DWBA), an adequate theory for transfer reactions 
to resonance states yet to be developed. By standard 
DWBA I mean the approach in which the one-step trans- 
fer matrix element is evaluated with incoming and out- 
going distorted waves calculated by fitting the deuteron 
and proton elastic scattering with local optical potentials. 
The transition operator contains finite range effects as 
well as the full complex remnant term. The main idea 
of the DWBA is that the transition matrix element is so 
small that one can use the first order perturbation theory. 
Since the nuclear potential is quite large by itself (~ 100 
MeV), the smallness of the transition operator can be 
fulfilled only if the reaction is peripheral enough, so that 
the non-diagonal matrix element, representing the trans- 
fer reaction amplitude, becomes small. However, since 
the resonance wave function is large in the nuclear inte- 
rior and different channels are coupled in the nuclear in- 
terior, the character of the stripping to resonances can be 



quite different from the stripping to bound states. Nowa- 
days the standard DWBA is gradually being replaced 
by more advanced approaches like continuum dicretized 
coupled channels (CDCC) adiabatic distorted wave 

(ADWA) [1], coupled reaction channels (CRC) and the 
coupled channels in Born approximation (CCBA) avail- 
able in FRESCO code There are two main reasons for 
absence of the practical theory of stripping to resonance 
states which could be used by experimental groups. First 
one is the numerical problem of the convergence of the 
DWBA matrix element when the full transition operator 
is included. However, it is only a technical problem. The 
second pure scientific unsolved problem is what spectro- 
scopic information can be extracted from the analysis of 
transfer reactions populating the resonance states. Be- 
sides, since the standard DWBA is deficient to more ad- 
vanced methods like CDCC or ADWA, a new approach 
should go beyond of the DWBA. 

Majority of theoretical works devoted to the develop- 
ment of the theory of single-nucleon stripping into un- 
bound states of the residual nucleus have been published 
in 1970-s [6U2l|. Great interest in these reactions at that 
time stemmed primarily from the fact that they allow one 
to extract reliable information on the properties of nu- 
clear resonant states by means of the combined analysis 
of the data on stripping and elastic resonant scattering 
of nuclcons from the target nucleus || [TH, 
most of the cited works the theory of stripping into reso- 
nant states was developed within the standard DWBA 
by analogy with usual stripping to bound states. In 
this case the expression for the reaction amplitude ob- 
tained instead of the bound-state wave function for the 
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captured nucleon (form factor) contained a continuum 
wave function which leads to slow convergence of the ra- 
dial integrals or even to their divergence depending on 
the choice of this wave function. In Refs. [a ©, [HI the 
form factor was taken to be a scattering wave function, 
which described the resonant scattering of the nucleon 
from the target nucleus. This wave function was calcu- 
lated using a single-particle potential whose parameters 
were adjusted to give a resonance with the correspond- 
ing properties. The Gamov decaying-state wave function 
and the Weinberg wave function which are regular at 
the origin and purely outgoing at infinity were used in 
Refs. 53 and [14|, respectively. Various methods were 
suggested to calculate radial integrals practically with 
the above-mentioned form factors: (i) the introduction of 
the convergence factor exp(-ar) into the integrand [6|; 
the integral obtained was calculated for various a > 
and then its values were extrapolated numerically to the 
limit of a = 0; (ii) the method of contour integration 
in the complex r-plane (complex scaling) [3; (hi) the 
method based on the correct account of the boundary 
conditions in the three-body scattering problem [ll| ; (iv) 
the Zeldowich-Berggren method [13 of the regularization 
of integrals containing the Gamov function in which the 
convergence factor exp(— ar 2 ) was introduced [13; (v) 
the pseudo-bound-states method [T3|. The methods (ii) 
and (hi) were most convenient for numerical calculations. 
Although the above methods allow one to avoid formal 
difficulties, nevertheless all the methods are rather com- 
plicated because of cumbersome numerical calculations 
and carry on the shortcomings of the standard DWBA 
for stripping to bound states. 

Even if we put aside the technical problem of conver- 
gence of the matrix element for stripping to resonance 
states, there is more important question remains: the 
spectroscopic information which can be extracted from 
analysis of deuteron stripping reactions (and other trans- 
fer reactions) into resonant states. This is really a crucial 
question because the answer determines the reason why 
we measure nuclear reactions. For more than 50 years 
transfer reactions to bound states, and deuteron strip- 
ping in particular, have been used to determine the spec- 
troscopic factors, which measure the weight of the single- 
particle state in the overlap function of the initial and fi- 
nal nuclei. That is why there was always a temptation to 
develop a theory of stripping into resonant states which is 
fully similar to stripping to bound states. For example, in 
13] it was assumed that the spectroscopic factor could be 
extracted from deuteron stripping into resonance states. 
In this case the spectroscopic factor is the ratio of the 
observable and single-particle resonance widths. How- 
ever, the spectroscopic factor is not observable and de- 
pends on the single-particle potential used to calculate 
the single-particle width. In 22] it has been shown that 
spectroscopic factors are not invariant under finite-range 
unitary transformations and, hence, in exact approach 
nuclear reactions cannot be a tool to determine spectro- 
scopic factors. In [22| it was called separation of nuclear 



reactions and spectroscopic factors. However, there is a 
model-independent information, which can be extracted 
from deuteron stripping reactions. I mean the asymp- 
totic normalization coefficients (ANCs), which are the 
amplitudes of the tails of the overlap functions [13 and 
are invariant under finite range unitary transformations. 
The most model-independent definition of the ANC is 
that it determines the residue of the elastic scattering 
S matrix in the pole corresponding to bound, virtual or 
resonance states. For the resonance state the ANC and 
partial resonance widths are related [2J, [25J : 

[c F Ajl ] 2 = (-r/e 2 ^ fc ^«> ^-r xAjl . (i) 

fe a;A(0)j'i 

Here I and j are the orbital and total angular momentum 
of particle x in the resonance state F = (Ax), fi xA is the 
reduced mass of x and A, k x A(o)j i is the real part of the 
resonance relative momentum of x and A, <fij i(k xA ) is the 
non-resonant scattering phase shift, and T xA ji are 

the ANC and partial resonance width in the channel x+A 
with the quantum numbers I and j. Eq. ([!]) stands for 
narrow resonance, i.e. for k xA mji << k x A(o)ju where 
k x A(i)j i is the imaginary part of the resonance momen- 
tum k xA (R)ji = k xA (o)jl ~ik xA (i)ji, which determines 
the location of the resonance pole in the momentum 
plane. Due to relation ([T]), the resonance width is also 
invariant under finite-range unitary transformations and 
can be determined from the experiment. 

Nowadays, it is quite well understood that the ANCs 
can be determined from peripheral transfer reactions, 
see [26l - [3~fl ] and references therein. However, the ANC 
method has been applied only for transfer reactions pop- 
ulating bound states. It is well known that from bi- 
nary resonance scattering and reactions using the con- 
ventional i?-matrix approach one can determine the res- 
onance partial widths, which, as we have underscored are 
related to the ANCs. i?-matrix method is one of the most 
popular tools among the experimental groups worldwide 
because the approach is comparatively simple even for 
many-body, many-channel cases and deals with the for- 
mal partial resonance widths determined from the fit to 
the experimental data. These formal widths can be eas- 
ily related with the observable partial widths. Using the 
i?-matrix approach one can fit simultaneously data for 
all available channels. It allows one to control the consis- 
tency of the obtained physical parameters. The question 
is whether the theory of stripping to resonance states can 
be formulated in terms of the same parameters which are 
used in the i?-matrix analysis of the binary resonance re- 
actions. 

It is the purpose of this paper to deliver a theory of the 
deuteron stripping, which will solve all the above men- 
tioned problems for the deuteron stripping into resonant 
states. This theory is based on the post CDCC formalism 
going beyond of the DWBA and surface integral formula- 
tion of the reaction theory (32| . The CDCC wave function 
takes into account not only the initial elastic d + A chan- 
nel but also its coupling to the deuteron breakup chan- 
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nel p + n + A missing in the DWBA. The convergence 
problem is also resolved in this formalism. The reaction 
amplitude is parametrized in terms of the reduced width 
amplitudes (ANCs), inverse level matrix, boundary con- 
dition and channel radius, that is the same parameters 
which are used in the i?-matrix method. Thus the the- 
ory provides a consistent tool to analyze both binary res- 
onant reactions and deuteron stripping in terms of the 
same parameters. 

The theory is based on the surface-integral formula- 
tion of nuclear reactions and valid for stripping to both 
bound and resonance states. First, just for demonstra- 
tion of the formalism, the transformation of the DWBA 
amplitude for stripping to the bound state is presented. 
The reaction matrix element is split into two parts: inter- 
nal (over the relative coordinate between the transferred 
nucleon and target) and external. The idea of such sep- 
aration is based on the fact that in the post formalism 
the main contribution to the stripping amplitude comes 
from the nuclear exterior while the prior form amplitude 
is dominated by the internal region. It will be shown 
that the dominant external post (internal prior) ampli- 
tude using the Green's theorem can be written as the 
dominant surface integral encircling the internal volume 
plus small addition from the prior external (post internal) 
part. Thus, both post and prior forms lead to the same 
reaction amplitude given by the sum of small internal 
post form, small external prior form and the dominant 
surface integral. The contribution of the post internal 
part can be minimized by a proper choice of the final- 
state optical potential, and the other two amplitudes are 
parameterized in terms of the reduced widths amplitudes 
(ANCs). After that the theory is extended to the CDCC 
formalism. Then the theory is applied for stripping to 
resonance states. First it is developed for the standard 
DWBA and then the post CDCC formalism based on the 
surface integrals is developed. One of the most impor- 
tant results of this paper is that the post CDCC form 
for stripping into resonant states can be written as the 
sum of the small internal (over the coordinate r n A) post 
form and the dominant surface part. The absence of the 
diverging (or poor converging) external part solves the 
problem of convergence of the matrix element for strip- 
ping to resonance state. 

In the developed approach the information about the 
resonance subprocess is contained in the scattering wave 
function of the fragments formed by resonance decay. 
This wave function is written in a standard i?-matrix 
form using its separation into the internal and external 
parts. It allows us to generalize the i?-matrix method 
for binary reactions to stripping reactions. Since the 
deuteron stripping into resonant states is 2 — > 3 parti- 
cles reaction, the excitation of the resonance occurs in 
the subsystem, while the third particle causes the distor- 
tion. The extracted partial resonance widths can be used 
for calculation of the (n, 7) processes. If the cross section 
for (n, 7) resonant capture is available, the simultaneous 
fit to the deuteron stripping and (n, 7) resonance cap- 



ture can be done. The method can be also applied for 
analysis of the Trojan Horse reactions [33]. Concrete cal- 
culations and the application of the theory for deuteron 
stripping and Trojan Horse reactions will be presented 
in the following up papers. In what follows we use the 
system of units in which H = c = 1. We also neglect the 
spins of the particles if not specified otherwise. 

II. SURFACE INTEGRAL FORMULATION FOR 
DEUTERON STRIPPING TO BOUND STATE. 

Before the theory of the deuteron stripping to reso- 
nant states will be outlined I will present a surface in- 
tegral formulation of the theory for stripping populating 
bound states. First, just for demonstration, I consider 
the DWBA and then extend it by including the CDCC 
wave functions. As it has been explained in Introduction, 
the transfer reaction matrix element will be split into two 
parts in the subspace determining the relative motion of 
the transferred nucleon and target: internal and external 
parts. After that replacing the potentials in the transi- 
tion operators by the kinetic energy operators and using 
the Green's theorem the matrix element in terms of the 
surface integral will be obtained. 

A. Stripping to bound state. Post form of DWBA. 

In this section we consider the post form DWBA am- 
plitude, which we split into the internal and external part 
in the subspace over the relative coordinate between the 
transferred n and A. Due to the choice of the transition 
operator in the post form, the internal part turns out to 
be small. The external part, which is parameterized in 
terms of the ANC, will be transformed into the dominant 
surface integral encircling the internal volume and small 
external prior DWBA amplitude. 

We start consideration from the exact reaction ampli- 
tude for the deuteron stripping to bound states 

d + A^p + F, (2) 

where F = (An) is the bound state. The post form of 
the exact reaction amplitude 

M^ os *)(k pF , k dA ) =< $(- } | AV pF |*| +) >, (3) 

where is the exact scattering wave function in the 

initial state with the two-body incident wave d + A, 
&f ^ = XpF V*f is the channel function in the exit state 
p + F, (fi is the bound-state wave function of nucleus 
i, x\j = Xk + ^( r y) ^ s tne distorted wave describing the 
relative motion of particles i and j with the relative mo- 
mentum k^; A V p f — V p a + V pn — UpF is the transition 
operator in the post form, is the microscopic interac- 
tion potential between nuclei i and j, Uij is the optical 
potential between nuclei i and j; r,y is the radius- vector 
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connecting the center of mass of particles i and j. I 
remind that the exact wave function is fully anti- 

symmetrized but the channel wave function $^ ' is not 
antisymmetrized with respect to exchange of the exiting 
proton and nucleons in F. However, the internal wave 
function of F ipp in $^ ^ is fully antisymmetrized. The 
reason why we can drop the antisymmetrization in the 
channel wave function is the presence of the fully anti- 
symmetrized exact wave function in the initial state and 
fully symmetric transition operator what can be seen be- 
low when the transition operator is expressed in terms of 
the kinetic energy operators. 

To obtain the post form of the DWBA from Eq. © 

we replace by the channel wave function = 

,(+) 



'■Pd <Pa XdA m the initial d + A state: 



Then wc use approximation 

tfF ~ I A <PA, 



s(+) 



> 



(4) 



(5) 



where 1% (r ra ^) is the overlap function of the bound state 
wave functions of nuclei F and A: 



Ia^ua) = (<Pa\ <Pf)- 



(6) 



Note that the integration in Eq. ^ is taken over all the 
internal coordinates of nucleus A. Then the transition 
operator in Eq. @ takes the form < <pa\& V pF \(fiA > — 
< <Pa\V p a\<Pa > +V pn - U pF . Potential < (pa\V p a\<Pa > 
is replaced by the optical potential U p a and we obtain a 
standard post form of the DWBA amplitude: 



M DW ( post \k pF , k, 



dA 



-< ^\AV pF \^ +) >, 



(7) 



where A V pF — U p a + V pn — U pF . Now we will transform 
this volume integral into the surface one. First, we adopt 
r n A and r pF as Jacobian variables and split the configura- 
tion space over r n A into the internal and external regions, 
while the integral over the second Jacobian variable, r pF , 
is taken over all the coordinate space. Splitting the reac- 
tion amplitude into internal and external amplitudes we 
get 



M DW ^ ost \k pF ,k dA ) 



■pF 



k dA ) 



■M. 



DW(post), 



<-pF 



(8) 



where the internal amplitude M?^^" 



is given by 



M 



DW(post) 
,(-) 



(k P F, Kia) 



-<X p F J lA 7 \&V pF \< PdXd + 2> 

Correspondingly, the external amplitude is given by 



(9) 



M^ ost \k pF ,\, dA ) 



<4- ) 7j|Ay pF |^xS > 



(+) 



r nA >Rn 



(10) 



Here, R n A is the channel radius similar to the one in- 
troduced in the i?-matrix approach, which separates the 
internal and external regions. 

The splitting of the amplitude into the internal a nd 
external parts in the subspace over the Jacobian variable 
r n A is natural and evident. The overlap function I^(r n A) 
is the only object in the reaction amplitude which pro- 
vides spectroscopic and structure information. In the ex- 
ternal region the overlap function has a standard radial 
shape given by the spherical Hankel function (for neu- 
trons) with the amplitude called the ANC (see below). 
To determine the behavior of the overlap function in the 
nuclear interior, which bring one of the main uncertain- 
ties in the analysis of the deuteron stripping, microscopic 
calculations are required 34] . In a standard approach the 
internal part of the overlap function is approximated by 
the single-particle bound state wave function calculated 
in the adopted mean field. The proportionality coeffi- 
cient is the square root of the spectroscopic factor. Due 
to the structure of the transition operator the external 
matrix element M®^^ ^ in the post form is dominant 
compared to a small contribution coming from the inter- 
nal part . This simple observation stems from 
the following. 

In the internal matrix element, r n A < RnA, due ab- 
sorption of the protons inside nucleus F, effective r pn ~ 
fpA ~ r pF > Rf, where R F is the radius of nucleus 
F. For the protons outside of F and neutrons inside or 
on the surface of A each nuclear interaction in the op- 
erator A V pF — U P A + Vpn — U pF is small. Potential 
U pF is arbitrary and often U pF is chosen to compensate 
for U p a so that the transition operator reduces to V pn - 
Since the DWBA is the first order perturbation theory, 
the minimization of the whole transition operator A V pF 
provides smaller higher order terms and, hence, better 
serves the theory. This choice is more preferable in the 
formalism presented here and we adopt U pF , which min- 
imizes A V pF = U p a + V P n — U pF at r n A < R n A making 
the contribution from the internal matrix element small 
compared to the external one. 

In the external matrix element (r n A > RnA), which is 
dominant, the overlap function 1^ can be replaced by 

its asymptotic tail. Although M®^^ ^ can be eas- 
ily calculated for stripping to the bound state, here we 
transform this matrix element into an alternative form, 
which has clear advantage in case of stripping to reso- 
nance states discussed below where convergence becomes 
a main impediment. 

Now we proceed to the transformation of the volume 
integral defining the external matrix element in terms 
of the dominant surface integral encircling the sphere at 
r n A = R n A and a small, due to the structure of the tran- 
sition operator in the prior form (see Eq. (]17jl). external 
volume integral in the prior form. Note that the trans- 
formation is exact within the DWBA formalism. 

To transform the external volume integral to the sur- 
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face one, we rewrite the transition operator as 



A V pF = U p A + Vp, 



u, 



pF 



-u dA ] 

-(U pA - 



- [U p f] 
UdA). (11) 



The bracketed operators are the right- hand- side opera- 
tors in the Schrodinger equations for the initial and final 
channel wave functions in the external region: 



(E - T) <p d Xd + 2 = (V P n + U dA ) <p d Xd + 2 (12) 



and 



(E-T)lZ XpF h = Up F lZxpF 



(-)- 



(13) 



To derive Eq. (|T3|) wc took into account that at r nA > 
RnA I a satisfies the asymptotic Schrodinger equation 
{snA — T u a) I a — 0, where is the binding energy of 
the bound state (ij) and Tij is the kinetic energy oper- 
ator of the relative motion of i and j. These equations 
imply the following connection between the external post 
form DWBA amplitude and the matrix element Mg W 
containing the surface integral: 



M. 



DW(post) 



(k P F, k, 



d A 



M° w (k pF: k 



d A 



-M. 



DW (prior) 



(V> k dA), (14) 



where 



M, 



DW (prior) 



(kpf, k d A) 



=<xU^\AV d A\'p d X d + 2 > 



r„A>RnA 



(15) 



the spherical harmonics and f = r/r is the unit vector. 
We assume that only one value of l n A contributes to ex- 
pansion ([T8|) . If the channel radius is taken larger than 
the range of the nuclear interaction, the radial overlap 
function can be replaced by its asymptotic term, 



*A j nA l nA ( H nA) ~ L 'Aj nA l ■« 



X K n A tl, (l KnA r nA ), (19) 



where h\ \ A (i K n AT n A) is the spherical Hankel function 
of the first order, C A ■ ; is the ANC of the overlap 
function, n n A — \J1 Jhui £nA is the bound state wave 
number. 

It is also useful to introduce the reduced-width ampli- 
tude used in the i?-matrix a ppr oach, which can be ex- 
pressed in terms of the ANC [25[ : 



JnA j nA l nA 



1 A j nA l„ A \ n nA) 



2fi 



ii A 



= \ ^^ A + 1 ^AC F A MlnA h^K nA RnA). (20) 

Correspondingly, the reduced width is 

RnA 



lnAj n A In 
RnA 



2/J>nA 



I 1 A j nA l„ 



,{RnA)Y 
a lnA hfJ A (iK n ARnA)] 2 ■ 

(21) 



and 



M 



g W (k P F, k d A) 



-< Y ( ~ ] I F 



-^\ VdXd + 2 > 



n A >Rn 



(16) 



Here, the transition operator in the prior form A V d A in 
the external region, where the nuclear n — A interaction 
disappears, takes the form 



A V d A = U p a — U d A- 
The overlap function is given by 



(17) 



I F (r nA ) = ^2 < J a M A jnA m jnA \J F M F > 

jnA m jnA m ^ A 

X < J„ M n l nA minA \JnA m jnA > 



x Y, 



l nA mi nA (*nA) J Aj nA U A ( r nA)- 



(18) 



Here, < ji mi jimi\jj,m^ > is the Clebsch-Gordan co- 
efficient, l n A {fni nA ) is the orbital angular momentum 
(its projection) of the relative motion of n and A, j n A 
(rrij nA ) is the total angular momentum (its projection) 
of n in the bound state F — (nA), Ji {Mi) is the spin 
(its projection) of nucleus i; I A i nA j nA { r nA) is the radial 
overlap function, which is a real function [23|, Yj m (f) is 



It is worth mentioning that, due to the presence of the 
channel radius R n A, the reduced width, in contrast to the 
ANC, is mo del- dependent. The dependence on the chan- 
nel radius becomes crucial with increase of the binding 
energy. We are going to use also the boundary condi- 
tion, which is the logarithmic derivative of the overlap 
function at r n A = RnA'- 



B n A = 



1 



d[r n Ahf\ (iK nA r n A)} 



hf' {iK nA RnA) 



dr 



Due to Eq. ([19]). the amplitude M { 



(22) 

DW (prior) 
ext Can De 

parametrized in terms of the ANC. We note that this am- 
plitude is also small. In the external region, r n A > RnA, 
the nuclear n — A interaction can be neglected. Besides 
in this region the overlap function exponentially fades 
away. Also, if the proton absorption is strong in the in- 
ternal region of A, the dominant contribution comes from 
TpA > Ra, where Ra is the radius of nucleus A. If the 
adopted radius channel R n A is larger than the n — A nu- 
clear interaction radius we can neglect n— A nuclear inter- 
action in the external region. In this region each nuclear 
potential U^ A and U^ A and their difference U p a — U d A are 
small. The Coulomb part U^ a ~U^ a k Z a e 2 R d /(2 R\), 
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where R d is the deuteron size and Za e is the charge 
of nucleus A, is also too small compared to the nuclear 
potential. Thus the dominant contribution to the post 
DWBA amplitude M®J (post) , Eq. JH]), and, hence, to 
the total post form DWBA amplitude M DW< - post ^ comes 
from the surface integral Mg W . Here and in what follows 
all the amplitudes with the transition operator 
are assigned the subscript S, which is abbreviation of 
" surface" , because the volume matrix elements of these 
amplitudes can be transformed into the surface ones in 
the subspace over variable r n A while over the second Ja- 
cobian variable r pF we always keep the volume integral. 

Now we express Mg in terms of the surface integral 
over variable r n A and the same technique will be used 
throughout the paper. The kinetic energy operator can 
be written as T = T pF + T n A ■ T pF is a Hermitian oper- 
ator in the subspace spanned by the bra and ket states 
in Eq. ([16]) . It can be proved if we take into account 
that at r p p —> oo the integrand in this equation van- 
ishes exponentially due to the presence of the bound state 
wave function ip d (r pn ) and the overlap function 7^(r n ^). 
Hence, integrating by parts twice the integral over r pF 
we obtain 



<X pF ] lA\T pF -^ pF \VdX { d + 2 > 
=<X { pF ] Ia\^ p f-T* p i 



r n A>R„ 



r„A>RnA 



(23) 



Then Mg w reduces to 



M 



--<X P F T A I TnA - ~^nA \ (fid X d + A > 



r nA >R„ 



(24) 



We apply now Green's theorem to transform the vol- 
ume integral into the surface one, which encircles the 
inner volume over the coordinate r: 



dr/(r) 



S(r) 



r<R. 



= -i- j dS [g(v) V r /(r) - /(r) V r g(r)\ 



r=R 



1 



= R 2 I dn T 

2 /i 



Or 



dr 



r=R 



(25) 



Here, d S = R 2 d il r f , where Sl r is the solid angle. Note 
that the unit vector f is the normal vector to the sphere 
directed outside of the restricted by the surface volume. 
The integration in Eq. (|2"4")l over r n A is taken over the 
external volume restricted by two spherical surfaces: the 
inner surface with the radius R n A and the external sur- 



face with the radius i?_ 



A I 



DW 



oo, that is 



(k pF , k dA ) = -M^ nA (k pF , k dA ) 
+M^(k pF , k dA ). 



<dA). (26) 

The first term in this equation is the surface integral 
encircling the inner surface of the external volume at 
r n A = R n A while the second term is the surface integral 
taken at r n A = R n A ~ * 00 ■ A negative sign in front of 
the first term appears because the normal to the surface 
is directed inward to the center of the volume, i.e. op- 
posite to the normal to the external surface (at infinitely 
large radius). The second term vanishes because of the 
presence of the overlap function JjT, which decreases ex- 
ponentially at r n A — > oo. Then for Mg W we get 



M° w (k pF , k, 



<IA 



(k pF , k 



d A 



2[i 



ti A 



dr pF X- + kL( r P f ) 



tPd(rpn)x kdA (*dA 



d[I F A {VnA)Y 



dtp d {r P n)x ( ^ A ){*dA) 



d 



TnA 



[Ia(^a)] 



'r n A 



- T nA =R nA 

(27) 



Here we took into account that xL ( r ) = X-k( r )- 
voking Eqs. (fl8|) and (p~9|) we can rewrite Mg W in the 
form explicitly showing parametrization in terms of the 



reduced width amplitude (ANC) and boundray condi- 
tion, the quantities used in the i?-matrix approach: 
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M° w (k pF , k dA ) 



1 K nA R n A hfJ A (i K nA RnA) ^ < Ja Ma i nA I j f M f > 



jnA ™.j nA m i nA M n 



x < J n M n l nA m lnA \j nA m jnA >< JpM p J n M n \J d M d > C A • K 



x J dr pFX % F {r pF ) J dn rnA Y* iAminA (r nA ) 



d¥d{Y P n)x^ A ){*dA) 



Vd{r P n) xi dA {*dA) (B nA - 1) - R n/ 



>r nA 



r n A = RnA 

(28) 



' RnA 
2^nA 



^2 < J A M A j nA m jnA \J F M F >< J n M n l nA mi nA \j nA mjnA ' > 



jnA m jnA mi nA M„ 



X < JpM p J n M n \J d M d > ~jnAj nA l nA J drpFX- 4 k pF ( r P^) / d ^r nA Y t * 



A^l nA 



JnA) 



<Pd(r P n) x'ClirdA) (B nA - 1) - Rnj 



d<p d {r pn )x^ A ){*dA) 



dr nJ 



r nA =R„ 



(29) 



by 



Finally, the total post form DWBA amplitude is given B. Prior form of DWBA. Stripping to bound state. 



M DW ^(k pF , k dA ) = M™ {post \KF, k rf A) 



+M, 



DW (prior) , 
ext 



(kp F , k dA ) + M° w (k pF , k dA ). (30) 



M. 



DW{post) 

DW( P nor) wg rewrite Eq. (|30[) in a different form: 



Taking into account that Mg w = M" x " t 



M DW( P ost) {kp ^ ^ = M DW( P ost) {kp ^ k ^ 



M^ p " or \k pF ,k dA ) 



+ 



M. 



DW( P ost) 
ext 



(k pF , k dA )-M ( 



DW (prior) 



(31) 



Thus, the main result of this section is that the post 
form of the DWBA amplitude can be written as the sum 
of the peripheral parts, y[®w{ P nor) j^dw ^ anc j sma n 

internal term M^^^ 1 " 3 ^. The peripheral part itself con- 
sists of the dominant surface amplitude Mg W and small 

external prior form M^^^ prlor \ The peripheral part is 
parametrized in terms of the ANC (reduced width ampli- 
tude) , channel radius R nA and the logarithmic boundary 
condition, that is in terms of the parameters used in the 
i?-matrix fitting. The model dependence of these two pe- 
ripheral amplitudes is caused by the ambiguity of the op- 
tical potentials and channel radius R nA . The strongest 
model dependence comes from M^^^ post \ because, in 
addition to the ambiguity of the optical potentials, to 
calculate it one needs to know the behavior of the overlap 
function in the internal region. For peripheral reactions 



contribution of M, 



DW( P ost) 



can be neglected. 



In subsection A the post form of the DWBA amplitude 
has been considered. However, all the results hold also 
for the prior form 

M DW{p " or) (k pF , k dA ) = < Xp -J /H AV dA \ VdX f2 > 



{ DW (pn or ){kpFj kdA) + M °J(Prior) {kpFj ^ 



(32) 



where 



M. 



DW (prior) 



(k p _p, k d A) 



=<x { i F ] I 1 A \&V dA \vdX d + A ) > 



(+) 



r nA <RnA 



(33) 



and 



M. 



DW (prior) , 



=< x y i F A \^v dA Wdx { Zl> 



r nA >R„ 



J,, -i ■ k,/ t 

<.p F T 

with the transition operator 

AF dj4 = u p A + v nA - U dA 



(34) 



(35) 



The n — A interaction potential V nA —< ip A \V nA \ip A > is 
the mean field real potential supporting the bound state 
(nA). The splitting of the amplitude into the internal 
and external terms in the subspace over the coordinate 
r nA helps us to further transform the prior DWBA ampli- 
tude. Due to the structure of the transition operator the 
external matrix element M^^^^ in the prior form is 
small (see the discussion in subsection lll A[) and the main 
contribution in the prior form comes from the internal 
part '( prlor ) _ Since the internal part is given by the 
volume integral, its calculation requires the knowledge of 
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the overlap function in the internal region. The model 
dependence of the overlap function in the nuclear inte- 
rior (r n A < RnA) brings one of the main problems and 
main uncertainty in the calculation of the internal matrix 
element. However, using the surface integral we can re- 
distribute the internal contribution in terms of dominant 
the surface term (over variable r n ^) plus small internal 
part written in terms of the volume integral in the post 
form. With reasonable choice of the channel radius R n A 
the contribution from the internal volume integral in the 
post form can be significantly decreased compared the 
surface matrix element. The latter can be expressed in 
terms of the i?-matrix parameters - the observable re- 
duced width amplitude (ANC), boundary condition and 
channel radius. To transform j^P^\P rlor > j n ^ ^he sur- 
face integral in the subspace over variable r n A we rewrite 
the transition operator in the internal region as 

A VdA = U pA + V nA ~ U dA 

= [V nA + U pF ] + (UpA + V pn - U pF ) - [Vp„ + UdA]. 

(36) 

The bracketed transition operators are the potential op- 
erators in the Schrodinger equations for the initial and 
final channel wave functions. Hence, for the internal prior 
form of the DWBA we obtain 

<r (p ™ r) (v,k dA ) 

= MZT {p ° st) (KFi Ku) + M^ pHor \kp F , k dA ), 

(37) 

where 

M° w (k pF , k dA ) = - < *y /ji ^ - f i % Xd + 2 > 

= -< XpF /£ I ^nA - T*nA\ Vd X d + 2 > 

= -M^Jk pF ,k dA ). (38) 

Note that here is the surface integral encircling 

the border of the internal volume at r n A = R nA with the 
normal directed outward. Thus we have demonstrated, 
what should be expected from the very beginning, that 
M DW( P nor) = M DW{ P ost)_ Hence all the equations ob- 
tained in the previous subsection III Al are also valid in 
the prior formalism. 

It is worth mentioning that in the post formalism, in 
contrast to the prior one, we have obtained two surface 
integrals (in the subspace over r n A) with the radii r n A — 
R n A and r n A = R n A ~^ 00 an d then proved that the 
second integral is zero. From the equality of the post 
and prior DWBA amplitudes we could conclude that the 
surface matrix element over infinitely large sphere r n A — 
R nA oo, which appears only in the post formalism, 
vanishes. 

There is another interesting point to discuss which ex- 
plains the advantage of the above outlined formulation 
of the stripping. As we have discussed, due to differ- 
ent structure of the transition operators in the post and 



prior forms, the main contribution to the post (prior) 
form comes from the external (internal) part (in the sub- 
space over variable r n A)- Since both forms give identical 
amplitudes, that is, describe the same reaction mecha- 
nism and the same physics, such redistribution of the 
main contribution is possible only if the main contribu- 
tion to each form comes from the border between external 
and internal parts. In the post (prior) form this border 
attributed to the external (internal) form and can be ex- 
pressed in term of the surface integral. Let us rewrite 
equality M DW ^ nor ^ = M DW ^ ost ^> in the following form: 

M°J^"° r \k pF , k dA ) + M^ rior \k pF , k dA ) 

= M°J^{k pF , kcu) + M°J^(k pF , k«u). 

(39) 

In this form the dominant terms are ^/[P^ r (.P rior > anc [ 
M^Y^ 05 ^ while the rest two terms, M^t vt " prlor ^ and 
M i ^^ pos *' ) are smaller. From Eq. (j39| we get 

M°Y ipost) (kp F , k dA ) - M™ (p "° r) (k pF , k dA ) 
= M™^ ior \kp F , k dA ) - M^ {P ° St) (KF, ^a) 
= Mg w (k pF , k dA ) = ~M^ a (k pF , k dA ). (40) 

Thus the difference between the post and prior external 
amplitudes (or the prior and post internal ones) is the 
surface integral in the subspace over r n A- 

There is one more point left to discuss. When deriving 
the post form of the DWBA amplitude from Eq. ((4]) we 
used approximation ip F « I a 'PA neglecting the contri- 
bution from the channels n + A n , n > 0, where A n is the 
excited state of A. However, I will show now that the 
surface integral formulation doesn't require this approxi- 
mation. To this end let us split M(P° st '> into the internal 
and external parts in the subspace over variable r n A- In 
the internal part we use a standard DWBA approxima- 
tion (p F w /J[ ip a to arrive to the standard internal post 
DWBA amplitude. In the external part we rewrite the 
transition operator as 

A Vp F = V p A + Vp n - UpF 

= -[V A + U pF ] + [Vpn + V A + U dA ] + (VpA - U dA ) . 

(41) 

The bracketed operators are the right-hand side opera- 
tors of the Schrodinger equations 

(E - T) $| +) = (Vpn + V A + U dA ) *| +) (42) 

and 

(E - T) <S> { f } * = (V A + U pF ) $^ } *. (43) 
Hence, the external part of M^ post ^ reduces to 
M^ st \k pF ,k dA ) 

= M s{ext) (k pF , k dA ) + M^W, k^), (44) 
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where 



r nA >Rr, 



and 



Ms(ext)O i pF, k 



dA 



(45) 



r n A>R n A 

(46) 



In the matrix element M, 



(prior) 



we can use a standard 



DWBA approximation cpp rs I^fA which leads to the 
standard external prior DWBA amplitude. The matrix 
element Mgr ext \ can be rewritten as 



= < X p ~F <Pf\ T nA ~ ~^ nA \fd H> A X d + A > 
=<X ( Jll\^nA-^nAW d X { d + 2 > 



r n A>R n A 

r nA >R n j 
r nA >RnA 



MJT (k pF , k dA ), 



We took into account that < $ 



-L A 



(+) 



/ 



pF + T / 



> = 0, where Ta is the internal motion kinetic 



energy operator of nucleus A, and T n A <-Pa = <PA T n A ■ 
Thus Mgiext) can be transformed to the surface inte- 
gral over variable r n A encircling the inner volume with 
the radius r n A = RnA without invoking approximation 
(pp w /JT if a- It means that, when deriving the post form 
of the DWBA amplitude, the approximation ipp rs I A tp A 

is required only to obtain two small terms, M®^^ post ^ 
and M^ x Y^ pnor \ but not the dominant surface term 
—Mg^ . In this sense the surface integral formalism 
is an improvement of the DWBA. 



Deuteron stripping to bound states. Post 
CDCC formalism 



In the previous sections we succeeded to parametrize 
the DWBA amplitude in terms of the ANC except for a 
small term, M®^^ post \ The most serious shortcoming 
of the DWBA is that it neglects the coupling to open 
reaction and breakup channels. This coupling can be 
taken into account if an exact wave function in the initial 
or final states is used. However, the exact wave functions 
are not yet available (if they would be available in the 
whole configuration space, we don't need to calculate the 
matrix element because the asymptotic terms of the exact 
wave functions provide the reaction amplitudes in all the 
open channels). Here we use the CDCC formalism, which 
takes into account the elastic d + A and the deuteron 
breakup channel p + n + A in the initial state. 



In this subsection the surface integral formulation of 
the reaction theory will be applied to the post form of the 
CDCC amplitude for deuteron stripping to bound states. 
It will allow us to parametrize the stripping amplitude in 
the CDCC approach in terms of the i?-matrix parameters 
- the reduced width amplitude, boundary condition and 
the channel radius. To obtain the CDCC wave function 
describing the initial state of the stripping reaction, first 
the exact initial scattering wave function ^ is replaced 

by the three-body wave function ty? B ( + ' , which takes into 
account the coupling of the initial channel d + A and the 
deuteron breakup channel p + n + A [J-Q and satisfies 
the Schrodinger equation (in the three-body p + n + A 
model space) 



{E — T — U p a — U n 



3B(H 







(48) 



with the outgoing waves in the elastic channel d + A and 
the breakup channel p + n + A. A general solution of 
this equation with the d + A incident wave has outgo- 
ing waves in the elastic, breakup and two rearrangement 
channels, n+(p A) and p+ (n A) . To damp rearrangement 
channels in the asymptotic behavior of the wave function 
the optical potentials U p a and U„a with strong 



3B(+) 



(47) imaginary terms can be used [35j. 



3B(H 



^ B{+ \r d A, r pn ) = Vd(j pn )x^ A {rdA) 



is given by 



J dpp„^^(rp„)xp(+)(p p 



)( r dA). 



(49) 



Here, fd(^pn) is the deuteron bound state wave function, 
x l ji ppL{ Y pn) the- p — n scattering wave function with the 
relative momentum p pn , X^] A {?dA) and Xp( Ppn ) { v <ia) 
are the expansion coefficients, EdA — £ pn — P 2 /(2 ^dA) + 

P 2 pn/( 2 fJ-pn)- 

In practical application the wave function ty^ 3 ^ is 
replaced by the CDCC wave function, which is a solution 
of the projected Schrodinger equation 

(E — T — Ufr ] - U { n P A pn) - V pn ) ttf DCC{+) = 0. (50) 

Here, U^ pn) = P pn U pA Ppn, and 



P„ 



x y; 



E E 

Z p „=0 m lpn =-l p „ 



d ^r p „ Y[ pn m lpn (fpn) 



l l P ,^n lpn \ L pn) ( 51 ) 

is the projection operator, which truncates the number of 
the spherical harmonics Y\ m , n {^pn) in the coordinate 
r pn . Application of this operator to the three-body wave 
function suppresses the rearrangement channels in the 
asymptotic wave function. The CDCC wave function is 
taken in the form 



CDCC(+) 



(rpn,T dA ) = Pp n £ ^{v' pn ) X ^ + \vdA), 



n=0 



(52) 
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where ippn( r P n) ~ Vdi^pn) is the deuteron bound state 
wave function, i[>pn (T pn ), n > 1, is the n-th discretized 
continuum state of the p — n pair obtained by averaging 
continuous breakup states in the n-th bin, Xi \ r dA) 
are the functions, which describe the relative motion of 
the center-of-mass of the p — n pair in the n-th state and 
A. Note that x\° ( r dA) asymptotically behaves as the 
incident Coulomb distorted d — A plane wave plus outgo- 
ing scattered wave, while x l ( " )(+) (i - dA) for n > asymp- 
totically do not contain any plane wave having only the 
outgoing scattered wave. 

To derive the post form of the CDCC amplitude from 
the exact one, first we replace the initial exact scattering 
wave function <fj +) by <£u#f B(+) . Note that $fB(+) ig 
the three-body model {p + n + A) wave function which 
treats nucleus A as a constituent particle leaving its inter- 
nal degrees of freedom intact. That is why the wave func- 
tion \f^ + ^ is approximated by the product of the bound 
state wave function if a and ^ B< > + \ Correspondingly, 
the transition operator AV p f — V v a + V pn — U pF is re- 
placed by A V p f = U p a + V pn — UpF ■ This replacement 
of the microscopic potential V p a in the exact post form 
amplitude by U p a is evident because the p— A interaction 
potential in the transition operator should be the same 
as the one in the Schrodinger equation for the initial scat- 
tering wave function ^\ B ^ + \ Potential V pn remains the 
same when we approximate the initial exact scattering 
wave function ny the three-body one. The final state 
optical potential U p p is arbitrary and we discuss the op- 
timal choice of this potential later on. These approxima- 
tions lead to the expression for the post form stripping 
amplitude in the three-body model in the initial state: 



M 3B ^ ost \k pF , k dA ) 

=<*y <p F \Av PF \<p A *i B (+) > 



< X pF 1 A 



I F A \AV pF \ i S> 



> 



the post form of the CDCC amplitude takes the form: 

M CDCC( P ost ){KFj ^ 

=<X ( p- F ] ll\AV P pT\^ DCC{+) >- (54) 

Now we split M CDCC into the internal and external parts 
in the subspace r n A- 



M CDCC ^ ost \kp F ,k dA ) = M i 



_ A *CDCC(post) 



+ M, 



CDCC(post) 



(kpF) k^). 



(55) 



The internal amplitude M^ GG ^ post ^ is given by 



M 



CDCC(post) 



in I 



(k P F, k 



dA 



=<x p f^|a<f^ cc(+) > 



(56) 



Correspondingly, the external amplitude is 



M. 



CDCC{post) 



(kpF, Kia) 



ext 

=<xtf/j[|AT^|*f DOO W> 



r„A>R n 



(57) 



I remind that the integral over the second Jacobian vari- 
able, r p F, is taken over all the coordinate space. Sim- 
ilarly to the DWBA case, the internal part is small if 
the channel radius R„a is not too large. Due to the 
strong absorption of the proton inside A, which is con- 
trolled by the imaginary part of the optical potential 
U the effective distances are t v a > Ra- Besides, 
in the internal region, r n A < RnA, and large r p A, where 

tP-o: 



TpA ~ r V n = \vp A - r nA \, U p l n + V P n can be well ap- 
proximated by a properly chosen optical potential U pF 
p 

minimizing A V pF and the internal matrix element. The 
next step is to transform the external matrix element to 
the surface one. To this end we rewrite the transition 
(53) operator in the form 



Thus, even if we treat the d + A collision in the initial 
channel in the three-body approach, the final state con- 
tains the overlap function, which is essentually many- 
body object. Eq. (|53[) is impractical to use because it 
requires the knowledge of the three-body wave function 
Eq. (|49p . which contains unknown expansion co- 
efficients Xk dA (rdA) and XP(p pn )( T dA)- In practical appli- 
cations the is approximated by the CDCC wave 

function ^ CDCCi * + \ which requires the knowledge of the 

finite number of the expansion coefficients. They can be 

found from the coupled equations. Correspondingly, the 

transition operator AV pF = U p a + V pn — U pF in Eq. (|53")l 
p p 

is replaced by AV pF = U + V pn — U pF . Note that 
only the potential U p a{y p a), where r p A = r d A + 1/2 r pn 
is affected by the projector P pn - Then the expression for 



u pF = [-U pF ] + [u^ n + v, 



v pn\ 



(58) 



The bracketed operators in (f58|) are the right-hand-side 
potential operators in the Schrodinger equations in the 
external region r n A > RnA, where the nuclear n — A 
interaction vanishes: 

(E - T) tff DCC{+) = {U p l n + V^) tff DCC{+) (59) 

and 

{E-T)j$*I%* = U pF £->lZ*. (60) 
Note that the second equation follows from 

(s nA - T nA ) I A =< <PA\V n A\ Vf > ■ (61) 
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In the external region, r n A > R n A, the source term 
on the right-hand-side disappears and Eq. becomes 
evident. Taking into account Eqs (|59l> and (f60|) we get 



M, 



CDCC(post) 



(k pF , k dA ) = M s 
<X pF 'lZ\T-T\V 



_ ,rCDCC(post) 



(kpF, k d A) 



> 



r„A>R„ 



(62) 



where T = T pF + T n A- Here, as in the previous section, 
for the surface integral we use the subscript " S" . Since 
the CDCC wave function doesn't propagate into the final 
state (its asymptotic terms have only elastic and breakup 
terms) the operator T pF is Hermitian, i.e. 



<x { Ji. 



pF 



pF 



pF 
2* 



pF\ 
T CDCC(+) 



r nA >R. 



r nA >R nA 
= 0. 

(63) 



It can be also shown explicitly taking into account that 
the volume integral over r p p can be transformed into 
the surface integral over the sphere with the radius 
r p p = R pF — > oo. Since the overlap function decays 
exponentially at r n A — > oo, the integration over r n A is 
limited. Hence, at r pF — > oo using Eqs (|B6[) we get 



l/r pF , see Eq. (|B14|) . Hence the surface integral vanishes 
at R nF -> oo as R^p/Rpp -> 0. 



IpF 

Then M 



CDCC(post) 



takes the form 



CDCC(post) 



M 



= <X P F ] lA^nA-^nA^ 



CDCC( + ) 



> 



r nA >R„ 



-M. 



CDCCipost) 



(k pF , k dA )+M^ CC{post) (k pF , k dA ). 

(64) 



Thus, the volume integral at r n A > R n A in the ma- 



trix element M, 



CDCC( P ost) 



can be written as the sum of 



two surface integrals encircling the external volume, the 
sphere with the radius r n A = R n A and the sphere with 
r n A — R n A ~~ ^ 00 • Note that the integral over r p p is 
taken over all the coordinate space. Evidently that the 
integral over the infinitely large sphere vanishes because 
the overlap function I~[ exponentially decreases. Hence, 



M, 



CDCC{post) 



(k pF , k dA ) = -M 



CDCC{ P ost) 



{k p Pi k, 



dA 



(65) 



that r d A ~ r p p — > oo and r pn ~ r pF — > oo. The first 
term of the CDCC wave function decays exponentially 
at t p f — > oo because of the presence of the deuteron 
bound state wave function. The terms with n > 1 decay 
as l/Vpp [36|. The distorted wave X p ~p*(Ppf) decays as j^CDCCipost) 



The negative sign in front of the inner surface integral 
appears because the normal vector to the inner surface 
is directed to the center, i.e. opposite to the direction of 
the normal to the external surface at r n A = R n A 00 • 
Now we can use equations from subsection A replacing 
the initial channel wave function by the CDCC one. For 



we get 



r?2 

n nA 
2 UnA 



R 



2/J.nA 



M 



CDCC{post) 



(kpF, k dA ) = -Mg®° c ( post \k pF , k dA ) 

[lI(rnA)}*&r nA -^r nA )*? DCC{+ \rp F , r nA ) 

CC + \r pF , r n A)~ 



dr pF£ + k ' M / dS!,„„ [*f DCC<+) (r,f> 



dr n 



r nA >R„ 



r nA =R nA 

(66) 



Natural Jacobian variables for ^ DCC ^ +S > are VdA anc j 
r„„, but here we use another set of Jacobian variables, 



r pF and r n A- Taking into account Eq. (JT5J) and (|19p wc 

get 
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M 



CDCC(post) 
S 



(k pF , k dA ) = -M Sr 



CDCC(post) 



k p F, k dA ) 



Rn 



2 fJ>nJ 



E 



< J a M A j nA m jnA \J F M F > 



3nA "lj„A mi „A 



x < J n M„ Z„a m/ nA | j„a mj nj4 > 7„a j„ a f nJl / d r pF X^ pf , (r pF ) / d fi r „ A Y", 



CDCC(+) 



{v pF , Y nA ) (B nA - 1) - RnA 



d^ DCC ^(r pF ,r nA ) 
dr nA 



rnA=Rn/ 



(67) 



Note that the CDCC wave function itself also depends on 
quantum numbers of p — n and d — A subsystems, which 
we don't specify here. It will be done in the following up 
paper where concrete calculations will be presented. 

Thus we have obtained a remarkable result: the post 
form of the CDCC amplitude, in contrast to the DWBA 
one, is given by the sum of only two terms: 



M CDCC( P os t){kp ^ ^ = M CDCC {p os t){K ^ ^ 

0. 

(68) 



. M CDCC (post){kpF)kdA) 



where the first term, which is the internal post form of the 
CDCC amplitude, can be minimized by a proper choice 
of U pF and the channel radius R nA , while the second 
term, which is dominant, represents the surface integral 
with the radius R nA , which encircles the internal volume 
in the subspace over the coordinate r nA . If the channel 
radius is larger than the n — A nuclear interaction radius 
the second term is parametrized in terms of the reduced 
width amplitude (ANC of the projection of the bound 
state wave function of F on the two-body state n+A) and 



the boundary condition at r nA = R nA . 
is small enough, 



If M 



CDCC(post) 
ini 



M CDCC ^ ost \k pF , k dA ) « -AT 



CDCC(post) 



(69) 



Thus we succeeded to parametrize the post form of the 
CDCC amplitude in terms of the i?-matrix parameters. 
Eq. (|68[) and parametrization of the surface term of the 
post CDCC amplitude in terms of the i?-matrix param- 
eters, Eq. (1671) . are one of the main results of this paper. 



r CDCC(post) 



can be min- 



Although it is assumed that M int 

imized so that the second term in Eq. becomes 

dominant, I would like to present a different form for 

Mint >CCt " P ° St \] i p F , kd A ), which leads to a different form 

for the whole amplitude M CDCC< -P ost ^{k pF , k dA ). To 

p 

this end, let us rewrite the transition operator A V pF " 



in M 



CDCC(post) 



k pF , k dA ) as 

UpF 



PpA + U nA + Vpn] ~ [VnA + U pF ] + V nA 



u 



(70) 



Here, V 71 a is the mean field potential supporting the 
bound state (nA) while U n v A is the projected optical 
potential describing the n — A interaction in the initial 
state of the reaction and entering the Schrodinger equa- 
tion for the projected CDCC wave function in the initial 
state. The bracketed potential operators are the right- 
hand-side operators of the Schrodinger equations in the 
internal region, r nA < R nA , 



and 



CDCC{+) _ (TT P pn 
— \ U pA 



(71) 



(E T) XpF } * I F A * = (V nA + U p p) Xp t I a * ■ (72) 



Replacing the bracketed potential operators [U pA n + 
UnT + V pn ] and [V nA + Upp] by E - f and E - T, 



correspondingly, we get for M 



CDCC(post) 



a new form: 



M anaa {post)(kpFjkdA) 



M 



CDCC{post) 



(k pF , k dA ) 



+ Mg° cc ^(k pF ,k dA ), 



M f„^«(k pF ,k dA ) 



(73) 



V pF 



Ppn^CDCC(+) 



> 



r nA <R n A 



AV^ A =V nA -U^ A 



(74) 



(75) 



Then the total post form of the CDCC amplitude can be 
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written as 

M CDCC( P ost) {KF) ^ = M CDCC {p ost) {K ^ ^ 



dA) 



rt nA 



CDCC(post) 



(k pF , k d A) 

M CDCC( P ost Hkp ^ ^ = M C S CC( pOS t) (kpF) kdA) 



r nA <R n A 



(76) 



Thus, we obtained another important result. The CDCC 
amplitude in the post form is equal to the inner volume 
integral over variable r n A with the transition operator 



u 



nA ■ 



This transition operator is the difference 



between the bound state potential V n A supporting the 
final bound state (n A) and the projected optical poten- 
tial describing the n — A interaction in the initial state. 
It is worth mentioning that Eqs (|68p and (f?6"|) are ex- 



act within the CDCC approach. If j^^ 00 ^ ^ [ s sma vj 
enough, then 



M CDCC( P o st){k 



■pF, K dA) 



M 



CDCC(post) 



(77) 



However, I prefer Eq. 



rather than (1751 . To calcu- 



late Mam CC ^ P ° st ^ one needs to know the overlap func- 
tion in the internal region, where the overlap function is 
model-dependent and requires microscopic calculations. 
In contrast, in Eq. (|68l) the dominant part is the surface 
integral, which is parametrized in terms of the reduced 
width amplitude (ANC). The model dependence of the 
surface part is related with the ambiguity of the optical 
potentials and the value of the cut-off orbital angular mo- 
mentum in the p — n subsystem in the CDCC approach. 
Comparison with experiment allows one to extract the re- 
duced width amplitude. The model-dependent internal 
part in Eq. (155)) is small. Eqs (f6"5f and (I75|) is prelude 
to the theory of the stripping to resonance, where the 
convergence problem of the external part is one of the 
main issue. As we have demonstrated in the post CDCC 
formalism the external part doesn't appear at all. It re- 
solves the convergence problem related with the external 
part. 



D. Deuteron stripping to bound states. Prior 
CDCC formalism 



A priori, the amplitudes in the post and prior forms of 
the CDCC formalism are not equal. That is why the ob- 
tained equations using the surface integrals are expected 
to be different in both formalisms. The prior form of the 
CDCC stripping amplitude is 



-< * 



CDCC{- 



(78) 



where 



A v dA - U pA 



+ V n A — UdA- 



(79) 



The projected CDCC wave function in the final state is 
a solution of the three-body Schrodinger equation 



{E-T-U^ 



V, 



nA 



Vp^)^ 000 ^* =0. (80) 



Here, 



PnA = 



E E 

1„a=0 mi =—l nJ 



r„A J lnA mi 



(fnA) 



X Y, 



InA mi 



(Ka)- 



(81) 



is the projection operator, which truncates the number of 
the spherical harmonics Yi nAmi (r n ^J in the coordinate 
r n A- 

Now, as usually, we split the amplitude M CDCC( - prior "> 
into the internal and external parts in the subspace over 
variable r n A- 



(V, k dA ). M ODaa <* rt °*Q hF , k dA ) 



= M 



CDCC(prior) 
int 



(k pF} k 



dA) 



M. 



CDCC(prior) 



(k pF , k d A), 

(82) 



where 

M CDCC ^° r \k pFl k dA ) 



=< f 



CDCC(~ 
f 



\UX+VnA-U d AWdX d + 2 > 



r n A<R n A 

(83) 



and 



M, 



CDCC(prior) 



{k p F, k d A) 



-< * 



CDCC(~) 
f 



\U^a A -UdAWdX^ > 



( + ) 



r„A>RnA 



(84) 



The external part of the prior amplitude (see discussion 
in subsection III Bp , due to the structure of the transition 
operator, is small and the dominant contribution comes 
from the internal amplitude. We will rewrite this ampli- 
tude singling out the surface integral over variable r n A- 
To do it we rewrite the transition operator 



pA 



Av p d r = u p p r 

+ v nA + V%*]- 



+ VnA — UdA 
[V pn +UdA 



]-+(Vpn-V p ^) 

(85) 



The bracketed operators are the right-hand-side opera- 
tors of the Schrodinger equations 



[E - T) * 
and 



CDCC(-)* 
/ 



pA 



V P nA ^iuCDCC(-)* 
' pn J 



(86) 



(E - T) ifd x ( d + 2 = [Vpn + U dA ) Vd XdA ■ (87) 
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Taking into account these equations we can rewrite 

CDCC(prior) , 
int 



MT„Y^^ yl " lu ' '(k pF , k dA ) in the form 



M, 



CDCC{prior) 



(k pF , k dA ) = M 

+Ae cc(pnor) (v,^), 

where 



CDCC{prior) 



(k pF , k, 



dA 



M CDCC( P rior) {KF) ^ 



-< 



CDCC(-) 



V. 



pn 



VdXdA > 



r nA <R nA 



(89) 



and 



CDCC(pr 



= - < * 



(k pF , k. d A 

CDCC( 

f 



T -T 



(+) 

Vd XdA 



> 



\<R n 



(90) 



Here, the kinetic energy operator T — T pF + T n A- In 
j^jCDCC(pnor) ^ e vo i ume integral over r pF can be trans- 
formed into the surface one taken over the sphere with the 
infinitely large radius: r pF = R p p — > oo. For r nA < R n A, 
due to the presence of the deuteron bound state wave 
function, the integrand goes to zero exponentially, that 
is this surface integral vanishes. Hence, only the sur- 
face integral encircling the inner volume with the radius 

TnA = RnA- 



M 



CDCC(prior) 



(k pF , k t 



dA 



< f 



CDCC(- 
f 



T nA - T nA | (p d x ( d + 2 > 



r nA <Rr, 



M CDCC { pos t){KF ^ dA) _ 



M, 



is 



given by Eq. 



(91) 

(EH). 

(k p F, k dA ) is an auxiliary internal part, 
which is small because at r n A < R n A and r p p > Rp due 
to the proton absorption in the nuclear interior, p — n 
nuclear interaction is significantly depleted, and so the 



CDCC(post) 

Sl *nA 

CDCC{prior) 



difference V, 



y pn 



Then 



M CDCC ^" or \k pF , k dA ) 

= MZ D x CC ^° r Hk pF , k dA ) - Ml DCC ^° st \k p p, k dA ) 

+ M^ CC{prlor \k pFl k d A), (92) 

Thus the total prior form CDCC amplitude consists of 
three terms, small auxiliary internal part, small external 
prior form and the dominant surface term. We can see 
that post and prior CDCC formalisms are not equivalent. 
In the approach used in the paper the configuration space 
over variable r n A was split into the internal and external 
parts. As it has been discussed in Introduction, such a 
splitting is natural because the main object of interest in 
the analysis of deuteron stripping is the overlap function 
I A of the bound states wave functions of the target A 
and final nucleus F. Its external part (r n A > Rua) is 



parametrized in terms of the observable ANC while the 
internal part is mo del- dependent. 

In the post formalism the external part is domi- 
nant. Invoking the post CDCC formalism allows us to 
rewrite the external CDCC matrix element in the form 
of the surface integral over variable r n A, which can be 
parametrized in terms of the parameters used in the R- 
matrix method for binary reactions, while the model- 
dependent internal part gives small contribution. Thus 
the volume part of the matrix element over variable r n A 
is transformed to the surface integral. For transfer to 
bound states such a transformation doesn't bring any sig- 
nificant advantages because the volume matrix element 
converges. However, for stripping to resonance states 
(see subsection MI C|) this transformation provides a de- 
cisive benefit because it solves the convergence problem 
of the matrix element. Here, the transformation of the 
post CDCC matrix element has been presented mostly 
for demonstration but the results will be used below in 
subsection MI CI for stripping to resonance states. 

The prior CDCC formalism would be preferable if we 
split the matrix element into the internal and external 
parts over variable r pn to separate the internal and pe- 
ripheral parts of the deuteron bound state wave function. 
But this wave function is well known and is not an ob- 
ject of study. That is why below, when considering the 
stripping to resonance states, we use only the post CDCC 
formalism. 



III. DEUTERON STRIPPING INTO 
RESONANCE STATES 

Now we proceed to the main goal of this paper, for- 
mulation of the deuteron stripping into resonance states 
using the surface integrals what will lead us to the gener- 
alized i?-matrix approach for the stripping into resonance 
states. Let us consider the deuteron stripping 



d + p + b + B. 



(93) 



We assume that the resonance formed in the system 
F = A + n can decay into channel B + b, which can 
be different from the entry channel A + n. We start from 
the post form and transform it to the surface integral 
following the method applied for the stripping to bound 
states. Now the application of the i?-matrix approach 
looks natural. Although we consider the deuteron strip- 
ping leading to a specific final channel d + A^p + b + B, 
there can be a few open channels coupled to the chan- 
nel n + A, which is formed after neutron is transferred 
to the target A. As in the previous sections, follow the 
i?-matrix approach, we split the integration region over 
r n A into two regions: internal and external. Internal re- 
gion is determined as the one where all open channels are 
coupled with each other, so that the transition from one 
channel to another can occur only in the internal region. 
The external region is the one where all the channels are 
decoupled. We obtain new forms for the DWBA and 
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then for the post form of the CDCC amplitude. For the 
DWBA both post and prior approach will lead to the 
same final expression. In the standard approach the post 
form of the DWBA amplitude is mainly contributed by 
the external part in the subspace r n A, where the con- 
vergence question of the DWBA matrix element, which 
contains the integration over r p p and r nA , becomes a 
main issue. In the prior form the main contribution to 
the DWBA matrix element mainly comes from the inter- 
nal region in the subspace r n ^, where a strong coupling 
between different open channels becomes an issue. In a 
new approach formulated below the DWBA amplitude 
(in the post and prior forms) is written as the sum of 
three amplitudes: small internal post and external prior 
forms, and the dominant surface integral in the subspace 
over r n A- This surface term is parametrized in terms 
of the reduced width amplitudes, resonance energies and 
boundary condition, that is the quantities used in a stan- 
dard i?-matrix approach. In the post CDCC approach 
the amplitude is given by the sum of the small inter- 
nal post form and the dominant surface term, that is, in 
contrast to the DWBA, no external prior form appears 
in the CDCC method. This resolves the issue of the con- 
vergence for stripping into resonant states. 



A. Stripping to resonance states. Post form of 
DWBA. 



The post form of the DWBA amplitude can be ob- 
tained by generalizing the corresponding equation for the 
deuteron stripping to the bound state. As a starting 
point, we use Eq. (J7]) in which, to get the amplitude 
for the deuteron stripping to resonance states, we should 
replace the overlap function 1% by the exact scattering 

wave function with the incident wave in the channel 



b + B: 



M DW(post)^ k 



dA) 



,(+) 



where AV p f = U p a + V pn — U p f and 



OtS k bB —K-bB 



(94) 



(95) 



Since we consider the stripping to the resonance state, 
which decays into two fragments b and B, there are three 
particles, p, b and B, in the final state. Hence, the kine- 
matics of the final state of the reaction depends on two 
Jacobian momenta, for which we adopt the relative mo- 
mentum of two fragments b and B and by the momen- 
tum corresponding to the relative motion of the exit- 
ing proton and the center of mass of the system b + B. 
Thus the deuteron stripping reaction amplitude depends 
on the momentum P — {k p F, k^s}, which is the six- 
dimensional momentum conjugated to the Jacobian co- 
ordinates of the system p + b + B Y = {r p p, r^}. 



Then repeating the steps used in derivation of the ex- 
pression for the post form of the DWBA amplitude for 
deuteron stripping to the bound state we get 



M DW( P ost)^ ^ 

_l_ M DW (P rior ) ( p 1, 
+ M ext l-n k dA 



M^ w (P,k dA ). 



(96) 



Here, internal post amplitude M®Y^ post \p, k^) and 
external prior amplitude ^ w%or ^ {P , k^) are given 



by 

M™ [p ° St) {P, *dA 



^ X p F 



\- 1 J 

Tlf^\AV pF \^l> <n 07) 



and 



DW (prior) 



(P, kdA) 



z < X pF 



(-) r^(ext)(-) ] 



AV pF \<PdXdA > 



r„A>RnA 



(98) 



Here, T 



(*«*)(-) 



•.(ext) (- 



(r„A) 



\r nA ) =< ip A \y 



=< 

(exi)(-) 
bB 



(»"*)(-) 
bB 



> and 



>. 



The last term of Eq. (|96|) . which will be transformed 
to the surface integral, is 



M DW 



dA 



ApF L nA 



-?\<Pdx'&> 



l>Rn 



(99) 



Let us discuss the advantage of this new form of the 
DWBA amplitude for the deuteron stripping to reso- 
nance state(s). Since the internal part is 
given by the volume integral, its calculation requires 

the knowledge of ^S> bB in the internal region. The 
model dependence of this function in the nuclear inte- 
rior (r n A < Rua), where different coupled channels do 
contribute, brings one of the main problems and main 
uncertainty in the calculation of the internal matrix el- 
ement. However, as it has been discussed in subsection 
III Al this matrix element gives a small contribution to 
the total post form amplitude M c ^ p ° st ! due to the 
structure of the transition operator A V p f and constrain 
r nA < RnA- These arguments are also valid when con- 
sidering the stripping into resonance states. A proper 
choice of the optical potential U p f and the channel ra- 
dius R n A may significantly reduce the contribution from 
the internal post form DWBA amplitude. Due to the 
structure of the transition operator AV dA, which has 
been also discussed in subsection III Al the external ma- 
trix element M® x Y^ p " or ^ in the prior form is also small 
and in some cases, with reasonable choice of the channel 
radius R n A, even can be neglected. Note that, in order to 
keep M^ V< " post ' > small, the channel radius R n A cannot be 
too large and, in order to keep ^^^(p" 01 ") small, cannot 
be too small. Thus with an optimal choice of the chan- 
nel radius the dominant part is the surface part Mg W , 
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which contains only one volume integral over r p p. Eq. 

which presents a new form of the DWBA amplitude 
for stripping to resonance states, is quite important for 
analysis of the stripping to resonance. 

Using the i?-matrix representation of the scattering 
wave function we are able to express the to- 

tal DWBA amplitude in terms of the reduced width 
amplitudes, level matrix, boundary condition and the 
channel radius, that is parameters used in a standard 
i?-matrix method to analyze binary resonant reactions 
n + A — > b + B. Since the reaction under consideration 
is the deuteron stripping, the presence of the deuteron 
in the initial state and exiting proton causes the distor- 
tions. That is why the reaction amplitude, in addition to 



the i?-matrix parameters describing the binary subpro- 
cess, contains additional factors - distorted waves in the 
initial and the final states. That is why we can call the 
obtained expression for the DWBA amplitude a general- 
ized i?-matrix for deuteron stripping to resonance states. 

Now we proceed to the derivation of the expressions for 
each amplitude in the right-hand-side of Eq. (IM1) and the 
total post form DWBA amplitude. Since the stripping 
into resonance states can lead to rearrangement, the exit 
channel b + B may differ from the entry channel n + A. 

To proceed further we now use the equations for 
obtained in Appendix [A] Taking into account Eqs. (|{J5)) 
and (|A1|) we get 



M°7 ipost \P,k dA ) 



2?r k, 



N 



^2 <sm s l mi\J F Mp >< J n M n J p M p \J d M d > 

Jf M f , I mi M n 



e- i ^-i i F i * TOi (-k 6B ) F»»B s ljAE b B)} 1/2 [A-% T <X p - f ) KIa P \&Vpf\^X { £ > 



r nA <R n A 



(100) 



In this equation we assume that the channel spin s and 
its projection m s in the exit channel c = b + B are fixed 
41]. Jf is the resonance spin (Mjr its projection) in the 
subsystem F = n+A = b+B and I is the b+B orbital an- 
gular momentum in the resonance state. The sum over 
Jf and I assumes that a few resonances with different 
spins may contribute to the reaction. The subscript c 
used in Appendix [A] for the channel b + B is replaced 
here by bB. Also ^ J T F nA F =< ip A \X J " MF > is projec- 
tion of X^ F Mf introduced in Appendix [A] on the bound 
state ip a- The bound-state like wave function X JfMf 



describes the system F = n + A = b + B in the internal 
region. A priori, it can be calculated using, for example, 
the shell model approach [13] ■ In Appendix [A] X? F Mp is 
written as a nonorthogonal sum of coupled channels, see 
Eq. (IA4[) . If we neglect the contribution from the chan- 
nel c, then ^ J T F nA. F can be approximated by the internal 
part of the overlap function, see Eq. <p~8|) . Taking into 
account this equation (rewritten in LS'-coupling scheme) 
we get 



M 



DW{post) 



(P, k, 



dA 



2 7T / k bB 
kbB V fJ-bB 



E 



i < s m s lmi\JF Mp > < s' m s > I' mv \ Jf Mp > 



Jp- Mp s' I V m s f m-i M n 



< J n M n J A M A \s'm s , > < J n M„ J p M p \J d M d > e - t&bBl Y z * (-k 6B ) 



N 



E lr»bBsij F (E b B)} 1/2 [A-V <x pF ) Y l t mi ,(r nA )I A 7 s , l , jF (r nA )\AV 

z/,r — 1 

i 



pFlVdXdA > 



r nA <Rr, 



(101) 



Here we added the sum over the channel spin s (its pro- 
jection m s i) in the entry channel c! = n + A of the reso- 
nant subreaction n + A—>F—>b + B and over the n + A 
orbital angular momentum The sum over M n and 
s' (m s i) appears because the transferred neutron is in- 
termediate (virtual). It is important that with a proper 



choice of the optical potential U p f the matrix element 



M 



DW(post) 



can be minimized so that its model depen- 



dence wouldn't have impact on the total matrix element 

M DW(post)^ 



To obtain the expression for M t 



jDW(prior) 



we use for 
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the external part i which can be obtained from 

Eq. (|A3ip , assuming that the resonance contribution to 
this wave function is dominant. In the sum over Jp in 
Eq. (|A31[) we keep only those total angular momenta at 
which resonances contributing to the reaction occur. Let 
us consider two possible cases. 



(i) The exit channel c = b+B in the resonant sub-process 
n + A — » b+B is different from channel d = n + A. In this 
case the external resonant wave function is given by Eq. 
(|A36|) and its projection on the bound state £ C ' = f>A is 
determined by Eq. [[A37]) . Then M™ [pnor) reduces to 



M 



DW (prior) 



hi, D \/ 7>„ A t— < 



kbB V v nA 



i l < Imi s m s \ Jp Mp > < I' my s' m s i \ Jp Mp > 



J F s' 1 1' m s i mj my M n 

x < J n M n J A M A \s'm s , >< J n M n J p M p \J d M d > Y? mi (-k 6B ) S^ s , 

X < XpF ~^— r Y l* mi ,( r nA)\ AV d A\^PdX d A > 



r nA >R„ 
I 



,A s' I' 



(102) 



Here, V dA is given by Eq. (|35l) . In the external region 
V n A = and V 'dA = U p a — UdA- Also has been added 
the sum over the orbital angular momentum I and its 
projection mi (I' and me) in the exit (entry) channel 
c = b + B (d = n + A) of the resonant sub-reaction n + 
A — > b + B, the sum over the channel spin s' and its 
projection m s < in the entry channel d = n + A of the 
resonance subprocess n + A — > b + B and the sum over 



M n because the neutron is the transferred particle. The 
projections of the spins of the incident deuteron Md, the 
exiting proton M pi the channel spin s and its projection 
m s of the exiting particles b and B are fixed. We also use 
the symmetry of the S matrix: S^F., — S Jp 



The matrix element S*/r 



bB s l: nA s' I' 



l';csl c s I; c' s' I' ' 

is given by Eq. (|A45|) . 



Substituting it into Eq. (|102[) gives 



M. 



DW (prior) , 



ex I 



kbB V v nA 



i < I mi s m s I Jf Mf > <l' mv s f m s f I Jf Mf > 



Jf Mf s' I V m s i mi m ; / M n 

x < J n M n J A M A \s'm s > >< J n M n J p M p \J d M d > Y t * mi (-k bB ) 



x e 



iStir „-iS h " „ \T (V Ml/2rA-ll \T- IV Ml/2 Ol>(k n A, R 
bBi e nAl > ^ [1 v bBsl J f (Hb)\ ' [A \ UT [1 mAs'V J F \bnA)\ ' p 



(-) OjtjknA, r nA ) 
X < XpF 



Rr, 



Y*(v nA )\AV dA \fdX d + 2 > 



r nA 0*,(k n A, RnA] 



r nA >R nA 



(103) 



Now we take into account that 



Offe Rc) = xjFf(k £ , R e ) + G?(k £ , R 5 ) 



% 0_ r 



(104) 



which for the channel c = c' = n + A and / = Z' takes the 
form 

Ol>(k nA , RnA) = \J Fp\%U, RnA) +Gf,(k nA , RnA) 
z &rct3,ri — 7-5 ^ r- 

X e °i' ( k n A . R 7i A > 



where in the absence of the Coulomb interac- 
tion Fi{p) = (^P/2) 1/2 Ji+iM and Gfo) = 
(-1) 1 (n p/2) 1 /' 2 J_ {l+1/2) (p), J±(i+i/2)(p) are Bessel 
functions. 



= ^F^(k nAl RnA) + Gf,{k nA , R n A)e lS ^'', (105) Then using Eqs. (fAHl and (fTU5j) we get 
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M. 



DW (prior) 



(P, k dA ) = 2 7T 



MfcB fcbB PnA 



E 



£ <lmi sm s \JF Mf > 



Jf Mf s' I V m s > mi m ( / M n 

x < l'm v s'm s ,\J F M F > < J n M n J A M A \s'm s , >< J„M n J p M p \J d M d > ^* mi (-k 6B ) 



N 



e -is b g, [T vbBs ij F {E hB )]V 2 [A-% TlTnAs , V j 

RnA 



U,T — 1 



(_) OMknA, r nA ) 
X < X P F — 



Of, (knA, R nA ) 



^ mi ,(w)|Ay dA |^x^ > 



(106) 



(ii) If c = c', that is 6 = n and B = A. Here two cases 
are possible: non-diagonal transition for which s =/= s' 
or/and I =/= I' and diagonal transition with / = I' and 
s = s' . The amplitude for the nondiagonal transition can 



be obtained from (|102|) . Here we present the expression 
for the diagonal transition (elastic scattering) amplitude, 
which can be obtained taking into account Eq. (|A33p : 



DW (prior) 



ext 



(p, k dA ) = i 



2tt 



k n A RnA 



E 



i < I mi s m s | Jf > < I m i' s m 8 > \ Jf Mf > 



Jf Mf I m s / mi mp M n 

x < J n M n J p M p \J d M d >< J n M n J A M A \sm s , > Y,* mi {-k nA 

RnA TA » m A T7 i (+) 



i _ oJf 

nAslinAsl 



Oi(k nAl R nA ) 



(-) OT(k nA , r nA ) 



T nA Oj~(knA, RnA) 



Y x * (r nA )\AV dA \<p dX Y A > > 



r nA >R„ A 



(107) 



Substituting the expression for the elastic scattering S- matrix element S nA s j. nA s x given by Eq. (|A43[) we obtain 



M. 



DW (prior) 



(P, k dA ) = i 



2n 



k n A RnA 



E 



i < I mi s m s I Jf Mf > < I m\> s m s > I Jf Mf > 



Jf Mf I m s t mi m^ M n 
X < J n M n J p M p \J d M d > < J n Mn J A M A \s m a > > Y{, '{-k n A) 



N 



l-e 



(l + » FvnAsljAEnA)} 1/2 {A-% T [r T nAsU F (EnA)} 1/2 ) O t {k nA , R nA ) 



„ (-) Of (k n A, r nA ) R n A \i a TT i (+K 

I nA \rfnA: n-nA) 



r nA >R„ 



(108) 



One-level, one channel case is the simplest one for which 

M DW{prior)^ ^ ^ ^ 
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M OW^ior) {p ^ dA)=i 



2tt 



k n A RnA 



E 



i < I mi s m s | Jp Mp > < I m\t s m s > \ Jp Mp > 



Jf Mp I m s i mi m ; / M n 



x < J n M n J p M p \J d M d >< J n M n J A M A \sm s , > Y t * {-k nA ) 



1 - e - 2lS *Asu F e 



2i<5„ 



(-) OfjknA, r nA ) 
x ^ Aj,_F 



Rn 



r nA Oi(k n A, RnA) 



Y*(v nA )\/\V dA W d X d + 2 > 



r nA >R nA 



Oi(k n A, RnA) 

(109) 



where 



JnA s I J F 



arctan 



nA s I J F 



(E nA ) 



2{E n A(o)si j F 



E 



nA(0) si J F 



> E n / 



(110) 



is the resonant phase shift, E n A(o) sij F is the real part of 
the complex resonance energy of the resonance with the 
quantum numbers s I J F in the channel n + A. Now we 



derive the equation for Mg by transforming it into the 
surface integrals over variable r n A- We can repeat the 
discussion in Section II A. The integration in Eq. (|9"9")l 
over r n A is taken over the external volume restricted by 
two spherical surfaces: the inner surface with the radius 
R nA and the external surface with the radius R nA — > oo . 
As it has been shown in Appendix [B] after regularization 
the integral over the infinitely large sphere vanishes (see 
Eq. dEp) and 



Mg (P, k, 



dA 



-M^(P,-k dA )=Rl A 



1 



2 [InA 

(+) 



o rry. (e;ct)( — ) / 

dr pF / dn rn AMr P n)xi + jAr dA )xit F M- nA 



i) 



r n A 



(+) , UT (ext)(-) rr w 9Mr P n)xl d ' A )(r d A) 

X-k pF \ r pF) [L nA \ v nA)\ JT- , 

U ' nA 



(111) 



Here, —MP W (P, k d A) is the surface integral encircling 
the inner surface of the external volume at r n A = R n A- 
A negative sign appears because the normal vector to the 
surface is directed to the center of the volume, i.e. op- 
posite to the normal vector to the external surface (at 
infinitely large radius). For simplicity, we dropped the 
quantum numbers in Eq. (|111|) but they will be recov- 



ered below. Note that Eq. (|111[) can be obtained from 
Eq. (|2T|) by substituting T^*^ _ - ) (r„ y i) for the overlap 
function I^(r nA ). 

For the exit channel c = b + B in the resonant sub- 
process n + A — > b + B different from channel c' = n + A 
using Eq. (|A37[) we get 



Ms° w (k pF , k dA ) = ~M™ (P, k dA ) 



2 7T IVbB „9 1 



R 



■bB V v nA 



nA 



'AWL 



E 



< I mi s m s \JF Mp > 



Jf Mf I I' mi m-i' s' M n 



x <l'm v s m s , \JM>< J n M n J A M A \ s m s > >< J„ M n J p M p \ J d M d > Y t * mi (-k bB ) s , v . bBsl J d r pF x { ^ pF ( r pF ) 

m v i- \\ , ^ (+)/ d ° l ' [k ;^ rnA) () .O v {k nA ,r nA ) dip d {r pn )x^ A ){*dA)- 

^T nA yi' mi ,{r nA ) (p d (r pn ) xl dA (rdA) jr^ X-J f {^ p f) 

l ■ o r n A 

27T IVbB 1 



r n A 



dr n A 



r„A=R r , 



kbB V V nA 2 \l n / 



y~] i l <lmi sm s \J F M > <l' m v s' m s >\J F M F > < J n M n J A M A \s' m s i > 



J F Mil' mi ra i' s ' M n 



x <J n M n J p M p \J d M d > Y l * mi {-k b B)S^ sl . nAs , ll O v (k nA , R nA ) / dr pF x { ^ pF (r pF ) / dCl InA Y v mi ,(r nA ) 



Vd(rpn) Xkl ( T dA) (BnA ~ I) ~ Rj, 



9^(r p „)Xkti)( r ^)l 



>r nA 



r nA = R nA 



(112) 
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Here. 



B n A — RnA 



dO t i (k n A,r„A) 




dr nA 


r n A=R,iA 



mal because M d and M p are fixed. The coefficient 
< J n M n J p M p \J d M d > appears from the vertex d — > 
p + n and the product < I' mi' s' m s >\JF Mf >< 
(f f 3) Jn M„ J a Ma I a' m s > > from the vertex n + A — > F. The 
matrix element S^£ sl As , v is given by Eq. (|A45|i . Sub- 



Ol'(k n A, RnA) 

is the boundary condition. Sum over M n is a for- stituting it into Eq. (HT2J) gives 
I 



M^{k pF , k dA ) = —Mg™(P, k dA ) = — ^ 

kbB V V nA VnA 



7T IVbB 1 



V] i l <lmi si7i s \JfMf > 

Jf Mf 1 1' s' mi m^ m s i M n 



x < I'mv s'm s ,\J F M F ><J n M n J A M A \s' m s , ><J n M n J p M p \J d M d > Y t * mi {-k bB ) e~ ' e"* 5 '"' 

N 

x E t r 

v bB s I Jf {E bB )} 1/2 [A- 1 ] 

vt [I-V nA s' V Jf m,i (^nA) 



<Pd{*pn) Xk dA ( v dA) {BnA ~ 1) ~ RnA 



dVd{r pn )xCX){ v dA) 



dr n A 



r nA = Rn 



(114) 



r 



Taking into account Eq. (|A41[) and Eq. (I104[) we arrive at the final form for Mg W {k pF , k dA ) 
I 



M° w (k pF , k dA ) = -Mjj? (P, k dA ) = 7T 



DW 



2 R nA 



HbB UnA kbB 



E 



% <lmi s m s \Jjr Mjr > 



Jf Mf 1 1' s' mi m t / m s t M n 

x < I' m v s' m s , | J F M F > < J n M n J A M A \ s' m s , >< J n M n J p M p \ J d M d > F z * mi (-k 6B ) e~ 1 ' 

^2 [ T vbBsl.J F ( E bB)} 1/2 [A^ 1 ]^ JmAs'l'.J / d Y pF X^k pF ( r pF ) / d ^r nA Y l' m,/ (? nA) 



<Pd(r pn ) Xktl ( TdA ) { B nA ~ 1) - RnA 



dVd(Tpn)Xkl)( r dA) 



dr n A 



r nA =R n 



(115) 



Now let us consider the diagonal transition csl csl, be replaced by ^cilm-ca im > + "^csT*^ ' 



where c = c' = n + A. To get Mg W once again we start 



m s ;csim s 



given by 

Eqs. (|A30[) and (|A33I) . Then the equation for the surface 
from Eq. (jlllj) . Now in this equation Y^*'*- - ' should matrix element for the diagonal transition takes the form 



Mg w (k pF , k dA )=i 



/J-nA k n A 



E 



i < I mi s m s | Jf Mf > < I rni>> s m s >> \ Jf Mf > 



Jf Mf I m/ m^// m s rr M n 

x < J n M n J a M A \s m s /r ><J n M n J p M p \J d M d > r z * mi (-k„ A ) 



iV 



1 -e 



-i2 5" 



v nA s I Jf {E nA )} 1/2 [A" 1 ] 

vr J- r nA s I Jf (£„A)] 1/2 ) \Ol{k nA , RnA) 



x / iVX-k sy (v) / dfi r „ A Yi mi „(f nA ) Vd(r p „)Xk^(r dA ) (B nA - 1) - i?„ A 



^yd(r p »)xLtI)( r ^) ' 



(116) 
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Summing up all three amplitudes M int (pos ' (P, k dA ), 

M°J {pn ° r \P, k dA ) and M° w (k pF ,k dA ) 
-M^ A (k pF ,k dA ) we get the total post DWBA 

for the (d,p) stripping. 



(i) Resonant reaction n + A -+ b + B, that is 
c = b + B ^ c' = n + A. Then the total post 
form of the DWBA deuteron stripping amplitude is 



M W "(P, k dA ) = 2nJ V- 



E 



i < s m s lmi\JF Mp > < s' m s > V vtiy \ Jp M f > 



Jf Mf s' I V m s i mi M n 

< J n M n J A M A \s'm s ,><J n M n J p M p \J d M d > e~ iS ™< ^(-kja) £ [T vhBs i j p (^6b)] 1/2 [A" 1 ], 



N 



x | <X ( pF } lAs>i>j F (r nA )\A V pF | xff > 
(-) Of'jknA, r nA ) R nA 



r nA <R nA 



X < 



X P F 



RnA 

(+) 



IAT— 1 



7r nA s' V J F 



r nA OT,(k nA , R Tl 



^ mi Xr nA )\AV pF \ip d xY A ' > 



r„A>Rn 



R 



11 A 



It nA s'l'j J d r pF x^ pF (i" p f) J d £l rnA Y v „ h , (r n/ 



<Pd(r pn ) Xu/A^dA) (BnA - 1) - RnA 



dip d (r pn )x < C,\)(r dA ) 



dr n/ 



r nA =R n 



(117) 



Assuming in this equation b = n and B = A, that is 
c = d but I ^ V and/or s ^ s' we get the expression 
for the DWBA deuteron stripping for the non-diagonal 
transition in the resonant subprocess (n + A)i s — > F — > 
(n + A)i> s f. 

Equation (|117[) is very instructive for understanding 
the difference between the stripping to resonance states 
and on-shell binary resonant reactions. As we can see, the 
transfer reaction amplitude contains the resonance fac- 
tors determining the resonant subprocess n + A — > b + B, 
the partial width amplitude [T^bBsij F {E bB )Y' 2 of the 
level v for the decay to the exit channel b + B, the matrix 
elements of the inverse i?-matrix level matrix [A _1 ]^ T 
and the reduced width amplitude j T nA s i \> j F of the level 
t for the entry channel n + A rather than the correspond- 
ing partial width amplitude which would present if we 
consider the corresponding on-shell binary resonant re- 
action n + A -+ b + B. The difference is crucial because 
the partial width amplitude \T u bBsi ,/ F (^bs)] 1 ^ 2 contains 
the penetrability factor, see Eq. (|A41|) . which doesn't 
present in the reduced width amplitude j T nA s' l' J F & n d, 
hence, in Eq. (|117[) . The lower is the energy of the reso- 
nance, the stronger is its suppression due to the barrier 
penetrability in the entrance channel in the on-shell bi- 



nary resonant reaction n + A — > b + B. Besides, if a few 
resonances do contribute with the different I', then the 
higher I', the stronger its suppression. However, it is not 
the case if one tries to populate low-energy resonances 
with different I' using transfer reaction. Missing pene- 
trability factor in the entry channel of the subresonancc 
reaction n + A — > b + B in the transfer amplitude makes 
it possible to populate low-lying resonances. Moreover 
for the same reason, the resonances with higher /' are 
not suppressed in the stripping. Hence, when a few res- 
onances are populated in the transfer reaction, the mea- 
sured experimental spectrum of the fragments b and B 
can be quite different from the one measured using the 
on-shell binary resonant reaction. The missing penetra- 
bility factor in the entry channel n + A of the resonant 
subreaction n + A — > b + B in the transfer reaction ex- 
plains the power of the Trojan Horse method as indirect 
technique in nuclear astrophysics (see [H, [H| and refer- 
ences therein). 

(ii) Diagonal transition in the resonant subprocess (n+ 
A)is -> F -> (n + A)i s , that is, c = d , I = I' , s = s'. 
The total post form of the deuteron stripping DWBA 
amplitude is 
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M DW( P ost) (p ; k dA ) = 2 7T il < s TO « 1 m i I J F M F ><s m s „ I m V > \ J F Mp > 

Jf Mp I m s n nil m ;" M n 

x < J n M n J A M A \sm s „ ><J n M n J p M p \J d M d > e -' lS ^< y* (-£„a) 



\ J ^7— Y,^^ sl j F (EnA)} 1/2 [A- 1 U <X { ; F ] I F A sl j F {r nA )\&V pF \v d X d + 2 > 

V VnA K n A y 
N 

i _ e -i 2 , ^ + , ^ [r, „^ s ; Jf (£„a)] 1/2 [A- 1 ],. r r nA , , Jf (£„a)] 1/2 ' 



r„A<-R„ 



■ (M^nAj tt-nA) < X pF 



Rn 



k n A RnA r nA 0*(k nA , R nA ) 

J dr ?FxL"k pf ,(r p F) J dn rnA Yi mill (r nA ) 



Y*{inA)\^V pF \ VdXd + 2 > 



r„A>RnA 



2 [InA 

' i \ (+) i wo -i \ D <9^(r p „)x^)( r ^) 

<fd(r p n) Xk dA (rrfAj (-BnA - 1) - 



r„,=,R 




(118) 



r 



B. Stripping to resonance states. Prior form of 
DWBA. 



Here we show that starting from the prior form we are 
able to obtain the generalized DWBA i?-matrix ampli- 
tude for the deuteron stripping to resonance states, Eq. 
(|96p . much easier than from the post form. The prior 
of the DWBA amplitude for deuteron stripping to reso- 
nance states is 



M DW(prior)^ j^j 



=< x pF 



(-) t(-)| I A T/ 



&V dA \(p d tp AXd + 2 > 



(119) 



where AV dA is defined by Eq. and =< 



(p A \t> 



(-) 

bB 



> 



As usually, we split the amplitude into 



internal and external parts 



M 



DW (prior) 



(P, k dA ) = M 



DW (prior) 



+M™ [pnor) 



(P, kdA) 

(P, k dA ) (120) 



with 



M DW( P r,or ){ ^ ^ 



X pF 



r nA <R n 



(121) 



and 



M, 



DW (prior) 



--<X^T[:f^\AV dA \^ d x d + 2> 



r nA >Rn 



(122) 



The splitting of the amplitude into the internal and ex- 
ternal parts in the subspace over the coordinate r nA is 
necessary to rewrite the prior DWBA amplitude in the 



generalized i?-matrix approach for stripping to resonance 
states. As we have discussed in subsections III Al and 



IIII A[ the external matrix element M 



DW (prior ) 



the 



prior form is small and in some cases, with reasonable 
choice of the channel radius R n A, even can be neglected. 
It is important for analysis of the stripping to resonance 
states because the external part in the post form doesn't 
converge. In this sense the usage of the prior form in 
the external part has clear benefit. The main contribu- 
tion to the prior form amplitude M DW ^ rior ^ comes from 
the internal part M®^(P rlor \ Since the internal part 
is given by the volume integral, its calculation requires 

(int)( — ) 

the knowledge of in the internal region. The 

model dependence of this function in the nuclear interior 
(r n A < RnA), where different coupled channels do con- 
tribute, brings one of the main problems and main uncer- 
tainty in the calculation of the internal matrix element. 
Using the surface integral we can rewrite the volume inte- 
gral of the internal matrix element in terms of the volume 
integral in the post form and dominant surface integral 
taken over the sphere at r nA = R nA . With reasonable 
choice of the channel radius R n A the contribution from 
the internal volume integral in the post form can be min- 
imized to make it significantly smaller than the surface 
matrix element. The latter can be expressed in terms of 
the i?-matrix parameters - the observable reduced width 
amplitude (ANC), boundary condition and channel ra- 
dius. Repeating the steps outlined in subsection III Bl we 
get 



M, 



DW (prior) 



(P, k dA ) 



„ r DW(post) 
= M int 



(P, k dA ) + M° W (P, kdA). 



(123) 
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Here, M^^^° at ^ has been previously considered and is 
given by Eqs. ([TOT) and ([TUT]) while Mg W takes the form 



{P, k d A 

(-) 



< X 



pF 



T (int)(-)i fc 



Trt I (+) . 



r nA <R„ 



(124) 



where T^"*)^ =< >• The fact that the 

volume integral in this equation is the internal one makes 
transformation of this volume matrix element to the sur- 
face one much easier than for the post form. The tran- 
sition operator T — T p p + T n A- Since r n A < R n A at 
r p p — > oo the integrand in Eq. (|99[) vanishes exponen- 
tially due to the presence of tp d . Hence, the operator T p p 
is Hermitian, that is, applying the integration by parts 
over t p f twice we get 



(-) T (int)(-)| 7f$ I (+) 

_. v H T M)Hi^ 7$ I (+) ^ 

- < XpF 1 nvi M pF - J pF I XdA > 



= 0. 



r nA <R nA 



r nA <R nA 

(125) 



Thus Mj? w {P, k dA ) reduces to 
"(P, k dA ) 



Mg W ' 



2* 



(126) 



Using the Green's theorem we can transform this vol- 
ume integral into the surface one. Note that the volume 
integral over r n A is constrained by the sphere with the 
radius r nA = R n A- Hence, only one surface integral ap- 
pears with r n A = R n A- Here we see an important ad- 
vantage of using the prior form versus the post one. In 
the post form transformation of the external volume in- 
tegral to the surface one led to two surface integrals at 
r n A = R n A and r n A = R n A - > oo. It required an elabo- 
rate proof, which included regularization, to demonstrate 



that the surface integral at r nA = R n 



oo vanishes. 



After transformation to the surface integral we get 
Mg W (P, k dA ) = -M™ (P, k dA ), 



(127) 



Eqs (jllip . (|112p and (|115l) determine this surface inte- 
gral. 



Stripping to resonance states. Post CDCC 
formalism. 



The CDCC approach for stripping to resonance states, 
which takes into account the deuteron breakup in the 
initial channel, definitely has advantage compared to a 
standard DWBA. The application of the surface formu- 
lation of the reaction theory for the DWBA has been 



done mainly for demonstration, but our main goal is the 
CDCC. 

Here we present the derivation of the post form CDCC 
amplitude using the surface integral formulation. This 
amplitude is 

M CDCC( P ost ){kpFj ^ 



-<x { Jri-]\Av:^? DCC ^> 



ipF x nA I " v pF 

This equation is an extension of the post CDCC ampli- 
tude for stripping to bound states, see Eq. (|54|) . obtained 
using replacement 1% — > T, lA \ A V pF " is defined by Eq. 
([TO)) . Now, as usually, we split M CDCC( - post ^ into the 
internal and external parts in the subspace r n A- 

M CDCC( P o S t ){Pi ^ = M CDCC(post)^ ^ 

+ Mg° cc{post \p,k dA ). (129) 



The internal amplitude j^ G J^ GG< J> ost ) j s gj ven by 



M. 



CDCC(post) 



(P k 



■dA) 



z < X P F 



(-) T (mt)(-)| 



^v P p7\^ DCC{+) > 



r nA <R nA 

(130) 



Correspondingly, the external amplitude is 



MZ? CC(p ° st \P k dA ) 



-< X 



pF nA 



pF 



,yCDCC(+) > 



r nA >R nA 



(131) 



Now we repeat the steps outlined in subsection lll CI Tak- 
ing into account Eqs (f58|) . ([59]) and (|60|) we arrive at 



M. 



CDCC { post ){p ^ kdA) = M CDCC( P ost) { ^ R 



dA) 



_ . _.(-) T («Et)(-)|fe _ 7^, V CDCC( + ) 
— X p F L nA I ± ± I *i 



r nA >R„ 



(132) 

where T = T p p + T n A- It is shown in Appendix [Cl that 

DC C (post) | ;)e recm ced 



j^CDCC{post) 
XpF L nA l 1 



{P, k dA ) 

" ~f nA\ ^' 



CDCC(post) 



M 

CDCC(+) 



(P k dA ) 



> 



r nA >R nA 

(133) 



This integral can be directly transformed into the surface 
integral with r n A = RnA encircling the internal volume, 
while the integral over r p p is taken over all the coordi- 
nate space. Thus we have shown that the post CDCC 
amplitude for stripping to resonance states is given by 
the difference of two terms, internal post CDCC ampli- 
tude and the surface integral: 



M CDCC(post)^ kdA) = M 



-M 



CDCC(post) (p , s 
int i K dA) 

CDCC(post), 



'(k pF ,k dA ). (134) 
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The internal amplitude ]\/[CDCC(post) can ^ e m j mm j zec j 
by a proper choice of U pF and the channel radius R nA , 
while the surface integral is dominant. If the channel 
radius is larger than the n — A nuclear interaction ra- 
dius the second term is parametrized in terms of the re- 
duced width amplitude and the boundary condition at 
T n A = RnA- Thus we succeeded to parametrize the post 
form of the CDCC amplitude in terms of the i?-matrix 
parameters. It is one of the main results of this paper. 
Eq. (|134[) is the most important result of this paper. 
Due to the absence of the external term, which is present 



in the DWBA and which causes the convergence issue, 
the convergence problem in the post CDCC approach is 
resolved: the integration in the surface matrix element 
is performed over the full coordinate space only over one 
coordinate r pF rather than over two coordinates, r pF and 
r n A- 

Expression for ^[Gi)CC(%iost) ^ different cases can 
be obtained from Eq. (jlOip by replacing the 
initial channel wave function fd{^pn) Xk~l)( r <w) by 
V i ; (r pF , r nA ): 



A , r CDCC(post) /r> , \ 2 7r k bB 



i < sm s I mi\Jp Mp > 



J f Mp s' I V m s i mi m,/ M n 

< s'm s , Vmv \J F M F > < J n M n J A M A \s> ' m s , > e -"«, Y^-k^) J2 l T »bB sij F {Eb B )] 1/2 [A" 1 ] 



N 



x < X P J Y; mi ,(r nA ) I* A s , , jF (r nA )\AV pF \^^^>(r P F, r nA ) > 



CDCC( + ), 



U,T=1 



r nA <R nA 



(135) 



Note that the CDCC wave function itself also depends on 
quantum numbers of p — n and d — A subsystems, which 
we don't specify here. It will be done in the following 
up paper where concrete calculations will be presented. 
Natural Jacobian variables for \$)Cdcc{+) arg r ^ 
r pn , but we use here another set of Jacobian variables, 



r pF and r„ 



A- 



To write down explicitly Mg^ )CCt ' post \'k pF ^ k^) in 
terms of the surface integral we can use Eq. (II 1 1 p re- 
placing the initial channel wave function by the CDCC 
one: 



M 



CDCC(post) 



RnA 



2fl 



ii A 



(V, ^a) = ~M^ DCC{post \k pF , k dA ) 



J dr pF xL + k pF (rpF) J dil r 



(r pF , r nA ) 



(ext)( — )* 



dr n/ 



■^{ext){ — )* 



-u 



CDCC{+) 

dr n/ 



(r pF , r nA ) 



(136) 



We can extend corresponding equations from subsection in the resonant subreaction d s 1 1' — > cs I, where c = b+B 
IIII Al bv replacing the initial channel wave function by the and d = n + A, we get from Eq (| 1 1 5[) 
CDCC one. In particular, for the nodiagonal transition 
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M 



, r CDCC(post) 



(k pF , k dA ) = — Af 



CDCC(post) 



(P, k dA ) = IT , 



2R. 



nA 



MfcB VnA k 



E 



» Jf Mf / / s m; m ( / m s / iVZ-, 

x < Imi s m s \ Jp Mp > < I' mi' s' m s > \Jp Mp > < </„ M n Af^|s m s ' > 

x y ; * m; (-k 6B ) e" 1 - ^ [T„ fcS , , Jf (E bB )] l ' 2 [A-\ T 7r nA s , v Jf J d r pF X % p (v pF ) J d Sl TnA Y v „ v (f„ A ) 

(137) 



*? DCCW (r pF , r nA ) (B nA - 1) - i?„, 4 



0<ltf DCC(+ W>'»A) 



dr nA 



r nA =R n 



r 



Correspondingly, the surface integral for the diagonal transition c si ^ c si can be obtained from Eq. (|116l) : 
I 



M 



CDCC(post) 
S 



(k pF , k dA ) 



VnA k nA 



E 



i < I mi s m s | Jp Mp > < / ra\n s m s n \ Jp Mp > 



x < J n M n J A M A \sm s „ > Y* (-k nA ) 



N 



1 - e~"" nAl 
dr pF X-L^pF ) / dn rnA Y lmi „(r nA ) 



''l +1 [ T -nAsljAEnA)] 1/2 [A- 1 ]urT TnA slJ F (EnAT l1 

1 LH 1 + '(r pF , r nA ) (B nA - 1) - R nA 



Ol(k n A, RnA) 

d^f DCC ^(r pF ,r nA ) 



dr nA 



r nA =R nA 

(138) 



I 

Summing up two amplitudes M^ l f CC ^ post \p, k dA ) (i) Resonant reaction n + A -» b + B, that is c = b + B ^ 
and Mj^kj^, k dA ) = -Mg™ (k pF , k dA ) we get the d = n + A. The total post form of the CDCC deuteron 
total post CDCC amplitude forthe (d,p) stripping. stripping amplitude can be obtained from Eq. {HZ}: 

I 



M CDCC(post)^ k ^ )(pj ^ = 2n 



1 



VbB kbB 



i l < s m s lmi\J F M F > < s m s / 1' mp \ Jp Mp > 

Jf Mf s' I V ra s t mi m^ M n 
N 



<J n M n J A M A \s'm s , > e " lS ^i Yl * mi (~-k bB ) £ [T v bB . I j F {E bB )] 1/2 [A" 1 ], 



x { < X^J I F As -u j F (r nA ) \AV pF \yf DCC{+ \r p p, v nA ) > 



r nA <R nA 



1 RnA 
2 VnA 



lr nA s< I'Jp J d V p p X^p ( r pF) J d f2 r „ A Y V m; , (v nA ) 

*? DCC ^(r pF , r nA ) (B nA - 1) - R nA ^£^!W^ 

O r nA 



r„A=Rn 



(139) 



Assuming in this equation b = n and B — A, that is c = d but I ^ V and/or s ^ s' we get the expression for 
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the post CDCC deuteron stripping for the non-diagonal (ii) Diagonal transition, c = c', I — I' , s = s'. The 

transition in the resonant subprocess (n + A)i s — > F — > total post form of the CDCC amplitude is 
{n + A) Va ,. 



M CDCC ^ ost \p : k dA ) = 2tt i l <sm s lmi\J F M F ><sm s „lmi„\J F M F > 

Jf Mf I m 8 n mi Ttii/t M n 

x < J n M n J A M A \sm s » > e~ iS ^ Y t * (-k nA ) 



\\l 1 , l r vnAsij F (E nA )} 1/2 [A^U < XpF } I F A sl j F (r nA ) | A V pF \9 

[ V Hn A K nA ' 

N 

v n A s I Jf ut T T n A s I Jf 



T7 ui,CDCC(+) 



(r pF , r nA ) > 



r nA <R„ 



l_ e -«0« 



l —— [ dr pF x^ pF (rpF) / dO rnA y imi „(f„A) 

^ ^n A K nA J J 

xSi cdcc(+) ( w„ i \ n> a*f I?C ' C(+) (r pi ,, r„ A ) 

* 4 l r P-F) r nAj ~ J- J ~ ^nA 



dr nA 



(140) 



r 



Eqs (|139p and (|140j) are the final main results of this 
paper. Both matrix elements consist of only two terms, 
the internal post CDCC and the surface term. The in- 
ternal term contains the integration over r nA in the in- 
ternal volume r nA < R nA . Hence, at r pF — ¥ 00 vari- 
ables Td A ~ T pF — > 00 and r pn ~ r pF —¥ 00. But 



then * 



CDCC(+) 



(rpF) r nA ) 



[361 ] and the integral 



over r p ^ does converge. The same conclusion is true 
for the surface integral in which r nA = R nA . Hence, 
in this matrix element also ^f DGG ^ + \r pF , r nA ) ~ r ~p 
and integral over r pF converges. Both amplitudes are 
parametrized in terms of the parameters used in the con- 
ventional i?-matrix approach and allows us to analyze 
the stripping into resonance states using generalized R- 
matrix approach. Finally, both amplitudes, (|139[) and 
(|140[) . don't have penetration factor in the entry channel 
n + A of the resonance formation in the resonant subre- 
actions n + A — >• b + B and n + A n + A. That is 
why stripping to resonant states provides a powerful tool 
to measure resonances in the subsystem n + A very close 
to the threshold, which can be suppressed in the on-shell 
binary resonance reactions but not in the stripping to 
resonance states. 



IV. SUMMARY 

The theory of the deuteron stripping populating bound 
and resonance states based on the surface integral for- 
malism is presented. To demonstrate the theory I first 
develop it for the DWBA. Since the DWBA is outdated 
and, definitely, deficient compared to the CDCC, the the- 
ory is extended to the CDCC formalism. The theory is 



applied for stripping to bound and resonance states. The 
eventual goal of this paper is to deliver the theory of the 
deuteron stripping to resonance states within the CDCC 
formalism using the surface integral formulation of the 
reaction theory [32] . Transformation of the volume ma- 
trix element to the surface one (in the subspace over r nA ) 
and i?-matrix representation of the scattering wave func- 
tion of the fragments formed by the resonance decay al- 
lows one to parametrize the reaction amplitude in terms 
of the i?-matrix parameters used in the analysis of the 
binary resonant reactions. Since the reaction under con- 
sideration is the deuteron stripping, the presence of the 
deuteron in the initial state and exiting proton causes 
the distortions. That is why the reaction amplitude, 
in addition to the i?-matrix parameters describing the 
binary subprocess, contains additional factors - CDCC 
wave function describing the d — A scattering in the ini- 
tial channel (coupled to the deuteron breakup channel) 
and the proton distorted wave in the final state. Hence, 
the approach can be called a generalized R matrix for the 
stripping to resonance states. The advantage of the ap- 
proach is that the reaction amplitude for stripping to res- 
onance states in the post CDCC formalism doesn't have 
convergence problem and is parametrized in terms of the 
same observables as binary resonant reactions. Hence, 
the formalism provides experimentalists a consistent tool 
to analyze binary resonant reactions and stripping re- 
actions populating resonant states extracting the same 
observable parameters, namely, reduced widths (ANCs). 
The power of the method has been demonstrated in the 
analysis of the Trojan Horse reaction 19 F(d, n a) 16 [33| . 
The numerical application of the method will be demon- 
strated in the following up papers. 
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Appendix A: b + B scattering wave function 



In this Appendix we consider the representation of the 
scattering wave function used in the i?-matrix ap- 
proach for binary resonance processes [3^, H(| 



We consider only two coupled channels: c = b + B 
and c' = n + A. Also X Jf Mp is an eigenfunction of 
the Hamiltonian describing the compound system F — 
n + A = b + B in the internal region excited to the dis- 
crete level r with the total angular momentum Jp and 
its projection Mf [42| . A separable form for tPcsm!'"" 1 "' 
reflects the fact that we consider the b + B interaction 
proceeding through resonance states. The entry channel 
of this scattering is the channel c = b + B. The inverse 
level matrix contains contribution from all N resonance 
levels. In a simple one level case it reduces to the well- 
known Breit-Wigner resonance propagator. All the open 
channels coupled to c contribute to X Jf Mf and deter- 
mine possible exit channel contributions into resonance 
scattering. Hence, in the internal region, where different 
open channels are coupled, X Jf Mf can be written as a 
nonorthogonal sum of these channels [39] : 



1. Internal scattering wave function "J" 



(+) 

bB 



A general equation for the internal wave function con- 
tains the sum over total angular momentum Jp and its 
projection Mp- Since we are interested in a wave func- 
tion ^ h + g describing a resonance in the system F = b+B, 
we consider only the internal wave function at given Jp, 
at which resonance occurs. In the internal region in the 
state with the total momentum Jp, channel spin s (its 
projection m s ) in the initial channel c = b + B the wave 



function vE^g can be written as [3£ 



X 



J F M F _ 



E 

c s I ms 



l 



— W T r, 

r, 



TCJ 



,J F M F . 
cslms csljp j 



}, 

(A4) 



where £g is the product of the internal bound state wave 
functions of the fragments in the channel c, c = c, c', 
(rg) is a set of the radial wave functions of the 



M; - 



c s I J F j 



relative motion of the fragments in the channel c with 
the channel spin s, orbital angular momentum I and to- 
tal angular momentum Jp in some adopted potential, 
is the channel spin-angular wave function (in LS- 



L J F M F 
ill 



i uJ F (int){+) = ^ kc 
csm = k \ ii 



E 

Ml mi 



N 

E 

ZV,T = 1 



il/2 



Here, E c = E^b and k c = k^e are the relative energy and 
momentum of particles b and B, [i c = ^b , ^vcsij F (E c ) 
is the formal (i?-matrix) partial resonance width of the 
level v in the channel csl Jf, c = b + B, A is the R- 
matrix level matrix, N is the number of the levels in- 
cluded, 8^1 is the hard-sphere scattering phase shift in 
the channel c given by 



coupling) and mg is the projection of s. Also A is the 
antisymmetrization operator between the nucleons of the 
< s m s I mi | Jf Mp fragments in the channel c. We consider only two cou- 
pled channels, c = b + B and c! = n + A. Thus the initial 
channel c can propagate into two final channels c and 
] vt X' JfMf c' via the intermediate resonances. Although Eq. (IA4|) 

contains the sum over all channel spins s and its projec- 
(Al) tions in each open channel, in what follows consider the 
contribution to X^ F Mf only from the channel with fixed 
channel spin and its projection. 

First, let us consider the contribution of the channel 



„// into X^ F Mf . In this channel £ c = ipj, (fs and 



r ,T m = E < s " ms " l"mi"\J F M F > 



X Yl" m,// ( r c) 4>c s" m B n i 



(A5) 



<5' l f = — u>d + arctan 



Fi(k c , R c ) 
Gi(k c , R c )' 



(A2) 



where Fi(k c , r c ) and Gi(k c , r c ) are regular and singular 
Coulomb solutions of the radial Schrodinger equation, 



t>cs"m s „= E < Jb Ah Jb Mb\s" m s » > 

M b M B 



X 1pJ b M b lpJ B M B 



(A6) 



Ul = G c l — <J c0 



E^ arctan — , 



(A3) 



<j c i is the Coulomb scattering phase shift in the channel 
c and in the partial wave I, r] c is the Coulomb parameter 
for the scattering of the fragments in the channel c. 



Here, 4> CS " m ,, is the channel spin wave function in the 
channel c s" m s u , ipj t m { is the spin wave function of 
particle i, I" (mi") is the relative orbital angular mo- 
mentum (its projection) of the fragments in the channel 
c, r c = TbB is the radius- vector connecting b and the 
center-of-mass of B. We adopt the channel radius R c 
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large enough to neglect antisymmetrization between nu- 
cleons of b and B at r c = R c , that is 



A{& - kJ " Mf 

lJf M f 



cs" I" m B n u cs" I" J F j} 
1/2 



r c =R c 



r c =R c 



(A7) 



where JV c =(^il)" 

Assuming that the overlap of the channel c at the chan- 
nel radius R c with the channel c' is negligible we get for 
the component of X^^^ ;/ projected on £ c = ipbifiB at 
r c = RcM 



r c =R c 



(A8) 



where 



ires",/" J F (^ C ) = Wrq «cs",i" J F j(r c )- (A9) 



At r c = i? c by definition [3£ 



-cs"i"^(-Rc) = V 2 ^c RcItcs" i»j f , (A10) 



where 7 TCS » /" ,/ F is the reduced width amplitude of the 
level r in the channel c s" I" Jp. I remind that the system 
of units H = c = 1 is being used throughout the paper if 
not specified otherwise. Then 



iJ F M f 
c s" I" m a 



(All) 



Thus we can express the component Sj^^/^ „{ r c) 
taken at the channel radius r c — Rc in terms of the 
sum of the reduced width amplitudes, where the sum is 
taken over all allowed in the channel c partial waves I" 
at given Jp and s". Then the component of $csmf" + ' 
in the exit channel cs" m s » projected onto £ c = ipb ^>b at 
r r = R r takes the form 



~,JFyini)t,-t-) ,r> ~ \ _ < t | x jyJ F (mt)(+) 
c s m s ; c s" m // V ■'C c ) *^ Sc I c s m s ; c s" m // 



,J F (m*)(+) 
m s ; cs" r 

2 7T I k r 



kc Rc y 



E 

M F , I mi 



N 



< sm s l mi \J F M F > y; mi (k c ) [^csUf(Ec)} 1/2 [A-% T E'^ ms JR c r c ) 



2tt 



kc Rc 



^2 e~ i5 "' i l < sm s lmi\J F M F > < s" m s ,, 1" mi-,\J F M F > Y* mi (k c 



E Fvcs i j f (E c )} 1/2 [A- 1 ],, 7r c ,« ,« Jf y r m; „ (f c ) , 



AT 



(A12) 



r 



Here, s" is any channel spin value in the channel c = b+B 
allowed by the spin and angular momentum conservation 
law. In particular, s" may coincide with s, that is s" = s. 



The diagonal component, I" — I and s" = s, which is 
needed to determine the elastic scattering amplitude (see 
below) is 



* J>(mt)(+) 
' cs I m s ;c si m 

N 



„(i? c f c ) = 2tt 



e- i5 *> i l 



< s m s lmi\J F M F > < s m s » I mi>i \ Jp Mp > Y{" mi (k c ) 

Mf mi m L n 



X 



E c I Jf (Ec)] 1/2 [^Xrlrcsl J F V m ,„ (? 



(A13) 



I/.T— 1 



A similar consideration can be applied when we con- sider the contribution of the channel c f s f m s > , where 
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d = n + A, into X J T F M " . In this channel £ c / = ip A 
and 

<f> J /f/ m „ = E < s ' m *' l'mv\JFM F > 

m v 

xY Vmi ,{r c ,)<i> cl s > m3 ,, (A14) 



t>d s'm,, = E < J n M n J A M A \s'm s * > 

M n M A 



X 1pJ n M n i>J A M A ■ 



(A15) 



Here, <j> c i s < m , is the channel spin wave function in the 
channel d with the channel spin s' and its projection m s i , 
I' (mi') is the relative orbital angular momentum (its 
projection) of the fragments in the channel d, r c > = r nA 
is the radius-vector connecting n and the center-of-mass 
of A. We adopt the channel radius Rd large enough to 
neglect antisymmetrization between n and nucleons of A 



at r c i — R c >, that is 



Mi 



,J F M F 
d s f V r 



Uc's'V J F j} 

J F M F 



where N c , = 



u d s' I' Jf j 
-1/2 



r r t=R r 



(A16) 



= (A + l)-V2. 

Assuming that the overlap of the channel d at the 
channel radius R c < with the channel c is negligible we get 
for the component of X^^, F m i projected onto £ c / = ip A 



at r c i = R c > 

= E tiffim., "rC s' I' J F (Rc'), 



(A17) 



where 



c' »' I' 



J F ( r d ) = N c , E w r c' J «c' s' /' Jf j ( r d ) ■ (A18) 



At r c < = i? c / 

U T d s> V J F {R& ) = a/2/V Rc> It d s> V J F , ( A19) 

where /v = fJ, nA , Jtc's'I'Jf i s the reduced width am- 
plitude of the level r in the channel d s' I' Jp. Then 



(A20) 



that is it can be expressed in terms of the sum of the 
reduced widths amplitudes in all allowed partial waves I' 
in the channel d at given Jp and s' . Then the component 
^csnT^dV'm , projected on £ c / = (p A at r d = R c > takes 



the form 



yJ F (mt)(+) (R , r ,) 
c s m s \ d s' m s i V c C J 



2tt k c 



N 



kc Rc r y l^c 



M F , I mi 



2tt 



2/x c 



E e ld ^ 1 i l < sm s lmi\J F M F > < s' m s > I' m v \ J F M F > 



\x c k c R c > 

M H' mi m l i 
N 

x r ; * mi (k c ) E ;j F (^c)] 1/2 [A- 1 ]^7rc' S '/'7 F ^ m! ,(?c)' 

i/.r—l 

i 

The component T^^+l {Rd f c ') is given by 



(A21) 



,J F (mt)(+) 



2^c' 



e iS d. i l ^ <sm s lmi\J F Mp><s' m s > I' m v \J F Mp> 

M F mi m i' 



N 



xr; mi (k c ) E [TvcsiJ F {E c )] 1/2 [ATX^rd s > v j F Yt, mi ,{vd)< 

I/,T— 1 



>c' s' m a / • 



(A22) 
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2. External scattering wave function *t,B 



Now we proceed to the expression for the \E'i + ' ) in the 
external region, where r c > R c or r c ' > Rc' ■ In the 
external region the wave function ^csmi with fixed 
channel spin and its projection in the incident channel c 
can be written as 

tf,(ext)(+) jrf(ext)(0) , tf,(ext)(+) f Aoo"| 

csm s csm s ' csm s ;r' y^^u j 

where the first term is the incident wave and the second 
term is the sum of the outgoing waves in all the open 
channels. The incident term is 

* M 0) = 4 7T £ c E i l <sm s lmi\J F M F > 

Jf Mf I mi m 8 ti 

Fi(k c , r c ) 



x < sm s n lmi\J F M F > Y i * mi (k c )e I 

xY lmi (t c )<j> csmB „, (A24) 

where the subscript c means that the incident wave is in 
the channel c. The sum over m s » is a formal because 

< s m s I mi \J F M F > < s m s " lmi\J F M F > 

Jf M f 

— Sm s m s n ■ (A25) 

Note that here we use the incident wave with the unit 
amplitude rather than with the unit flux density. The 

component ^{^.^j m „ , which corresponds to the exit 
channel c s I m s >> and fixed J F , projected on £ c reduces to 

M mi 

x < s m s Imi \ J F Mp > < s m s » I mi\Jp Mp > Y* m (k c ) 



x e 



^^%^Y lmi (r c )^ 



Now we take into account that 

e iu "Oi(k c , r c ) - e-^> Ii(k c , r c ) 



Fi(k c , r c ) 



2i 



(A26) 



(A27) 



Here, Oi(k c ,r c ) and Ii(k c ,r c ) are the Coulomb Jost sin- 
gular solution of the Schrodinger equation with outgoing 
and ingoing asymptotic behavior (we follow the defini- 
tions used in 39]): 

Oi{k c ,r c ) rc Z°° e^^"-^ in(2*.r„)-W2+*„o] j ( A2g ) 
and 

h{k c ,r c ) rc 3°° e -^ c - Vc M2 fce r c ) 1 -iV2+.co]. (A29) 



Then we can rewrite T 



J F (ext)(0) 
csl m s ;c s I m s 
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in the form 



T J F (ext)(0) , s _ . . 

1 cslm B ;cslm s „\ L c) — 1 , l 
r*>e 1 c 

x < s m s I mi \Jp Mp > < s m s n lmi\Jp Mp > 

M mi 

x F ; * mi (k c ) [l,{kc, r c ) - e l2 ^'Oi(k c , r c )] Y lrril (r c )0 csm 

(A30) 



Thus the incident wave is the pure Coulomb scattering 
wave function in the incident channel c. The second term 
in Eq. (|A23|) is given by the sum of the outgoing waves 
in the open channels |40l |: 



csm s ;r j n / j 



Vc r~ c 



E 



Jf Mf 1 1 mj 



x < sm s lmi\J F Mp > F/* ro! (k c ) \e l2u "' 5 c - c S Ss 6 Il 



S 



Jf_ 

cshcsl 



x <f>- c 



c s mg • 



Oj(fcg, rg) < sm s lmj\J F M F > F|- m _(fg) 

(A31) 



Here, £g is the product of the bound state wave func- 
tions in the channel c — c, d ', S Jf ~ is the S'-matrix 

c s l; c s l 

element for transition csl — > csl. Note that we consider 
the outgoing waves in the channel with given total angu- 
lar momentum J Fl initial channel spin s (its projection 
m s ) and final channel spin s (its projection m$). Since 
only two open channels, c and c', are taken into account 
here, we will write explicitly the outgoing waves in both 
channels. First consider the elastic scattering, that is 
the outgoing channel c = c = b + B and the channel spin 
and orbital angular momentum coincide with the incident 
channel values, that is s = s and 1 = 1. The component 
of the outgoing elastic scattered wave (csl — > csl) is 



(ext)(+) 

csl m s ;cs I m" s 



■ 2n e 

?c 



E 

Jf Mp mi m^, 



x < s m s Ityii\Jf > < s m s I m/// 1 Jf Mp > 

x ^* mi (k c )[e i2 ^« - S J / al . csl ]Oi(k c , r c )Y lmi „{v c ) 

x 4> csm: . (A32) 



Hence, the projection of ^ , 



c s I m s : c s 



, „ on £ c leads to 



c s I m s :c s lm 



E 



Jf Mf mi mi>i 

x < s m s lmi\Jp Mp > < s m s I mi" \ Jp Mp > 
x S Y* mi (k c ) [e* 2 - 5 C J - ;c s ,] Oi (k c , r c ) ^ m; „ (f c ) 
x c , m «. (A33) 



Correspondingly, for the inelastic scattering, c = c but 
either s ^ s or / 7^ I or both differ from the entry values, 
we get 



(ext)(+) 
c s i m s ; c 5'' 



■ 2lT t 



E 



Jjr Mf mi m^/ 

x < s m s lmi\Jp Mp > < s m S " l" mi"\Jp Mp > 
x i l Y? mi fc)Sil, v ,. cal O l ,,{k e , r c )y,» TOl/ ,(f c )0 c 



(A34) 
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Then the projection of ^cslm\s" l" 



on £ c is 



c s I m s \ c s" I" ?n s // y c > 



2tt 



E 



Jj? Afp m; m l i, 



x < s m s lmi\JF Mp > < s m s >> I mi" \ Jp Mp > 

x i l Y? mi fe)S%, l „. cal O v ,{k , r c )Y v , mi „{r c )<j> cs „ mBl ,. 

(A35) 

Finally, for the outgoing scattered wave in the reaction 
channel c = d = n + A we have 



(ext)(+) 

csl m s ■ c' s' I' m B 



r^c I c' \ u c' j w 



(A36) 



x < s m s Z m; | Jp Mp > < s' to s ' Z' mi' \ Jp Mp > 
x YC mi {K)S J / slll . csl O v {k c r, r & )Y Vmi ,{i c ')4> c 

It leads to its projection on £ c /: 

(ext)(+) ( Tl \ = _i2jL [?£aI V 

c s L m s ; c' s' l' m / \ c J b t \ 11 ' J 

1 J f Mp nil 

x < s m s lmi\JF Mp > < s m s > I' m\> \ Jp Mp > 

x Y l * mi {K)S J / s , ll . csl O v {k cU r c ,)Y llmil (r c ,)0 c , s , ms ,. 

(A37) 

Now we can derive the expression for the matrix ele- 
ments of the 5* matrix. Since the wave function \&c is 
continuous using Eqs. (|A13|) . (|A30[) and (|A33|) we get 
the equality 



csl m.icslm.,, K-^d-c) — J- 



J F (ext)(0) ( f? r } 

cslm s ; cslrn,,, kP-c l c) 



T 



(ext)(+) 



- cslm s ;cslm s /i Fc )' 

which boils down to 

N 



(A38) 



-1 



[r„ cs ;,/ F (£ c )] 1/2 [A-\ T ^2k~R clTCsljF 



V.T—1 



= i [h{k c , R c ) - .s/iv. ^o, .:/.-,. R c )] . 

Taking into account that [39| 

Ijjkc, Rc) = Gi{k c , R c )-iFi{k Cl R c ) l2LOci = 
Oi{k c , R c ) G^kc, R c )+iFi(k c , R c ) 6 

and 



(A39) 



-2i«? 



r rCS / J F (Ec) - 2Pcl(E c , Rc)ll cs lJ F : 



where 



P cl (E c , R c )- ^ Rc) C +Gf(k c , R c ) 



(A40) 
(A41) 

(A42) 



is the Coulomb-centrifugal barrier penetrability, we get 
the elastic scattering S-matrix element: 



JY 



qJf 

csl\csl 



„-2i<5 



X [Trc.u,^)] 1 ^). 



(A43) 



From equality of Eqs. (|A22[) and (|A37[) at given Jp and 

r c > = R d 

y,J F (int) f jj ~ \ _ r T J F (ext)(+) fjj ~ \ 

1 csm s ;c' s' m 3 i K^ ' 1 c ' > ~ ^slm.jc's'I'm,, k""*' 1 c ' 7 

(A44) 

we obtain the reaction matrix element: 

N 



/2 



I/,T=1 



x [A- 1 ],, [r TC ^ FjF (£; c ,)] 



1/2 



(A45) 



Both obtained matrix elements coincide with the corre- 
sponding matrix elements from (40j . The only difference 
is in the definition of the solid scattering phase shifts. 
The obtained matrix elements of the S matrix confirm 
that the relative normalization of the internal and exter- 
nal wave parts of are correct and we can use them 
to calculate the reaction amplitude of the deuteron strip- 
ping proceeding through resonance states. 



Appendix B: Matrix element Mg W 

Let us consider the DWBA surface (in the subspace 
over r n A) matrix element, which appears in the post form 
(see subsection IIII A[) : 

M° W (P, k dA ) =< *y T^f «| f 2* i 



x <PdX$ > 



pF 

r„ A >R n A 



= M^ F) (P, k dA ) + M™ A) (k pF , k dA ), (Bl) 



where T 



'S(pF)V 

(ext)(-) 
nA 



(ext)(-) 
bB 



>, 



M° ( J F) (P,k dA )= J dr^jdvx^tv)! 

T n A>RnA 

x [T pF - ~f pF ]ip d {v pn )x d + 2^dA) > 



(ext)(-)*/ \ 
A \YnA) 



(B2) 



and 



M™ A) (P, k dA ) = J Ar pF J Ar nA Xp p h (r pF ) T^^*^) 

r nA >R n A 

x \ f~ nA - f nA ]Mr P n) X dA {rdA) > ■ (B3) 
^SCpF) can ^ e wr itten as 
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M° p v F) (P,k dA ) = J dr nA J dr pFX y*(r pf )T^ ) *(r„ A )[^ pi ,-7 t pF ]^(r p „) X | z + ) (r £iA ) 
- J dr nA y , dr pf x pF ) *(r p f)T| 1 ™ t)( ~ ) *(r„ A )[^ pF -^ pF ]^ d (r p „)xi^ ) (r ( iA) > 

r nA <R nA 

= J dr nA J dr^xy*^)^^^^)^-^]^^)^^^^). (B4) 



We took into account that for any finite volume r nA < 
R nA the matrix element containing ^T pF — 7^ pF vanishes 
as it has been discussed in Section II A for deuteron strip- 
ping to bound states. To estimate Mg^^ we need equa- 
tions connecting different variables: 

r dA = 1/2 v pn + r nA , (B5) 
r pF = A/ (A + 1) (B6) 



Replacing the variable r nA by r p „ and transforming the 
integral over r pF is transformed to the surface integral we 
get for for the matrix element 



M 



DW 



(P, k dA ) = ( 



A + l 



9 Xp F *( r pF) 



) 3 lim R 

R pF -^oo 



pF 



dr pF 



dr pn ^ d (r p „)T^f )( ) *(^— i r vF ~ r P n]) X d + 2 



dr pF 



A 

(-U,A+ 1, 



A 



A + 2 . 



dr p „ <p d (r pn ) T^jf )( >(2-±±[ TpF - r pn ]) X <+)(_±_1 r pF - r p „) 



r pF =R pF ^oc 

(B7) 



Due to the presence of the deuteron bound state wave 
function the integration over r pn is limited. At r pF — > 
co and r pn < co we can replace the distorted waves in 
the initial and final channels by their leading asymptotic 
terms: 



(+) / \ r dA - 
XdA ( r dA) ~ 



a ik dA r dA +ir) dA \n(k dA r dA -\n dA -r dA ) 



(B8) 



and 



y{— )* r pF^°° g-ikpF-rpF + i i) pF ln(fc pF r pF +k pF -r pF ) (B9) 



Here, 77^ is the Coulomb parameter of particles i and j 
in the continuum. Note that r dA = - r pF — 4^ r p „, 
and at r pF — > co and r pn < co — > co. Then 



Q e ik dA -r dA +ir] dA \n(k dA r dA -k dA -r dA ) 



d 



r p F 



r pF — ¥oo 



_^ j Jj^^ . £ gikdA-C^x 1 • r pf-4i 2 r P>i)+i»?£i.4 ln(fc dA w-kdATdA) (BIO) 

I 

and For T^*^~^*(^^[r pF — r p „]) we can take only the exter- 

, , . nal part, which contains the resonant S matrix element, 

Q -ik pF -r pF +ir] pF \n(k pF r pF +k pF -r pF ) £L . . . . ' 

see hq. (|A37|) . Neglecting all the spin-dependent and 

dr pF 

r pF ->oo _ j. e _ikp F -rp F +i?)p F ln(fc pF r pF +k pF -r pF ) 

(Bll) 



33 



angular parts and leaving only its radial part, which is 
O n A ( r nA ) / r nA , we get for its leading asymptotic term: 

OnA{k„A, T nA ) r p F->oo A 1 ^ A+l ^ rpF ^ knA i- pF . rpr 

r n A A + 1 r pF 

-i[lnA ln(2k nA r nA )+l„ A Tr/2-cr nA0 ] 



element over an infinitesimal bin in the momentum plane 

fDW . 
l S(pF)- 



leads to disappearance of AlP n 



x e 



(B12) 



The leading term of its derivative at r pF — > oo is 

dO nA (k nA , r n A)/r n A r pF ^oo . 1 
— > l K n A 

or pF r pF 



— (fenA r pF -k nA PpJ? •!■„„) 
i[77»A ln(2^ife„ A r pF )+/„ A 7r/2-cr„ A o] 



X g .l. M »rT"M'P^'-«»"r "^Uj^ (B13) 

We also need the asymptotic behavior of the plane wave 



iq-r pF 'p- 



27T 



iqr pF 

e -iir pF S(q + r pF ) 



[e^ r " F 5(q-r pF ) 



(B14) 



where q = k^A — k p j? . Then the matrix element 
reduces to 



M^ F) (P, k dA ) ~ lim fi(Rpp) e*««p' 



+ f2(R P F)e- iqR " 



(B15) 



Thus has no limit at P p f ~> oo but regular- 

ization of this matrix element by integrating the matrix 



q+e 

±- ( dqM™ F) (P,k dA )^ lim 



sin (e P p f) 



x e 



/i(iZpF) - e-***" h(R pF )] = 0, (B16) 



where e << q. 

Similarly we can estimate Mg^ A ^ given by Eq. (IB3|) . 
Since the integral over r n A is taken over external volume 
with r n A > R n A the transformation of the volume inte- 
gral into the surface one leads to two surface integrals: 

M™ A] ( Pl k dA ) = -M™ a (P, k dA ) + M™(P, k dA ). 

(B17) 

The first term is the surface integral encircling the in- 
ner surface of the external volume at r nA = R nAl while 
the second term is the surface integral taken at r nA = 
R nA — ► oo. A negative sign in front of the first term 
appears because the normal to the surface is directed in- 
ward to the center of the volume, that is opposite to the 
normal to the external surface (at infinitely large radius). 
The surface integral over the infinitely large sphere in the 
subspace over r nA is 



M£f (P, k dA ) = - hm R nAT j— 



dn dT^ ^4) 



dr n/ 



dr pn tp d {r pn ) X pF *(-£—j *nA + r pn ) X d + A {^ r P n + r nA ) 



a!sl T nA L n A 



dr nJ 



(iwOtt^— I dr pn ip d {r pn ) xij* (-r^— r n A + r P n) X d + A (^ r pn +?ua) 



1 p f y A + 1 VnA Vpn ' XdA ^2 Vpn VnAj 



r nA =R ,^ca 



(B18) 



Here, the Jacobian variable r pF is replaced by v pn . Due 
to the presence of the deuteron bound state wave func- 
tion the integration over r pn is limited. Hence, r pF — > oo 
and r dA — > oo at r nA — > oo. At r nA — > oo and r pn < oo 
we can replace the distorted waves in the initial and final 
channels by their leading asymptotic terms. The dis- 
appearance of the matrix element (|B18|) can be proved 
similarly to the proof of the disappearance of Mg^ F y 
Replacing the distorted waves by their leading asymp- 
totic terms (|B8|) and (|B9[) , singling out the plane wave 
containing r nA and using the asymptotic representation 
of this plane wave, see Eq. (IB14[) . integrating over Sl r „ A 



we eventually arrive at 

M£f(P, k dA )~ ,11m [e^' 9l (R nA ) 
+ e-^' g 2 {R nA )]. (B19) 

Regularization of this matrix element by integrating it 
over an infinitesimal bin in the momentum plane q' leads 
to disappearance of Mg W , that is 

Ms™A)(P> ^a) = -M™ a (P, k dA ). (B20) 



34 



Appendix C: Matrix element M g DCC(post) (P, k dA ) 



can be written as 



Here we show how to transform j\^^ DCC (p° st ) \ n ^ th e 
surface integral over the coordinate r nA . j\/[C DCC (p° st ) 



M^ CC (P, k dA ) = J dr nA J dr^x^(r^)Ti^(-)*(r nA )[^-^]^ CC W(r pF , r nA ) 



r nA >R n 

= Mg° cc (P, k dA ) - Mg™ c (P, k dA ), 



(CI) 



r 



where 



Mg D ot cc {P, k dA ) = J dr nA J dr pFX ^Jr pF ) 
x r ( -/(*nA)fr fM DCC{+ \r pF , v nA ) (C2) 



and 



M 



CDCC 
Sint 



(P,k dA ) = J dr nA J dY pF x { ^ pF {Y pF ) 

r nA <R nA 

T£f <->*(r„A) & - 2*] *f DCCW (r pF , v nA ). (C3) 



space while in Mg^f c the external region in the sub- 
space over r nA is excluded. Let us first consider Mg^f . 
The CDCC wave function is given by Eq. If 
we substitute the first term, n = 0, which contains the 
deuteron bound state wave function, the transformation 
leads to the surface integrals with r pF — R pF — > oo and 
r nA — R nA — ¥ oo. Both surface integrals vanish and 
the proof is similar to the one presented in the previ- 
ous section. For the rest of the CDCC wave function 
corresponding to the sum with n > 0, which we call 



CDCC(+) 



, transformation to the surface integrals gives 



Note that in the matrix element Mg®^ c the integra- 
tion is carried over r pF and r nA in all the coordinate 



M 



CDCC(c) 
Stot 



{P,k dA ) = J dv nA J dv pFX % F {v pF )T { -}\v nA )[^ pF -^ pF ]^ CC{+ \v pF ,v nA 



j dv nA J dvxit(v)Ti/(r„A)[r„A-^i]$ 



CDCC(+), n 



= lim 



R2 pF 



RpF— >oo 2 )J>p F 



dn r 



dr nj 



■i ii ,, T ,CDCC(+), A ^X-k pP ( r pP) 



L nA {TnA) 



'{r pF , Y nA ) ■ 



dr pF 



a iT ,CDCC(+) / s 







r pF 



r>2 

lim nA 



dO r 



-X%Sr vF )^->{v nA ) 



dr 



pF 



v (+) (r v\^ CDCC{+) (r * r A \ dT nA*{?nA) 

X-\ L „\ r pF) V ic (TpF, T n A) 



dr„ 



dr nA 



(C4) 



Let us, first, consider the first term, in which R pF — > oo. behavior of ^ff CC ^ + \r pFl r nA ) ~ r p j? and Eq. (|B14|) 

Let us divide the integration region over r nA into the we get that the first term goes to zer0 as R -2 In 

region r nA /R pF -> and the region where r„ A > i? pf -> the rema ined region r nA - i? pf -> oo and we consider it 

oo. In the first region we get that r dA ~ R pF -> oo and later xhe second term of Eq (JC4J), in which i? nA -> oo, 
r pn ~ i? p i? — > oo. Taking into account the asymptotic 
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we also separate into two regions: r p F/R„ 



and 



r pF ^ RnA 

and r dA ~ RnA 



oo. In the first region 
ooand ^?? CC{+) ( 



pn 



00 
-3 



(r pF , r„ A ) - , nA . 
Hence the matrix element goes to zero as R~ A — > 0. To 
consider the behavior of the first and second terms of 



Eq. (|C4j) in the second regions, where r nA , r p p —> oo, 
it is more convenient to introduce the hyper-spherical 
coordinates in the six-dimensional hyper-space: 



Here, m is the scaling mass parameter, for example, 
the nucleon mass. Then Mg t ®f c in the region, where 
r nA , t p f — > oo, can be written as the integral over the 
hypersphere encircling the volume integral with the ra- 
dius of the hyper-sphere p — > oo 32.]: 



PnA __2 

nA 



PpF 



'pF> 



r n A = P 

" PnA 

< a < tt/2. 



sm a, 



r p F = P 



PpF 



cos a, 



(C5) 



M 



CDCC 
Stot 



rn 



lim p 5 



2 (p n APpF) 3/2 P^°° 



v (-) Vr ^T^Wr „1 JLia ODCC M( r « r 1 
XpF l> r pFj J- n A (r p F, r n ^j 



dr 



pF 



7r/2 

r 

df raj 4 / da sin 2 a cos 2 i 


dp 



CDCC(+), x " (-)'/ \T ( -)*(r- 



r 



(C6) 



Here, hyper-radius p is the parameter going to infinity. 
The integrand contains highly oscillating (actually in- 
finitely oscillating) functions. The behavior of the inte- 
gral at p — > oo depends on the asymptotic behavior of the 
integrand. The integration over dr p F can be performed 

directly using the asymptotic form of X p f( t p f )- It is 
given by the Coulomb distorted plane wave, but for sim- 



plicity, what does not affect the final result, we neglect, 
as in the previous section, the Coulomb effects. Then 
the asymptotic form of the plane wave is given by Eq. 
(|B14p and, hence, integration over i p p using (5-functions 
is trivial leading to r p p = ±k pF . After performing the 
integration over <H p f only two integrals are left. From 
Eqs (dU, dHJ) and (Ell we get for 



r 2 

' pF 



A 
A+l 



r pF r nA 



A 2 

JaTT) 



2 nA 



/ 777 A I 777 / 777 A 2 777 

P \ / cos a T ~A 7 \ / \ / z sin 2 a + -r-. -r^ 

y PpF A + 1 y p p F V PnA (A+iy p nA 

I 



sm a 



(C7) 



and 



1 2 A + 2 (A + 2) 2 . 

UA = d 2 r vF + r pF r nA + r 



2 nA 



11 777 „ A + 2 rUT HrrT . n (A + 2) 2 m . , 
— cos^ a ± — — , / . / z srn 2 a + tti 7-77 sm a . (Co) 

4 Hp F 4(A + 1) V /XpF V PnA i(A + l) 2 /i„A 



r 



Here, z — r nA -k pF . We recall also that in Eq. (|52|) at 77 > e r |" pn , while x- (r<&i) ~ e ' ^ ' dj4 , where we, for 
^pn ( r pn) at r p „ — » 00 contains the asymptotic terms 



3G 



simplicity, neglected the Coulomb distortion. Then after 
integration over df pF the leading asymptotic form of the 
integrand with omitted Coulomb effects is a product of 
highly oscillating at p — > oo exponents: 



Let us first consider the first matrix element containing 
T p p. It is easy to show that this matrix element vanishes. 
After transforming it into the surface integral over r pF 
we get 



e ±ik pF r pF e ik n Ar n A g± i k pn r pn g i kdA r d A 



r P F r n A r^n r dA 

p-s-oo i_ Jpg(a, z) 



(C9) 



Thus we need to estimate a highly oscillatory integral: 

1 



Evidently that this integral and, hence, (P, k^) 

vanishes at p — > oo, whether a stationary phase point 
does exist or not, because the integration brings p to the 
denominator. 

Now we proceed to Mg^f c (P, kdA)- We rewrite is as 
Mg%f C (P, k dA ) = J dr nA J dr pF j$'(r pF ) 

r nA <R n A 

x T^ K ->*(r nA ) [^ pF - -? pF ] *? DCC ( + \r pF , r nA ) 
+ J dr nA J drpFXp'J* (r pF ) 

r nA <R n A 

x T^ K -)*(r nA ) fi nA - ^ nA ] *f DCC ^(v pF , v nA ). 

(Gil) 



J dr nA J ^ pF x pF ] *{v pF )T { ^ ){ ] *{r n/ 



r nA <R nA 

j p F — j p F 
i 



'{r pF , r nA ) 



Ji - lim [ dz [ dasin 2 a cos 2 a e lp9iz ' a) . (CIO) 2/z pF tf P ™°o 

p^°° J J 



lim R 2 F I dfi rpF 



J dr„^T^*(r n A) 



r nA <R„ 



X pF *{*pf) 



CDCC{+) 



(r pF , r nA ) 



dXpF*( v pF) 



dr 



pF 



, T ,cncc(+)/ \ 



dr 



pF 



r vF =R„ 



(C12) 



The matrix element containing n — term of the CDCC 
wave function vanishes because in the subspace r nA < 
R n A at TpF — > oo the deuteron bound state wave func- 
tion exponentially fades away. The terms of the CDCC 
wave function with n > 1 also produce vanishing matrix 
element because the CDCC wave function corresponding 
to these terms in the subspace r nA < R nA at t p f — > oo 
decays as l/r pF , that is the matrix element (|C12|) van- 
ishes as lim R 2 p /RIf — > 0. Thus we arrive at 

R pF ^oa p * vt 



M 



CDCC 



(P k 



dA 



Mg°? C (P, k 



dA 



dr r , 



dr pF T^->(r nA ) fi nA - ~f nA ] ^ DCC{+ \v pF , v nA ) 



r n A<R n A 



— — R 2 nA / drpFXpF ( y pf) *i ; (r P F, r nA )-— T,V (r nA ) - (r nA ) — - 



(r pF , v nA ) 
(C13) 



[1] G. H. Rawitscher, Phys. Rev. C 9, 2010 (1974). 

[2] M. Kamimura, M. Yahiro, Y. Iseri, H. Kameyama, Y. 

Sakuragi, and M. Kawai, Prog. Theor. Phys. Suppl. 89, 

1 (1986). 

[3] N. Austern, Y. Iseri, M. Kamimura, M. Kawai, G. Raw- 
itscher, and M. Yahiro, Phys. Rep. 154, 125 (1987). 

[4] R. C. Johnson and P. J. R. Soper, Phys. Rev. C 1, 055807 
(1970). 

[5] I. J. Thompson, Comput. Phys. Rep. 7, 167 (1988). 
[6] R. Huby and J. R. Mines, Roy. Mod. Phys. 37, 406 
(1965). 

[7] F. S. Levin, Ann. of Phys. 46, 41 (1968); Nucl. Phys. 



A115, 449 (1968). 
[8] H. T. Fortune and C. M. Vincent, Phys. Roy. 185, 1401 

(1969) . 

[9] C. M. Vincent and H. T. Fortune, Phys. Roy. C 2, 782 

(1970) . 

[10] J. Bang and J. Zim nyi, Nucl. Phys. A139, 534 (1969). 

[11] V. E. Bunakov, Nucl. Phys. AI40, 241 (1970). 

[12] V. E. Bunakov, K. A. Gridnev and L. V. Krasnov, Phys. 

Lett. 32B, 587 (1970). 
[13] V. E. Bunakov, K. A. Gridnev and L. V. Krasnov, Phys. 

Lett. 34B, 27 (1971). 
[14] R. Huby, Phys. Lett. 33B, 323 (1970). 



37 



[15] B. J. Cole, R. Huby and J. R. Mines, Phys. Lett. 33B, 
320 (1970). 

[16] B. J. Cole, R. Huby and J. R. Mines, Phys. Rev. Lett. 

26, 264 (1971). 
[17] R. Lipperheide, Phys. Lett. 32B, 555 (1970). 
[18] J. V. Noble, Phys. Rev. C 1, 1742 (1970). 
[19] E. I. Dolinsky, P. O. Dzhamalov and A. M. 

Mukhamedzhanov, Nucl. Phys. A202, 97 (1973). 
[20] Ya. B. Zeldovich, ZhETF (USSR) 39, 776 (1960). 
[21] T. Berggren, Nucl. Phys. A109, 265 (1968). 
[22] A. M. Mukhamedzhanov and A. S. Kadyrov, Phys. Rev. 

C 82, 051601(R) (2010). 
[23] L. D. Blokhintsev, I. Borbely, and E. I. Dolinskii, Fiz. 

Elem. Chastits At. Yadra 8, 1189 (1977) [Sov. J. Part. 

Nucl. 8, 485 (1977)]. 
[24] E. I. Dolinsky and A. M. Mukhamedzhanov, Izv. AN. 

SSSR, Ser. Fiz. 41, 2055 (1977) [Bull. Acad. Sci. USSR, 

Phys. Ser. 41, 55 (1977)]. 
[25] A. M. Mukhamedzhanov and R. E. Tribble, Phys. Rev. 

C 59, 3418 (1999). 
[26] S.A. Goncharov, J. Dobesh, E.I. Dolinskii, 

A.M. Mukhamedzhanov, J. Cejpek, Yad. Fiz. 35, 

662 (1982) [Sov. J. Nucl. Phys. 35 (1982) 383]. 
[27] S. A. Goncharov, A. M. Mukhamedzhanov, E. A. Ro- 

manovskii, G. E. Valiev, I. R. Gulamov, T. Iskhakov, G. 

Nie, N. K. Timofeyuk, R. Yarmukhamedov, V. Kroha, 

and V. A. Stepanenko, Czechoslovak J. of Phys., 37, 168 

(1987). 

[28] I. R. Gulamov, A. M. Mukhamedzhanov, G. K. Nie, Yad. 
Fiz. 58, 1789 (1995) [Phys. At. Nucl. 58, 1689 (1995)]. 

[29] C. A. Gagliardi, R. E. Tribble, A. Azhari, H. L. Clark, 
Y.-W. Lui, A. M. Mukhamedzhanov, A. Sattarov, and L. 
Trache, V. Burjan, J. Cejpek, V. Kroha, . Piskor, and J. 
Vincour, Phys. Rev. C 59, 1149 (1999). 

[30] A. M. Mukhamedzhanov, P. Bern, B. A. Brown, V. Bur- 
jan, C. A. Gagliardi, V. Kroha, J. Novak, F. M. Nunes, 



S. Piskof, F. Pirlepesov, E. Simeckova, R. E. Tribble, 1 
and J. Vincour, Phys. Rev. C 67, 065804 (2003). 

[31] A. M. Mukhamedzhanov, F. M. Nunes, and P. Mohr, 
Phys. Rev. C 77, 051601 (R) (2008). 

[32] A. S. Kadyrov, I. Bray, A. M. Mukhamedzhanov and A. 
T. Stelbovics, Ann. Phys. 324, 1516 (2009). 

[33] M. La Cognata, A. M. Mukhamedzhanov, C. Spitaleri, I. 
Indelicato, M. Aliotta, V. Burjan, S. Cherubini, A. Coc, 
M. Gulino, Z. Hons, G. G. Kiss, V. Kroha, L. Lamia, J. 
Mrazek, S. Palmerini, S. Piskof, R. G. Pizzone, S. M. R. 
Puglia , G. G. Rapisarda, S. Romano, M. L. Sergi, and 
A. Tumino, The Astrophysical Journal Letters 739, L54 
(2011), doi:10.1088/2041-8205/739/2/L54. 

[34] I. Brida, Steven C. Pieper, and R. B. Wiringa, 
arXiv:1106.312 [nucl-th] (2011). 

[35] N. Austern and M. Kawai, Progr. Theor. Phys. 80, 694 
(1988). 

[36] T. Sawada and K. Thushima, Progr. Theor. Phys. 76, 
440 (1986). 

[37] C. Mahaux and H. A. Weidenmiiller "Shell-Model Ap- 
proach to Nuclear Reactions" , North-Holland Publishing 
Company, Amsterdam, 1969. 

[38] C. Spitaleri et al, Phys. Rev. C 63, 055801 (2001). 

[39] A. M. Lane and R. G. Thomas, Rev. Mod. Phys., 30, 

257 (1958). 

[40] J. M. Blatt and L. C. Biedenharn, Rev. Mod. Phys. 24, 

258 (1952). 

[41] Note that when considering the wave function in Ap- 
pendix [X] the channel c = b + B was the entry chan- 
nel. But when considering the deuteron stripping reac- 
tion amplitude, this channel is the exit channel while the 
channel c' — n + A becomes the entry channel in the 
resonant subreaction n + A — > b + B. 

[42] It is shown in (37J how to calculate Xr F Mp in the shell- 
model approach. 



